Mixing
How do I merge RasterBrick columns and rename after extract from NetCDF data ? [R language]
everyone I am just a beginner of R,so following code or question might misunderstand R function. Hopes you guys don't mind it.
I got hundreds ERA Interim Daily Precipitation data from ECMWF.
Each monthly NetCDF file contain twice per day precipitation and extract them by Longitude & Latitude.
Finally convert to CSV file,but itself store (12:00 ;00:00) two columns for
same days rainfall.
Image of output CSV
Ideally ,I would like to obtain single day rainfall for whole year or month. But, don't know how to do it in R.
Q1: Is there any library or function support merge columns names and keep
same dates?Q2: However, I am not sure lubridate library can accomplish or others will be
better.
-----Thanks for any answers---------
If somebody want to view NetCDF file and .csv file ,you can get from my google driver:
https://drive.google.com/drive/folders/1ogQp_-MlvmJCDmzm38dXHbUmWAM6hi1v?usp=sharing
P.S In Google Drive SumMonth.csv is extract result
R Code:
library(raster)
library(ncdf4)
f <- list.files("C:/Users/asus/Downloads/extract",
pattern = "*.nc", full.names = TRUE)
ncloop <-lapply(f,nc_open)# Read whold folder nc.file
brick <-lapply(f,brick,varname = "tp")
brick2 <-stack(brick)# Combine all month
dim(brick2)#check dim right
##### input location file with Lon&Lat
pointCoordinates <- read.table("C:/Users/asus/Downloads/pointFile.csv",
header = TRUE, sep = ",",stringsAsFactors=FALSE)
coordinates(pointCoordinates) <- ~Long + Lat # Assign XY for extract
rasValue=extract(brick2, pointCoordinates)#extract
combinePointValue=cbind(pointCoordinates,rasValue)#construct database
## Outout as csv
cc <- as.data.frame(combinePointValue)
write.csv(cc,file ="SumMonth.csv")
r extract netcdf lubridate r-raster
edited Jan 2 at 21:25
lbusett
3,62921434
asked Jan 2 at 12:37
tannertanner
63
add a comment |
everyone I am just a beginner of R,so following code or question might misunderstand R function. Hopes you guys don't mind it.
I got hundreds ERA Interim Daily Precipitation data from ECMWF.
Each monthly NetCDF file contain twice per day precipitation and extract them by Longitude & Latitude.
Finally convert to CSV file,but itself store (12:00 ;00:00) two columns for
same days rainfall.
Image of output CSV
Ideally ,I would like to obtain single day rainfall for whole year or month. But, don't know how to do it in R.
Q1: Is there any library or function support merge columns names and keep
same dates?Q2: However, I am not sure lubridate library can accomplish or others will be
better.
-----Thanks for any answers---------
If somebody want to view NetCDF file and .csv file ,you can get from my google driver:
https://drive.google.com/drive/folders/1ogQp_-MlvmJCDmzm38dXHbUmWAM6hi1v?usp=sharing
P.S In Google Drive SumMonth.csv is extract result
R Code:
library(raster)
library(ncdf4)
f <- list.files("C:/Users/asus/Downloads/extract",
pattern = "*.nc", full.names = TRUE)
ncloop <-lapply(f,nc_open)# Read whold folder nc.file
brick <-lapply(f,brick,varname = "tp")
brick2 <-stack(brick)# Combine all month
dim(brick2)#check dim right
##### input location file with Lon&Lat
pointCoordinates <- read.table("C:/Users/asus/Downloads/pointFile.csv",
header = TRUE, sep = ",",stringsAsFactors=FALSE)
coordinates(pointCoordinates) <- ~Long + Lat # Assign XY for extract
rasValue=extract(brick2, pointCoordinates)#extract
combinePointValue=cbind(pointCoordinates,rasValue)#construct database
## Outout as csv
cc <- as.data.frame(combinePointValue)
write.csv(cc,file ="SumMonth.csv")
r extract netcdf lubridate r-raster
edited Jan 2 at 21:25
lbusett
3,62921434
asked Jan 2 at 12:37
tannertanner
63
2
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24
add a comment |
everyone I am just a beginner of R,so following code or question might misunderstand R function. Hopes you guys don't mind it.
I got hundreds ERA Interim Daily Precipitation data from ECMWF.
Each monthly NetCDF file contain twice per day precipitation and extract them by Longitude & Latitude.
Finally convert to CSV file,but itself store (12:00 ;00:00) two columns for
same days rainfall.
Image of output CSV
Ideally ,I would like to obtain single day rainfall for whole year or month. But, don't know how to do it in R.
Q1: Is there any library or function support merge columns names and keep
same dates?Q2: However, I am not sure lubridate library can accomplish or others will be
better.
-----Thanks for any answers---------
If somebody want to view NetCDF file and .csv file ,you can get from my google driver:
https://drive.google.com/drive/folders/1ogQp_-MlvmJCDmzm38dXHbUmWAM6hi1v?usp=sharing
P.S In Google Drive SumMonth.csv is extract result
R Code:
library(raster)
library(ncdf4)
f <- list.files("C:/Users/asus/Downloads/extract",
pattern = "*.nc", full.names = TRUE)
ncloop <-lapply(f,nc_open)# Read whold folder nc.file
brick <-lapply(f,brick,varname = "tp")
brick2 <-stack(brick)# Combine all month
dim(brick2)#check dim right
##### input location file with Lon&Lat
pointCoordinates <- read.table("C:/Users/asus/Downloads/pointFile.csv",
header = TRUE, sep = ",",stringsAsFactors=FALSE)
coordinates(pointCoordinates) <- ~Long + Lat # Assign XY for extract
rasValue=extract(brick2, pointCoordinates)#extract
combinePointValue=cbind(pointCoordinates,rasValue)#construct database
## Outout as csv
cc <- as.data.frame(combinePointValue)
write.csv(cc,file ="SumMonth.csv")
r extract netcdf lubridate r-raster
edited Jan 2 at 21:25
lbusett
3,62921434
asked Jan 2 at 12:37
tannertanner
63
everyone I am just a beginner of R,so following code or question might misunderstand R function. Hopes you guys don't mind it.
I got hundreds ERA Interim Daily Precipitation data from ECMWF.
Each monthly NetCDF file contain twice per day precipitation and extract them by Longitude & Latitude.
Finally convert to CSV file,but itself store (12:00 ;00:00) two columns for
same days rainfall.
Image of output CSV
Ideally ,I would like to obtain single day rainfall for whole year or month. But, don't know how to do it in R.
Q1: Is there any library or function support merge columns names and keep
same dates?Q2: However, I am not sure lubridate library can accomplish or others will be
better.
-----Thanks for any answers---------
If somebody want to view NetCDF file and .csv file ,you can get from my google driver:
https://drive.google.com/drive/folders/1ogQp_-MlvmJCDmzm38dXHbUmWAM6hi1v?usp=sharing
P.S In Google Drive SumMonth.csv is extract result
R Code:
library(raster)
library(ncdf4)
f <- list.files("C:/Users/asus/Downloads/extract",
pattern = "*.nc", full.names = TRUE)
ncloop <-lapply(f,nc_open)# Read whold folder nc.file
brick <-lapply(f,brick,varname = "tp")
brick2 <-stack(brick)# Combine all month
dim(brick2)#check dim right
##### input location file with Lon&Lat
pointCoordinates <- read.table("C:/Users/asus/Downloads/pointFile.csv",
header = TRUE, sep = ",",stringsAsFactors=FALSE)
coordinates(pointCoordinates) <- ~Long + Lat # Assign XY for extract
rasValue=extract(brick2, pointCoordinates)#extract
combinePointValue=cbind(pointCoordinates,rasValue)#construct database
## Outout as csv
cc <- as.data.frame(combinePointValue)
write.csv(cc,file ="SumMonth.csv")
r extract netcdf lubridate r-raster
r extract netcdf lubridate r-raster
edited Jan 2 at 21:25
lbusett
3,62921434
asked Jan 2 at 12:37
tannertanner
63
edited Jan 2 at 21:25
lbusett
3,62921434
asked Jan 2 at 12:37
tannertanner
63
edited Jan 2 at 21:25
lbusett
3,62921434
edited Jan 2 at 21:25
lbusett
3,62921434
edited Jan 2 at 21:25
lbusett
3,62921434
3,62921434
asked Jan 2 at 12:37
tannertanner
63
asked Jan 2 at 12:37
tannertanner
63
asked Jan 2 at 12:37
tannertanner
63
63
2
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24
add a comment |
2
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24
2
2
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24
add a comment |
0
active
oldest
votes
Your Answer
StackExchange.ifUsing("editor", function () {
StackExchange.using("externalEditor", function () {
StackExchange.using("snippets", function () {
StackExchange.snippets.init();
});
});
}, "code-snippets");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "1"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f54006534%2fhow-do-i-merge-rasterbrick-columns-and-rename-after-extract-from-netcdf-data%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
Thanks for contributing an answer to Stack Overflow!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f54006534%2fhow-do-i-merge-rasterbrick-columns-and-rename-after-extract-from-netcdf-data%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
2
Welcome to SO. Please consider that users are often unwilling of downloading data from unknown sources. To facilitate users helping you, it would be better to share a minimal part of your data within the question. In this case, you could share the results of dput(head(cc[,1:10])) (see stackoverflow.com/questions/5963269/…)
– lbusett
Jan 2 at 21:24