Hi Jeff,
What's wrong with the LCC image? :) Are you looking to use it in a GIS?
For example, this mapserver layer config will show LCC (awips211) images
just fine:
LAYER
NAME "east_vis_1km"
STATUS ON
DATA "/home/ldm/data/gis/images/awips211/GOES13_EAST_VIS.png"
TYPE RASTER
PROJECTION
"proj=lcc"
"lat_0=25"
"lat_1=25"
"lat_2=25"
"lon_0=-95"
"a=6371200.0"
"b=6371200.0"
END
END
Perhaps you could use gdalwrap (a part of the gdal rpm package) to
transform your LCC image into 'unprojected'. I use it all the time to
transform these types of images into spherical mercator or 'unprojected'.
The command is something like
gdalwarp -t_srs 'EPSG:4326' -s_srs 'proj=lcc lat_0=25 lat_1=25 lat_2=25
lon_0=-95 a=6371200.0 b=6371200.0' input_lcc.jpg output_unproj.jpg
I have examples of these LCC files ready for GIS here:
http://mesonet.agron.iastate.edu/data/gis/images/awips211/
daryl
On Fri, 23 Sep 2011, Jeff Lake wrote:
Greetings...
as the title says ...
the following script .. will produce the gem file
but will not produce the gini
I need the satellite images in a non-projected format
and I thought the following would do it
If I use any other projection (lcc, nps, mer) it works
img2gd << EOF
proj = ced
grdarea = 9;-114;57.289;-49.385
kxky = 2560;2560
grdnam = irsat
gdfile = east_ir.gem
glevel =0
gvcord =none
gdattim=0827/0331
cpyfil=
maxgrd=5000
imgtim=last
calimg=no
imgfil= $GEMDATA/EAST-CONUS/4km/IR/IR_${1}
r
e
EOF
gd2img << eoff
GDATTIM = last
GDFILE = east_ir.gem
GLEVEL = 0
GVCORD = none
SCALE = 0
GFUNC = irsat
proj = ced
grdarea = 9;-114;57.289;-49.385
kxky = 2560;2560
CPYFIL =
SATFIL = eastir.gini
CALINFO = 99/10/SAT,0,94,0,286
WMOHDR = TICZ99/CHIZ/
clear = y
r
e
eoff
--
/**
* Daryl Herzmann
* Assistant Scientist -- Iowa Environmental Mesonet
* http://mesonet.agron.iastate.edu
*/
