MPA names not matching up with data_CPCAD_areas

[j-harbin]

Hey @remi-daigle ,

I am working on pulling scrapping the MPA websites to obtain their site level objectives. I was specifically working with the Atlantic Ocean MPAs and I noticed that majority of them are being captured from the data_CPCAS_areas function. For example, the names of the MPAs on the DFO website are as follows:

But when I do the following:

MPAs <- data_CPCAD_areas(data_bioregion(),  zones = FALSE)
MPAs$NAME_E

I get the following:

[1] "Boot Island National Wildlife Area"                                          
 [2] "Caughey-Taylor, Private Land Component"                                      
 [3] "Emerald Basin Sponge Conservation Area"                                      
 [4] "Grand Manan Migratory Bird Sanctuary"                                        
 [5] "Gully Marine Protected Area"                                                 
 [6] "John Lusby Marsh National Wildlife Area"                                     
 [7] "Jordan Basin Conservation Area"                                              
 [8] "Kejimkujik National Park And National Historic Site Of Canada"               
 [9] "Machias Seal Island Migratory Bird Sanctuary"                                
[10] "McPhersons Point"                                                            
[11] "Musquash Estuary Marine Protected Area"                                      
[12] "Northeast Channel Coral Conservation Area (Restricted Bottom Fisheries Zone)"
[13] "Port Hebert Migratory Bird Sanctuary"                                        
[14] "Port Joli Migratory Bird Sanctuary"                                          
[15] "Sambro Bank Sponge Conservation Area"                                        
[16] "St. Anns Bank Marine Protected Area"                                         
[17] "Western/Emerald Banks Conservation Area (Restricted Fisheries Zone)"   

You'll notice that places such as banc-des-americains, basin head, east port, gilbert bay, and laurentian channel aren't there. Is there another way that I should be obtaining the MPA names?

[remi-daigle]

Part of the issue is that you looking at a "whole Atlantic" list of MPAs, but only looking on the Scotian Shelf with R. The other part of that is the CPCAD contains more than just Ocean's Act MPAs. To get the same list as the website (i.e. Atlantic Ocean's Act MPAs), you could do:

atlantic <- dplyr::bind_rows(data_bioregion("Scotian Shelf"),
                             data_bioregion("Newfoundland-Labrador Shelves"),
                             data_bioregion("Gulf of Saint Lawrence"))
MPAs <- data_CPCAD_areas(sf::st_make_valid(atlantic),  zones = TRUE)
unique(MPAs$NAME_E[MPAs$MECH_E=="Oceans Act ( 1996, c. 31 )"])

That being said, your approach is correctly returning all protected areas in the Scotian Shelf Bioregion which I think is appropriate for our "network approach"