组保留的标签 余下为需要删除的标签
unset($oldtag[$key]);
}
}
}
if (!empty($oldtag)) {
$tagids = array();
foreach ($oldtag as $tagid => $tagname) {
$tagids[] = $tagid;
}
well_oldtag_delete($tagids, $tid);
}
$r = well_tag_process($tid, $fid, $create_tag, $tagarr);
return $r;
}
// 删除标签和绑定的主题
function well_oldtag_delete($tagids, $tid)
{
$pagesize = count($tagids);
$arrlist = well_tag_find_by_tagids($tagids, 1, $pagesize);
$delete_tagids = array(); // 删除
$tagids = array();
$n = 0;
foreach ($arrlist as $val) {
++$n;
if (1 == $val['count']) {
// 只有一个主题
$delete_tagids[] = $val['tagid'];
} else {
$tagids[] = $val['tagid'];
}
}
!empty($delete_tagids) and well_tag_delete($delete_tagids);
$arlist = well_tag_thread_find_by_tid($tid, 1, $n);
if ($arlist) {
$ids = array();
foreach ($arlist as $val) $ids[] = $val['id'];
well_tag_thread_delete($ids);
}
!empty($tagids) and well_tag_update($tagids, array('count-' => 1));
}
// 标签数据处理 $arr=新提交的数组 $tagarr=保留的旧标签
function well_tag_process($tid, $fid, $new_tags = array(), $tagarr = array())
{
if (empty($tid)) return '';
// 新标签处理入库
if ($new_tags) {
$threadarr = array();
$tagids = array();
$i = 0;
$size = 5;
$n = count($tagarr);
$n = $n > $size ? $size : $size - $n;
foreach ($new_tags as $name) {
++$i;
$name = trim($name);
$name = stripslashes($name);
$name = strip_tags($name);
$name = str_replace(array(' ', '#', "@", "$", "%", "^", '&', '·', '<', '>', ';', '`', '~', '!', '¥', '……', ';', '?', '?', '-', '—', '_', '=', '+', '.', '{', '}', '|', ':', ':', '、', '/', '。', '[', ']', '【', '】', '‘', ' ', ' ', ' ', ' ', ' '), '', $name);
$name = htmlspecialchars($name, ENT_QUOTES);
if ($name && $i <= $n) {
// 查询标签
$read = well_tag_read_name($name);
if ($read) {
// 存在 count+1
$tagids[] = $read['tagid'];
} else {
// 入库
$arr = array('name' => $name, 'count' => 1);
$tagid = well_tag_create($arr);
FALSE === $tagid and message(-1, lang('create_failed'));
$read = array('tagid' => $tagid, 'name' => $name);
}
$tag_thread = array('tagid' => $read['tagid'], 'tid' => $tid);
$threadarr[] = $tag_thread;
$tagarr[$read['tagid']] = $read['name'];
}
}
!empty($threadarr) and tag_thread_big_insert($threadarr);
!empty($tagids) and well_tag_update($tagids, array('count+' => 1));
}
$json = empty($tagarr) ? '' : xn_json_encode($tagarr);
return $json;
}
?> $v = implode(",", $v);
$temp[] = $v;
}
// 去掉重复的字符串,也就是重复的一维数组
$temp = array_unique($temp);
// 再将拆开的数组重新组装
$output = array();
foreach ($temp as $k => $v) {
if ($stkeep) $k = $starr[$k];
if ($ndformat) {
$temparr = explode(",", $v);
foreach ($temparr as $ndkey => $ndval) $output[$k][$ndarr[$ndkey]] = $ndval;
} else $output[$k] = explode(",", $v);
}
return $output;
}
// 合并二维数组 如重复 值以第一个数组值为准
function array2_merge($array1, $array2, $key = '')
{
if (empty($array1) || empty($array2)) return NULL;
$arr = array();
foreach ($array1 as $k => $v) {
isset($v[$key]) ? $arr[$v[$key]] = array_merge($v, $array2[$k]) : $arr[] = array_merge($v, $array2[$k]);
}
return $arr;
}
/*
* 对二维数组排序 两个数组必须有一个相同的键值
* $array1 需要排序数组
* $array2 按照该数组key排序
* */
function array2_sort_key($array1, $array2, $key = '')
{
if (empty($array1) || empty($array2)) return NULL;
$arr = array();
foreach ($array2 as $k => $v) {
if (isset($v[$key]) && $v[$key] == $array1[$v[$key]][$key]) {
$arr[$v[$key]] = $array1[$v[$key]];
} else {
$arr[] = $v;
}
}
return $arr;
}
?>
全国数字地质图1:20W shp格式【保姆级别下载教程】编程频道|福州电脑网
全国数字地质图1:20W shp格式【保姆级别下载教程】
编程之家1220
更新时间:2026-04-03 20:25:37
若转发本文章,请添加引用
文章目录
前言
最近正好下载地层数据,刚开始以为需要手动矢量化,即下载tif后地理配准进行矢量化,研究区较小还好,若区域较大那将会是一个较大的工程,偶然有现成数据可以利用,因此整合网络资源教程给同学们节省时间。
一、地层、岩性
谈论到地层与岩性,需要知道这两个是不太一样的概念,虽然在一定程度上有关联。岩性一般诸如火成岩、沉积岩、变质岩、岩浆岩,地层的话一般按照时代划分,能够分为宇、界、系、统、阶,这其中就能包括岩石特性。
这里我们使用两个网站:
全国地质资料馆https://www.ngac/125cms/c/qggnew/index.htm
地质科学数据出版系统:http://dcc.ngac/cn/page/index
在全国地质资料馆可以查到不同比例尺的地质图,比如1:5W,1:20W,1:25W等,然后再地质科学数据出版系统中找到对应区幅进行下载(由于地质科学数据出版系统中公开了1:20W的资料,因此以1:20W地质图为例)
二、数据下载
1.确定数据地质图幅号
首先在全国地质资料馆—数据与服务—数字资料馆(地质图数据库)找到对应分比例尺地质图
然后打开网站:http://ngac/Distribute/20wanQD.htm,根据研究区经纬度确定图幅及图幅号,比如汶川地区图幅号为H-48-02。然后对其进行检索。检索方式记得点击地址图库。
然后点击进入即可在线游览,利用WTMS服务也可以加载到arcgis中游览,不过速度会比较慢。
2.下载数据
在地质科学数据出版系统中找到全国1∶20W数字地质图(公开版)空间数据库http://dcc.ngac/cn//geologicalData/details/doi/10.23650/data.A.2019.NGA120157.K1.1.1.V1,然后下滑找到对应地质图幅号,需要注意这里需要注册账号,且一个账号每天只能从该网站下载五个数据,若论文使用需要仔细阅读添加引用。下载解压后有三个文件夹。
3.数据转换
该数据格式为mapgis格式,下一步需要将其转化为我们常用的shp格式。
Mapgis转换方法:
(1)Mapgis软件转换(但是目前破解版Mapgis67较多,自己也尝试后发现转换后属性数据部分会有偏差,有些显示错误,改善方法是使用高版本Mapgis)
(2)使用FME软件转化。
这里用第二种方法进行转换,参考https://zhuanlan.zhihu/p/371786900,使用FME2020+MyFME插件进行转换(MyFME为必须,否则无法识别Mapgis格式)
按照文章及压缩包内教程后,打开FME Workbench:(记得修改时间!)
网盘链接
1.右上角文件—新建—空工作空间
2.工具栏选择读模块,数据集选择解压数据后MAPGIS文件夹下的所有数据。
3.工具栏选择写模块,格式选择esrishapeflie ,然后确定后勾选全选、排序。
4.连线,将对应的左右灰色箭头拉向红色箭头,最后点击运行即可,正常运行不会报错。
总结
最后打开arcgis查看效果,发现结果还不错,数据包括断层数据、地层数据、侵入岩数据等,较大程度上可以避免矢量化,但是转换时没有选择参考下,使用时需要定义投影CGS——WGS84。
本文发布于:2024-10-17,感谢您对本站的认可!
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
本文标签:地质图保姆级别格式数字
发布评论