[gdal-dev] Building a Mercator Image from Un-projected data..Major confusion and help requested.

[Cassanova, Bill]

Well, I sorta tried that:

 

gdalwarp -s_src '+proj=WGS84' -t_srs '+proj=merc'
hiradff_200910121800_f180_weather_radar.tiff foo.tiff

 

And the reply is:

 

Usage: gdalwarp [--help-general] [--formats]

    [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"]

    [-order n] [-tps] [-rpc] [-geoloc] [-et err_threshold]

    [-te xmin ymin xmax ymax] [-tr xres yres] [-ts width height]

    [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]

    [-srcnodata "value [value...]"] [-dstnodata "value [value...]"]
-dstalpha

    [-r resampling_method] [-wm memory_in_mb] [-multi] [-q]

    [-cutline datasource] [-cl layer] [-cwhere expression]

    [-csql statement] [-cblend dist_in_pixels]

    [-of format] [-co "NAME=VALUE"]*

    srcfile* dstfile

 

Available resampling methods:

    near (default), bilinear, cubic, cubicspline, lanczos.

 

The file is something GDAL should understand:

 

[vizuser at vizintqa exploded_visualizer]$ gdalinfo
hiradff_200910121800_f180_weather_radar.tiff

Driver: GTiff/GeoTIFF

Files: hiradff_200910121800_f180_weather_radar.tiff

Size is 5796, 2464

Coordinate System is:

GEOGCS["Coordinate System imported from GRIB file",

    DATUM["unknown",

        SPHEROID["unnamed",6367470,0]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433]]

Origin = (234.010419999999982,49.706487000000003)

Pixel Size = (0.010399000000000,-0.010399000000000)

Metadata:

  AREA_OR_POINT=Area

Image Structure Metadata:

  INTERLEAVE=PIXEL

Corner Coordinates:

Upper Left  (     234.010,      49.706) (234d 0'37.51"E, 49d42'23.35"N)

Lower Left  (     234.010,      24.083) (234d 0'37.51"E, 24d 5'0.06"N)

Upper Right (     294.283,      49.706) (294d16'58.89"E, 49d42'23.35"N)

Lower Right (     294.283,      24.083) (294d16'58.89"E, 24d 5'0.06"N)

Center      (     264.147,      36.895) (264d 8'48.20"E, 36d53'41.71"N)

Band 1 Block=5796x1 Type=Byte, ColorInterp=Red

  Mask Flags: PER_DATASET ALPHA

Band 2 Block=5796x1 Type=Byte, ColorInterp=Green

  Mask Flags: PER_DATASET ALPHA

Band 3 Block=5796x1 Type=Byte, ColorInterp=Blue

  Mask Flags: PER_DATASET ALPHA

Band 4 Block=5796x1 Type=Byte, ColorInterp=Alpha

 

 

Am I missing something here?  I would however prefer to do it with code
so that I can go directly from my data to an image without the
intermediate step

Of using gdalwarp.

 

-----Original Message-----
From: gdal-dev-bounces at lists.osgeo.org
[mailto:gdal-dev-bounces at lists.osgeo.org] On Behalf Of Christopher
Barker
Sent: Monday, November 02, 2009 4:39 PM
To: gdal-dev at lists.osgeo.org
Subject: Re: [gdal-dev] Building a Mercator Image from Un-projected
data..Major confusion and help requested.

 

Cassanova, Bill wrote:

> I am looking to gain a better understanding of how I can use GDAL to 

> build a Mercator Image from un-projected data.

 

Do you need to write your own code for this?, or could you just use the 

  gdalwarp utility? (maybe after converting to a format gdal knows).

 

If you do need your own code, you might look at the gdalwarp code.

 

Which brings up a question I could answer by reading the source, but 

since I'm here:

 

Is the core warping functionality in the GDAL API (and wrapped for 

Python?), or in the gdalwarp code -- I may need to do some warping in 

code as well.

 

-Chris

 

 

 

-- 

Christopher Barker, Ph.D.

Oceanographer

 

Emergency Response Division

NOAA/NOS/OR&R            (206) 526-6959   voice

7600 Sand Point Way NE   (206) 526-6329   fax

Seattle, WA  98115       (206) 526-6317   main reception

 

Chris.Barker at noaa.gov

_______________________________________________

gdal-dev mailing list

gdal-dev at lists.osgeo.org

http://lists.osgeo.org/mailman/listinfo/gdal-dev

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20091102/2bfcb0bc/attachment-0001.html