Hi, The gdalinfo in my post is the gdalinfo of the subdataset 37. I have also tried gdal_translate but when I tried to use gdalwarp using the gdal_translate output, it has error. I really need to use gdalwarp so that I could project the image into epsg:4326(+prj=longlat). If I'm going to use the gdal_translate I am ask to indicate the gcp which is not good in my case because I have many hdf files in my local directory and maybe they have different gcp and I need to batch process them later using the gdalwarp with the +proj=longlat command. What could be the solution?
On Fri, Sep 12, 2014 at 3:00 AM, <[email protected]> wrote: > Send gdal-dev mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > http://lists.osgeo.org/mailman/listinfo/gdal-dev > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of gdal-dev digest..." > > > Today's Topics: > > 1. Re: S57 and area reconstruction (Antti Castr?n) > 2. Re: Extract subdataset from netCDF in a given directory > (Andre Joost) > 3. Re: ERROR 1: Too many points (440 out of 441) failed to > transform, unable to compute output bounds. (Andre Joost) > 4. Re: Extract subdataset from netCDF in a given directory > (Andre Joost) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Thu, 11 Sep 2014 14:10:42 +0300 > From: Antti Castr?n <[email protected]> > To: vwf <[email protected]> > Cc: gdal dev <[email protected]> > Subject: Re: [gdal-dev] S57 and area reconstruction > Message-ID: > < > capcvoo+scqr+cco8jbyss5snaagwlgga22gx+kmk-xd0bpk...@mail.gmail.com> > Content-Type: text/plain; charset=UTF-8 > > Hello Frits, > > It is possible to (re)merge the objects across cell boundaries if they > share common Feature Object IDentifier (FOID), which is a combination > of AGEN, FIND and FINS subfields. Some ECDISes can hide the boundary > depending on masking. Usage of FOIDs vary, so it might be possible for > your datasets - if you are lucky. > > See Clause 3.1 of Edition 2.0 (November 2000) of the ENC Product > Specification (S-57 Appendix B.1) and ENC Encoding Bulletin Number 33. > > Cheers, > > Antti > > > 2014-09-11 8:33 GMT+03:00 vwf <[email protected]>: > > Hello, > > > > My question is about S57, not gdal, but one of you may be able to help > > me. > > > > When an area (like "prcare") is at the edge of a cell (map) there will > > be a vertex on the edge of the cel. With several cells the full area can > > have several lines going through the area. Is there a trick/technique to > > reconstruct the outer edge or remove the extra lines? > > > > Thank you, > > Frits > > > > _______________________________________________ > > gdal-dev mailing list > > [email protected] > > http://lists.osgeo.org/mailman/listinfo/gdal-dev > > > > -- > Antti Castr?n > +358505997538 > [email protected] > > > ------------------------------ > > Message: 2 > Date: Thu, 11 Sep 2014 16:39:11 +0200 > From: Andre Joost <[email protected]> > To: [email protected] > Subject: Re: [gdal-dev] Extract subdataset from netCDF in a given > directory > Message-ID: <[email protected]> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Am 10.09.2014 20:02, schrieb VictoriaH: > > Hello, I am using gdal_translate to extract a subdataset from a netCDF > file > > to a geotiff. If I run the command from the same directory where the > netCDF > > file is located, there is no problem: > > > > Filename: PM25_2001_02_average.nc > > > > Ex. gdal_translate -of GTiff HDF5:"PM25_2001_02_average.nc"://Data > > /var/www/pm25/netCDF/output.tif' > > > > However, if I try to run gdal-translate and point to a netCDF file in a > > different directory, I get an error: > > > > Ex. gdal_translate -of GTiff > > /netCDF/files/HDF5:"PM25_2001_02_average.nc"://Data > > /tiff/files/pm25/netCDF/output.tif' > > > > I am not sure about your syntax at ://Data. This worked for me: > > gdal_translate -a_srs EPSG:4326 -a_ullr -180 -90 180 90 -of "Gtiff" > HDF4_SDS:UNKOWN:"LISOTD_HRFC_V2.3.2013.hdf":0 test.tif > > with all files in the same directory. You have to run gdalinfo first to > get the exact name of the subdataset including the number after the : > The output filename needs to be separated by a blank. > > HTH > Andr? Joost > > > > > ------------------------------ > > Message: 3 > Date: Thu, 11 Sep 2014 16:46:05 +0200 > From: Andre Joost <[email protected]> > To: [email protected] > Subject: Re: [gdal-dev] ERROR 1: Too many points (440 out of 441) > failed to transform, unable to compute output bounds. > Message-ID: <[email protected]> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Am 11.09.2014 10:45, schrieb Love: > > I know I'm not the first to post this topic and encountered this problem. > > But I've search and read some posts related to this and tried their > > solutions but no luck. I really can't figure out what is causing this > > problem. How am I going to solve it? Badly need help! > > > > I have tried gdalwarp -wo SOURCE_EXTRA=100 -wo SAMPLE_GRID=YES and still > no > > luck. > > > > Here are the info: > > > > gdalwarp -t_srs "+proj=longlat +datum=WGS84" > > HDF4_SDS:UNKNOWN:"A2014042042000.L2_LAC.SeAHABS.hdf":37 hdf.tif > > maybe source crs is missing. > > > > *GDALINFO:* > > > > Driver: HDF4Image/HDF4 Dataset > > Files: A2014042042000.L2_LAC.SeAHABS.hdf > > Size is 563, 880 > > Coordinate System is `' > > That is for the whole dataset. Try gdalinfo on the subdataset 37 you want. > > I was successful with extracting the subdatset using gdal_translate, and > work on with that: > > gdal_translate -a_srs EPSG:4326 -a_ullr -180 -90 180 90 -of "Gtiff" > HDF4_SDS:UNKOWN:"LISOTD_HRFC_V2.3.2013.hdf":0 test.tif > > > In your case, a_ullr might be different. > > > HTH, > Andr? Joost > > > > ------------------------------ > > Message: 4 > Date: Thu, 11 Sep 2014 16:58:17 +0200 > From: Andre Joost <[email protected]> > To: [email protected] > Subject: Re: [gdal-dev] Extract subdataset from netCDF in a given > directory > Message-ID: <[email protected]> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > Am 10.09.2014 20:02, schrieb VictoriaH: > > Hello, I am using gdal_translate to extract a subdataset from a netCDF > file > > to a geotiff. If I run the command from the same directory where the > netCDF > > file is located, there is no problem: > > > > Filename: PM25_2001_02_average.nc > > > > Ex. gdal_translate -of GTiff HDF5:"PM25_2001_02_average.nc"://Data > > /var/www/pm25/netCDF/output.tif' > > > > However, if I try to run gdal-translate and point to a netCDF file in a > > different directory, I get an error: > > > > Ex. gdal_translate -of GTiff > > /netCDF/files/HDF5:"PM25_2001_02_average.nc"://Data > > /tiff/files/pm25/netCDF/output.tif' > > > > According to http://www.gdal.org/frmt_hdf5.html, the syntax should be: > > HDF5:file_name:subdataset > > So it has little use to put the path before the HDF5:. I'm not sure if > GDAL accepts the path before the file_name. Running the command in the > same folder should be a safe thing. The output file might go to another > folder. > > HTH, > Andr? Joost > > > > ------------------------------ > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev > > End of gdal-dev Digest, Vol 124, Issue 17 > ***************************************** >
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
