Thanks for the input Kyle. Just a quick question. In the second operation using gdal_translate, you indicate an -s_srs of EPSG:4326 for the *.vrt file. Wasn't the vrt file written in the Lambert projection in the first step? And if so, the second step could probably be done without having to specify either s_srs or t_srs?
Is this true? Thanks again, Bill On Tue, 2010-03-02 at 11:36 -0700, Kyle Shannon wrote: > Bill, > > It may be necessary to create an intermediate file that supports > creation (I think the (rw+) flag on gdalinfo --formats shows which > formats support create, AAIGrid supports CreateCopy() but not > Create()). Use gdalwarp to make a vrt file, then translate the file > to AAIGrid: > > /usr/local/gdal/bin/gdalwarp -of VRT -s_srs EPSG:4326 -t_srs > '+proj=lcc +lat_1=33n +lat_2=45n +lon_0=97w' mcd12_sep2109.asc > lambert_mcd12_sep2109.vrt > > then: > > /usr/local/gdal/bin/gdal_translate -of AAIGrid -s_srs EPSG:4326 -t_srs > '+proj=lcc +lat_1=33n +lat_2=45n +lon_0=97w' lambert_mcd12_sep2109.vrt > lambert_mcd12_sep2109.asc > > hope that helps, > > regards, > > -Kyle > > > On Tue, Mar 2, 2010 at 8:57 AM, Bill Hudspeth <w...@unm.edu> wrote: > Hello, > > Thanks for responding. When I use gdalinfo --formats, I see > that AAIGrid > is included in the list of acceptable formats. However, when I > use > > /usr/local/gdal/bin/gdalwarp -of AAIGrid -s_srs EPSG:4326 > -t_srs > '+proj=lcc +lat_1=33n +lat_2=45n +lon_0=97w' mcd12_sep2109.asc > lambert_mcd12_sep2109.asc > > I get, > > Output driver `AAIGrid' not recognised or does not support > direct output file creation. The following format drivers are > configured > and support direct output: > VRT: Virtual Raster > GTiff: GeoTIFF > NITF: National Imagery Transmission Format > HFA: Erdas Imagine Images (.img) > ELAS: ELAS > MEM: In Memory Raster > BMP: MS Windows Device Independent Bitmap > PCIDSK: PCIDSK Database File > ILWIS: ILWIS Raster Map > HDF4Image: HDF4 Dataset > PNM: Portable Pixmap Format (netpbm) > ENVI: ENVI .hdr Labelled > EHdr: ESRI .hdr Labelled > PAux: PCI .aux Labelled > MFF: Atlantis MFF Raster > MFF2: Atlantis MFF2 (HKV) Raster > BT: VTP .bt (Binary Terrain) 1.3 Format > IDA: Image Data and Analysis > RMF: Raster Matrix Format > > > ???? > > > > > > > On Tue, 2010-03-02 at 08:50 +0100, Markus Neteler wrote: > > On Mon, Mar 1, 2010 at 7:12 PM, Bill Hudspeth <w...@unm.edu> > wrote: > > ... > > > When I try to re-project the exported AAIGrid to another > projection, I > > > use: > > > > > > gdalwarp -s_srs EPSG:4326 -t_srs '+proj=lcc +lat__1=33n > +lat_2=45n > > > +lon_0=97w' MODIS_output_latlong.asc > MODIS_output_lambert.asc > > > > Note that without -of format] you generate a GeoTIFF file... > > > > > While the resultant projection information looks correct, > > > > > > Driver: GTiff/GeoTIFF > > ^^^^^^^^^^^^^^^^^^^^^^^ > > > > > The resultant text file has only the following on the top > line: > > > > ... it is binary, not a text file. check > > > > gdalwarp --formats > > > > for the available formats. > > > > Markus > > > _______________________________________________ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev > > > _______________________________________________ > gdal-dev mailing list > gdal-dev@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
