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