Re: [gembud] projection dilemma

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
 */