<div dir="ltr">Hello, <div><br></div><div>I have a question about an error I am not receiving when I believe I should.</div><div><br></div><div>Below are three datasets that with their SRID</div><div><div style="font-size:13px">igbp_2012 (6842), us_dem_prj (3410), and us_states (4326).</div><div style="font-size:13px"><br></div><div style="font-size:13px">When I run a query determining the SRID for these datasets in PostgreSQL they confirm their SRIDs</div><div style="font-size:13px">select st_srid(r.rast), st_srid(p.geom) from dasymetric.us_dem_prj r, dasymetric.us_states p limit 5<br></div><div style="font-size:13px"><br></div><div style="font-size:13px">When I run the next query, I do not receive an error, but no records are returned. The raster is in 3410 and the states geom column is in 4326</div><div style="font-size:13px"><div>SELECT state_label, ST_Clip(r.rast,p.geom) as rast FROM dasymetric.us_states p inner join dasymetric.us_dem_prj r on ST_Intersects(p.geom, r.rast) WHERE p.state_label = 'Minnesota'</div></div><div style="font-size:13px"><br></div><div style="font-size:13px">However, If I adjust the raster dataset to IGBP_2012, the error</div><div style="font-size:13px"><div>ERROR: Raster and geometry do not have the same SRID, shows up.</div><div><br></div><div>Why am I not receiving a conflict projection error for a raster in 3410 and a geometry in 4326? Any idea why that is occurring?</div></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">David Haynes, Ph.D.<div>Research Associate Terra Populus</div><div>Minnesota Population Center</div><div><a href="http://www.terrapop.org" target="_blank">www.terrapop.org</a></div></div></div>
</div></div>