<div dir="ltr">Paul,<div><br></div><div>I suppose POINTM is the solution to get the data somehow into PostGIS but I will have a handful of options to process them and millions of rows with points,</div><div>PostGIS has raster data type tailored to 2-d arrays (please, correct me if I am wrong):</div><div><a href="https://postgis.net/docs/RT_reference.html">https://postgis.net/docs/RT_reference.html</a><br></div><div><a href="http://postgis.net/docs/manual-2.2/using_raster_dataman.html">http://postgis.net/docs/manual-2.2/using_raster_dataman.html</a><br></div><div><br></div><div>Therefore, I seek a way to import a 3-d array and be able to apply the functionality described in the links above </div><div><br></div><div>Antonio</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-10-30 15:25 GMT+03:00 Paul Ramsey <span dir="ltr"><<a href="mailto:pramsey@cleverelephant.ca" target="_blank">pramsey@cleverelephant.ca</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">As others have noted, a POINTM or MULTIPOINTM will serve to store your data just fine, but what you plan to *do* with that data after will determine whether a relational database is really the correct tool for you.<div><br></div><div>ATB,</div><div>P</div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 30, 2017 at 4:56 AM, Antonio Rodriges <span dir="ltr"><<a href="mailto:antonio.rrz@gmail.com" target="_blank">antonio.rrz@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I suppose it has a bit different purpose (like an efficient handling of sparsity which is not the case with dense climate data) and I hope there is an easier solution (I just have 3 dimensions...)<span class="m_113721193486212023HOEnZb"><font color="#888888"><div><br></div><div>Antonio</div></font></span></div><div class="m_113721193486212023HOEnZb"><div class="m_113721193486212023h5"><div class="gmail_extra"><br><div class="gmail_quote">2017-10-30 14:52 GMT+03:00 Stephen V. Mather <span dir="ltr"><<a href="mailto:svm@clevelandmetroparks.com" target="_blank">svm@clevelandmetroparks.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="auto">
Ya, I’m not sure point clouds are at all the fix. They just address the dimensionality question well, though not the gridded data requirement.
<div><br>
</div>
<div>Cheers,</div>
<div>Best,</div>
<div><span>Steve <br>
<br>
<div id="m_113721193486212023m_4804648563658017774m_7238205789501690193AppleMailSignature">
<div><br>
</div>
<img alt="http://sig.cmparks.net/cmp-ms-90x122.png" align="left" height="122" hspace="12" width="90" src="http://sig.cmparks.net/cmp-ms-90x122.png"><span style="background-color:rgba(255,255,255,0)"><strong>Stephen V. Mather</strong><br>
GIS Manager<br>
<a dir="ltr" href="tel:(216)%20635-3243" target="_blank">(216) 635-3243</a> (Work)</span>
<div><a href="tel:(216)%20339-6347" value="+12163396347" target="_blank">(216) 339-6347</a> (Cell)<br>
<div><span style="background-color:rgba(255,255,255,0)">--sent from phone--</span></div>
<div><span style="background-color:rgba(255,255,255,0)"><br>
</span></div>
<div><span style="background-color:rgba(255,255,255,0)"><br>
</span></div>
</div>
</div>
</span><div><div class="m_113721193486212023m_4804648563658017774h5"><div><br>
On Oct 30, 2017, at 07:49, Antonio Rodriges <<a href="mailto:antonio.rrz@gmail.com" target="_blank">antonio.rrz@gmail.com</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div>
<div dir="ltr">Thank you for pointing to this tool.
<div><br>
</div>
<div>However, I thought that since PostGIS uses GDAL it may be easier to import such arrays, e.g. just split them onto individual 2-d grids (since PostGIS mainly understands 2-d grids).</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2017-10-30 14:34 GMT+03:00 Stephen V. Mather <span dir="ltr">
<<a href="mailto:svm@clevelandmetroparks.com" target="_blank">svm@clevelandmetroparks.com</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="auto">I don’t know if it’s the ideal tool for the job, as it’s more flexible than you need, not being a regularized grid but a point cloud, but you might look to the pgPointCloud extension: <a href="https://github.com/pgpointcloud/pointcloud" target="_blank">https://github.com/<wbr>pgpointcloud/pointcloud</a>
<div><br>
</div>
<div>Cheers,</div>
<div>Best,</div>
<div>Steve<br>
<br>
<div id="m_113721193486212023m_4804648563658017774m_7238205789501690193m_-5296256961561710695AppleMailSignature">
<div><br>
</div>
<img alt="http://sig.cmparks.net/cmp-ms-90x122.png" align="left" height="122" hspace="12" width="90" src="http://sig.cmparks.net/cmp-ms-90x122.png"><span style="background-color:rgba(255,255,255,0)"><strong>Stephen V. Mather</strong><br>
GIS Manager<br>
<a dir="ltr" href="tel:(216)%20635-3243" target="_blank">(216) 635-3243</a> (Work)</span>
<div><a href="tel:(216)%20339-6347" value="+12163396347" target="_blank">(216) 339-6347</a> (Cell)<br>
<div><span style="background-color:rgba(255,255,255,0)">--sent from phone--</span></div>
<div><span style="background-color:rgba(255,255,255,0)"><br>
</span></div>
<div><span style="background-color:rgba(255,255,255,0)"><br>
</span></div>
</div>
</div>
<div>
<div class="m_113721193486212023m_4804648563658017774m_7238205789501690193h5">
<div><br>
On Oct 30, 2017, at 07:09, Antonio Rodriges <<a href="mailto:antonio.rrz@gmail.com" target="_blank">antonio.rrz@gmail.com</a>> wrote:<br>
<br>
</div>
<blockquote type="cite">
<div><span>Hello,</span><br>
<span></span><br>
<span>Thank you for the reply, however my data is slightly different. Sorry</span><br>
<span>that I did not make it clearer at the very beginning.</span><br>
<span></span><br>
<span>Actually I would like to import a dense, 3-d array of wind speed (a</span><br>
<span>time series of grids, each grid point contains the wind speed value)</span><br>
<span></span><br>
<span>The array is stored as a NetCDF file</span><br>
<span>FIles are here <a href="https://www.esrl.noaa.gov/psd/cgi-bin/db_search/DBListFiles.pl?did=61&tid=59909&vid=4298" target="_blank">
https://www.esrl.noaa.gov/psd/<wbr>cgi-bin/db_search/DBListFiles.<wbr>pl?did=61&tid=59909&vid=4298</a></span><br>
<span></span><br>
<span>The size of the array and its dimensions are below</span><br>
<span></span><br>
<span>dimensions:</span><br>
<span> lat = 94 ;</span><br>
<span> lon = 192 ;</span><br>
<span> time = 1460;</span><br>
<span></span><br>
<span>The array</span><br>
<span></span><br>
<span>short uwnd(time,lat,lon) ;</span><br>
<span> uwnd:long_name = "6-Hourly Forecast of U-wind at 10 m" ;</span><br>
<span> uwnd:valid_range = -32765s, -8765s ;</span><br>
<span> uwnd:unpacked_valid_range = -120.f, 120.f ;</span><br>
<span> uwnd:actual_range = -38.2f, 38.07f ;</span><br>
<span> uwnd:units = "m/s" ;</span><br>
<span> uwnd:add_offset = 207.65f ;</span><br>
<span> uwnd:scale_factor = 0.01f ;</span><br>
<span> uwnd:missing_value = 32766s ;</span><br>
<span> uwnd:_FillValue = -32767s ;</span><br>
<span> uwnd:precision = 2s ;</span><br>
<span> uwnd:least_significant_di<wbr>git = 1s ;</span><br>
<span> uwnd:GRIB_id = 33s ;</span><br>
<span> uwnd:GRIB_name = "U GRD" ;</span><br>
<span> uwnd:var_desc = "u-wind" ;</span><br>
<span> uwnd:dataset = "NCEP/DOE AMIP-II Reanalysis (Reanalysis-2)" ;</span><br>
<span> uwnd:level_desc = "10 m" ;</span><br>
<span> uwnd:statistic = "Individual Obs" ;</span><br>
<span> uwnd:parent_stat = "Other" ;</span><br>
<span> uwnd:standard_name = "eastward_wind" ;</span><br>
<span></span><br>
<span>2017-10-30 11:04 GMT+03:00 Giuseppe Broccolo <<a href="mailto:g.broccolo.7@gmail.com" target="_blank">g.broccolo.7@gmail.com</a>>:</span><br>
<blockquote type="cite"><span>Hi Antonio,</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>2017-10-29 12:31 GMT+01:00 Antonio Rodriges <<a href="mailto:antonio.rrz@gmail.com" target="_blank">antonio.rrz@gmail.com</a>>:</span><br>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span></span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>Hello,</span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span></span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>Whether PostGIS allow importing 3-d, 4-d, etc. arrays or only 2-d arrays?</span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span></span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>Specifically, I have a 3-d array with axes (time, lat, lon).</span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>Does this mean that I need to split it onto 2-d bands (lat, lon) and</span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>import the number of bands that is equal to the number of time steps</span><br>
</blockquote>
</blockquote>
<blockquote type="cite">
<blockquote type="cite"><span>in the 3-d array?</span><br>
</blockquote>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>Which is the data source from which you import the data (e.g. textual,</span><br>
</blockquote>
<blockquote type="cite"><span>etc.)?</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>If I've correctly understood, you have arrays where geospatial and</span><br>
</blockquote>
<blockquote type="cite"><span>non-geospatial information</span><br>
</blockquote>
<blockquote type="cite"><span>is present, each one providing a "dimension" of the array.</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>Just FYI, in PostGIS is possible to define mixed, structured data with</span><br>
</blockquote>
<blockquote type="cite"><span>constructors like POINTM</span><br>
</blockquote>
<blockquote type="cite"><span>and POINT, that allow to add a further dimension to the 2D/3D (respectively)</span><br>
</blockquote>
<blockquote type="cite"><span>geospatial ones, that</span><br>
</blockquote>
<blockquote type="cite"><span>includes a scalar information.</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>Hope this can help in your import, otherwise provide more information about</span><br>
</blockquote>
<blockquote type="cite"><span>source data and how</span><br>
</blockquote>
<blockquote type="cite"><span>you'd like to import.</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>Giuseppe.</span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span></span><br>
</blockquote>
<blockquote type="cite"><span>______________________________<wbr>_________________</span><br>
</blockquote>
<blockquote type="cite"><span>postgis-users mailing list</span><br>
</blockquote>
<blockquote type="cite"><span><a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a></span><br>
</blockquote>
<blockquote type="cite"><span><a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a></span><br>
</blockquote>
<span>______________________________<wbr>_________________</span><br>
<span>postgis-users mailing list</span><br>
<span><a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a></span><br>
<span><a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a></span></div>
</blockquote>
</div>
</div>
</div>
</div>
<br>
______________________________<wbr>_________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a><br>
</blockquote>
</div>
<br>
</div>
</div>
</blockquote>
<blockquote type="cite">
<div><span>______________________________<wbr>_________________</span><br>
<span>postgis-users mailing list</span><br>
<span><a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a></span><br>
<span><a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a></span></div>
</blockquote>
</div></div></div>
</div>
<br>______________________________<wbr>_________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a><br></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailma<wbr>n/listinfo/postgis-users</a><br></blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/postgis-users" rel="noreferrer" target="_blank">https://lists.osgeo.org/<wbr>mailman/listinfo/postgis-users</a><br></blockquote></div><br></div>