ฉันต้องการโหลดไฟล์. adm ลงในอาร์ข้อมูลมาจากหน้านี้: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057
ฉันลองใช้รหัสต่อไปนี้ที่ฉันพบหลังจากการวิจัยในอินเทอร์เน็ต ปัญหาคือว่าในคลาส RasterLayer ฉันได้รับค่า negativ ที่ไม่ควรมี ฉันไม่รู้ว่าทำไมสิ่งนี้จึงเกิดขึ้นดังนั้นหวังว่าใครบางคนจะสามารถช่วยฉันได้!
รหัส:
library(rgdal)
library(RColorBrewer)
dpath<- path...
x <- new("GDALReadOnlyDataset", dpath)
getDriver(x)
getDriverLongName(getDriver(x))
xx<-asSGDF_GROD(x)
r <- raster(xx)ผลลัพธ์สำหรับ 'r' คือ:
คลาส r: RasterLayer ขนาด: 2160, 4320, 9331200 (nrow, ncol, ncell) ความละเอียด: 0.08333333, 0.08333333 (x, y) ขอบเขต: -180, 180, -90, 90 (xmin, xmax, ymin, ymax) coord อ้าง : + proj = longlat + ellps = WGS84 + towgs84 = 0,0,0,0,0,0,0 + no_defs แหล่งข้อมูล: ในชื่อหน่วยความจำ: ค่าชื่อ band1: -997, 16 (นาที, สูงสุด)
ค่า '16' ในค่าอ้างอิงถึง 16 คลาสของความยาวของระยะเวลาการเติบโต แต่ฉันสงสัยว่า '-997' เหล่านั้นมาจากไหน อาจมีบางอย่างผิดปกติกับ coord โทษ?
นี่คือข้อมูลสรุปของ 'xx':
สรุปข้อมูล: ขั้นต่ำ วันที่ 1 ค่าเฉลี่ยมัธยฐาน 3 Qu แม็กซ์ NA's -997 3 5 -9 8 16 7123158
และถ้าเราดูข้อมูลใน xx ให้ละเอียดยิ่งขึ้น:
ตาราง (xx $ band1)
-997 1 2 3 4 5 6 7 8 9 10 11 12 31711 429643 83011 166674 207228 270161 240958 183342 118608 98795 88473 73743 56022 13 14 15 16 30104 45521 52216 31832
มีเพียง '-997' สิ่งนี้อยู่ในนั้น ฉันคิดว่า NA เป็นมหาสมุทรดังนั้นมีอะไรผิดปกติในการโหลดข้อมูลหรือฉันแค่ไม่เข้าใจข้อมูลหรือไม่