bcgov / designatedlands

Python script to combine conservation related spatial data from many sources to create a single 'Designated Lands' layer for British Columbia
Apache License 2.0
9 stars 4 forks source link

north edge of 20k tiles does not match official BC Boundary #47

Closed smnorris closed 6 years ago

smnorris commented 6 years ago

The 20k grid/tile layer from DataBC https://catalogue.data.gov.bc.ca/dataset/bcgs-1-20-000-grid has a northern boundary of 60deg N.

The BC Boundary file that defines the provincial border has a northern boundary that is close but not exactly at 60degN, it is generally a few metres farther north.

Since inputs are cut by tile to speed processing, some area is lost. The effect is very minor (especially as there are few if any designated lands in this gap) but the difference is noticeable when calculating total areas.

To fix, nudge the northernmost tile boundaries up by ~250m in this step https://github.com/bcgov/designatedlands/blob/master/sql/create_tiles.sql