<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv=Content-Type content="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 12 (filtered medium)">
<style>
<!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:"OCR A Extended";
panose-1:2 1 5 9 2 1 2 1 3 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:10.5pt;
font-family:"OCR A Extended";}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"OCR A Extended";}
.MsoChpDefault
{mso-style-type:export-only;}
@page Section1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.Section1
{page:Section1;}
-->
</style>
<!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang=EN-US link=blue vlink=purple>
<div class=Section1>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>Well, I sorta tried that:<o:p></o:p></b></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>gdalwarp -s_src '+proj=WGS84' -t_srs '+proj=merc'
hiradff_200910121800_f180_weather_radar.tiff foo.tiff<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>And the reply is:<o:p></o:p></b></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Usage: gdalwarp [--help-general] [--formats]<o:p></o:p></p>
<p class=MsoPlainText> [-s_srs srs_def] [-t_srs srs_def] [-to
"NAME=VALUE"]<o:p></o:p></p>
<p class=MsoPlainText> [-order n] [-tps] [-rpc] [-geoloc]
[-et err_threshold]<o:p></o:p></p>
<p class=MsoPlainText> [-te xmin ymin xmax ymax] [-tr xres
yres] [-ts width height]<o:p></o:p></p>
<p class=MsoPlainText> [-wo "NAME=VALUE"] [-ot
Byte/Int16/...] [-wt Byte/Int16]<o:p></o:p></p>
<p class=MsoPlainText> [-srcnodata "value
[value...]"] [-dstnodata "value [value...]"] -dstalpha<o:p></o:p></p>
<p class=MsoPlainText> [-r resampling_method] [-wm
memory_in_mb] [-multi] [-q]<o:p></o:p></p>
<p class=MsoPlainText> [-cutline datasource] [-cl layer]
[-cwhere expression]<o:p></o:p></p>
<p class=MsoPlainText> [-csql statement] [-cblend
dist_in_pixels]<o:p></o:p></p>
<p class=MsoPlainText> [-of format] [-co
"NAME=VALUE"]*<o:p></o:p></p>
<p class=MsoPlainText> srcfile* dstfile<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Available resampling methods:<o:p></o:p></p>
<p class=MsoPlainText> near (default), bilinear, cubic,
cubicspline, lanczos.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>The file is something GDAL should understand:<o:p></o:p></b></p>
<p class=MsoPlainText><b><o:p> </o:p></b></p>
<p class=MsoPlainText><b>[vizuser@vizintqa exploded_visualizer]$ gdalinfo
hiradff_200910121800_f180_weather_radar.tiff<o:p></o:p></b></p>
<p class=MsoPlainText>Driver: GTiff/GeoTIFF<o:p></o:p></p>
<p class=MsoPlainText>Files: hiradff_200910121800_f180_weather_radar.tiff<o:p></o:p></p>
<p class=MsoPlainText>Size is 5796, 2464<o:p></o:p></p>
<p class=MsoPlainText>Coordinate System is:<o:p></o:p></p>
<p class=MsoPlainText>GEOGCS["Coordinate System imported from GRIB
file",<o:p></o:p></p>
<p class=MsoPlainText> DATUM["unknown",<o:p></o:p></p>
<p class=MsoPlainText>
SPHEROID["unnamed",6367470,0]],<o:p></o:p></p>
<p class=MsoPlainText> PRIMEM["Greenwich",0],<o:p></o:p></p>
<p class=MsoPlainText>
UNIT["degree",0.0174532925199433]]<o:p></o:p></p>
<p class=MsoPlainText>Origin = (234.010419999999982,49.706487000000003)<o:p></o:p></p>
<p class=MsoPlainText>Pixel Size = (0.010399000000000,-0.010399000000000)<o:p></o:p></p>
<p class=MsoPlainText>Metadata:<o:p></o:p></p>
<p class=MsoPlainText> AREA_OR_POINT=Area<o:p></o:p></p>
<p class=MsoPlainText>Image Structure Metadata:<o:p></o:p></p>
<p class=MsoPlainText> INTERLEAVE=PIXEL<o:p></o:p></p>
<p class=MsoPlainText>Corner Coordinates:<o:p></o:p></p>
<p class=MsoPlainText>Upper Left (
234.010, 49.706) (234d 0'37.51"E,
49d42'23.35"N)<o:p></o:p></p>
<p class=MsoPlainText>Lower Left (
234.010, 24.083) (234d 0'37.51"E, 24d
5'0.06"N)<o:p></o:p></p>
<p class=MsoPlainText>Upper Right (
294.283, 49.706) (294d16'58.89"E,
49d42'23.35"N)<o:p></o:p></p>
<p class=MsoPlainText>Lower Right (
294.283, 24.083) (294d16'58.89"E, 24d
5'0.06"N)<o:p></o:p></p>
<p class=MsoPlainText>Center
( 264.147, 36.895) (264d
8'48.20"E, 36d53'41.71"N)<o:p></o:p></p>
<p class=MsoPlainText>Band 1 Block=5796x1 Type=Byte, ColorInterp=Red<o:p></o:p></p>
<p class=MsoPlainText> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p>
<p class=MsoPlainText>Band 2 Block=5796x1 Type=Byte, ColorInterp=Green<o:p></o:p></p>
<p class=MsoPlainText> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p>
<p class=MsoPlainText>Band 3 Block=5796x1 Type=Byte, ColorInterp=Blue<o:p></o:p></p>
<p class=MsoPlainText> Mask Flags: PER_DATASET ALPHA<o:p></o:p></p>
<p class=MsoPlainText>Band 4 Block=5796x1 Type=Byte, ColorInterp=Alpha<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><b>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<o:p></o:p></b></p>
<p class=MsoPlainText><b>Of using gdalwarp.<o:p></o:p></b></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>-----Original Message-----<br>
From: gdal-dev-bounces@lists.osgeo.org
[mailto:gdal-dev-bounces@lists.osgeo.org] On Behalf Of Christopher Barker<br>
Sent: Monday, November 02, 2009 4:39 PM<br>
To: gdal-dev@lists.osgeo.org<br>
Subject: Re: [gdal-dev] Building a Mercator Image from Un-projected data..Major
confusion and help requested.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Cassanova, Bill wrote:<o:p></o:p></p>
<p class=MsoPlainText>> I am looking to gain a better understanding of how I
can use GDAL to <o:p></o:p></p>
<p class=MsoPlainText>> build a Mercator Image from un-projected data.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Do you need to write your own code for this?, or could
you just use the <o:p></o:p></p>
<p class=MsoPlainText> gdalwarp utility? (maybe after converting to a
format gdal knows).<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>If you do need your own code, you might look at the
gdalwarp code.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Which brings up a question I could answer by reading the
source, but <o:p></o:p></p>
<p class=MsoPlainText>since I'm here:<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Is the core warping functionality in the GDAL API (and
wrapped for <o:p></o:p></p>
<p class=MsoPlainText>Python?), or in the gdalwarp code -- I may need to do
some warping in <o:p></o:p></p>
<p class=MsoPlainText>code as well.<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>-Chris<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>-- <o:p></o:p></p>
<p class=MsoPlainText>Christopher Barker, Ph.D.<o:p></o:p></p>
<p class=MsoPlainText>Oceanographer<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Emergency Response Division<o:p></o:p></p>
<p class=MsoPlainText>NOAA/NOS/OR&R
(206) 526-6959 voice<o:p></o:p></p>
<p class=MsoPlainText>7600 Sand Point Way NE (206)
526-6329 fax<o:p></o:p></p>
<p class=MsoPlainText>Seattle, WA
98115 (206) 526-6317 main
reception<o:p></o:p></p>
<p class=MsoPlainText><o:p> </o:p></p>
<p class=MsoPlainText>Chris.Barker@noaa.gov<o:p></o:p></p>
<p class=MsoPlainText>_______________________________________________<o:p></o:p></p>
<p class=MsoPlainText>gdal-dev mailing list<o:p></o:p></p>
<p class=MsoPlainText>gdal-dev@lists.osgeo.org<o:p></o:p></p>
<p class=MsoPlainText>http://lists.osgeo.org/mailman/listinfo/gdal-dev<o:p></o:p></p>
</div>
</body>
</html>