Closed dfilatow closed 4 years ago
Interesting. Seems to work with the BCGW layer name:
bcdc_query_geodata("WHSE_FOREST_VEGETATION.BEC_BIOGEOCLIMATIC_POLY")
Querying 'WHSE_FOREST_VEGETATION.BEC_BIOGEOCLIMATIC_POLY' record
● Using collect() on this object will return 15267 features and 20 fields
● At most six rows of the record are printed here
───────────────────────────────────────────────
Simple feature collection with 6 features and 20 fields
geometry type: POLYGON
dimension: XY
bbox: xmin: 749212.5 ymin: 1244512 xmax: 843512.5 ymax: 1265594
CRS: 3005
# A tibble: 6 x 21
id FEATURE_CLASS_S… ZONE SUBZONE VARIANT
<chr> <int> <chr> <chr> <chr>
1 WHSE… 435 CMA un NA
2 WHSE… 435 CMA un NA
3 WHSE… 435 CMA un NA
4 WHSE… 435 CMA un NA
5 WHSE… 435 CMA un NA
6 WHSE… 435 CMA un NA
# … with 16 more variables: PHASE <chr>,
# NATURAL_DISTURBANCE <chr>,
# MAP_LABEL <chr>, BGC_LABEL <chr>,
# ZONE_NAME <chr>, SUBZONE_NAME <chr>,
# VARIANT_NAME <chr>, PHASE_NAME <chr>,
# NATURAL_DISTURBANCE_NAME <chr>,
# FEATURE_AREA_SQM <dbl>,
# FEATURE_LENGTH_M <dbl>,
# FEATURE_AREA <int>, FEATURE_LENGTH <int>,
# OBJECTID <int>, SE_ANNO_CAD_DATA <chr>,
# geometry <POLYGON [m]>
hmmm. bcdc_query_geodata("WHSE_FOREST_VEGETATION.BEC_BIOGEOCLIMATIC_POLY") works, but these do not... bcdc_query_geodata("bec-map") bcdc_query_geodata("f358a53b-ffde-4830-a325-a5a03ff672c3")
I think it might be because the WMS resource lists "Resource Storage Location" as "Web or FTP Site" rather than "BC Geographic Warehouse". We use this field to help determine if a resource has wms service...
https://catalogue.data.gov.bc.ca/dataset/bec-map/resource/3ec24cb4-f78d-48a9-87e7-516763f00cc8
Edit: I've emailed the catalogue team to find out more
The 3 BEC resources have now been changed to "BC Geographic Warehouse". Thanks for noting this.
I think we can close this now. Thanks for reporting @dfilatow and for the fix @pumapants!
is this the same issue as #107 that was closed?