CN104581000A - 一种视频关注目标的运动轨迹的快速检索方法 - Google Patents

一种视频关注目标的运动轨迹的快速检索方法 Download PDF

Info

Publication number
CN104581000A
CN104581000A CN201310475493.XA CN201310475493A CN104581000A CN 104581000 A CN104581000 A CN 104581000A CN 201310475493 A CN201310475493 A CN 201310475493A CN 104581000 A CN104581000 A CN 104581000A
Authority
CN
China
Prior art keywords
video
camera
target
database
information
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201310475493.XA
Other languages
English (en)
Inventor
郝文娟
王安
杜海涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Changfeng Science Technology Industry Group Corp
Original Assignee
China Changfeng Science Technology Industry Group Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Changfeng Science Technology Industry Group Corp filed Critical China Changfeng Science Technology Industry Group Corp
Priority to CN201310475493.XA priority Critical patent/CN104581000A/zh
Publication of CN104581000A publication Critical patent/CN104581000A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Closed-Circuit Television Systems (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Television Signal Processing For Recording (AREA)

Abstract

本发明涉及一种视频关注目标的运动轨迹的快速检索方法,包括以下步骤:(1)整合海量视频,关联地理空间信息;(2)异构数据转换,生成统一格式的视频流进行视频查看;(3)数据库预先设定的摄像头拓扑关系表,遍历所有出现目标的摄像头;(4)对所有出现目标的摄像头的视频图片进行要素编辑记录,保存到数据库;(5)根据要素记录,按时间排序拼接关注目标的整个运动轨迹视频;同时在GIS地图上按时间顺序形成一条目标运动轨迹路线;(6)导出统一格式的目标运动轨迹视频路线图,便于存档及以后查看。本发明可快速发现目标,并直观地获取其空间运动轨迹,提高检索效率,节省查看时间。

Description

一种视频关注目标的运动轨迹的快速检索方法
技术领域
本发明涉及一种视频关注目标的运动轨迹的快速检索方法,当目标在一定监控区域内连续运动时,可通过同时审看若干个摄像头视频来直接掌握关注目标的运动轨迹,属于视频监控领域。
背景技术
随着“平安城市”、“智慧城市”等系统建设陆续开展,视频监控得到了大规模应用,积累了海量的视频数据资源,为相关部门查找关注目标提供了有效资源,在维护公共安全方面起到了一定作用。但是当需要掌握某个关注目标的运动轨迹时,如何在海量的视频数据资源中,快速地发现关注目标的实际空间运动轨迹,成为了视频监控行业当前面临的一个巨大挑战。由于目前视频监控设备生产厂家众多,以及人工查看海量视频资源的局限性,从海量视频中有效地提取出感兴趣目标的运动轨迹还是存在一定的困难,如:
(1)监控视频分散且格式不一,视频的统一导入导出困难:目前社会上所使用的监控视频来自不同厂家,它们所采用的数据编码、封装格式并不一样,有时需要使用不同的播放器来查看视频,而且监控网络分布比较分散,需要到现场逐一导出视频。
(2)人工查看视频工作繁重,效率低下:大量的视频资源需要采用“人海战术”,人工长时间专注于同一视频资源,容易疲劳,注意力无法集中,有可能忽略视频中的某些重要目标。
(3)视频资源数据归档困难,线索与视频对应混乱:目前案件归档主要采用的文件目录管理方式,这需要通过其他方法将视频的重要内容作记录,造成了线索与视频无法对应。
(4)视频线索时空信息不足,无法直观获取目标运动轨迹:通常关注目标的行为过程会涉及到多个摄像头,在分析目标运动轨迹时都需要人工参与,缺乏有效地理分析工具及监控分布点示意图。
发明内容
本发明的目的是结合海量视频数据与地理信息技术,提供一种视频关注目标的运动轨迹的快速检索方法,帮助相关人员快速发现目标,并直观地获取其空间运动轨迹,提高检索效率,节省查看时间。
本发明的技术方案如下:
一种视频关注目标的运动轨迹的快速检索方法,其特征在于,包括以下步骤:
(1)整合海量视频,关联地理空间信息:建立数据库,将不同监控网络的摄像头位置信息,包括摄像头ID、视频储存器IP地址及端口号信息,及其关联的视频数据信息储存在数据库中,并建立每一个摄像头的东、南、西、北、东南、东北、西北、西南八个方向的相邻摄像头的拓扑关系表;在GIS地图上显示相应位置的摄像头;
(2)异构数据转换:从数据库获取相关视频资源的摄像头ID、视频储存器IP地址及端口号信息,与相应的视频存储装置建立连接,读取指定视频,并对此视频加以解码生成统一格式的视频流进行视频查看;
(3)获取关注目标的摄像头视频:在视频查看时,将首次发现关注目标的摄像头作为当前关注摄像头,通过读取数据库预先设定的摄像头拓扑关系表,获取当前关注摄像头的视频及其八个方向的相邻摄像头的视频资源同时进行查看,九路视频同步播放,当发现某个方向的摄像头的视频中出现目标时,则将该方向的摄像头作为新的当前关注摄像头,再获取该新的当前关注摄像头及其八个方向的相邻摄像头的视频资源同时进行查看,如此循环,直至遍历所有出现目标的摄像头;
(4)对所有出现目标的摄像头的视频图片进行要素编辑记录,保存到数据库;
(5)根据要素记录,按时间排序拼接关注目标的整个运动轨迹视频;同时在GIS地图上按时间顺序形成一条目标运动轨迹路线,实现目标运动轨迹在实际空间中的映射。
所述方法进一步还可包括以下步骤:
(6)导出统一格式的目标运动轨迹视频路线图,便于存档及以后查看。
上述步骤(4)中所述的要素,包含时间、地点、人物、车辆、事件,对所述要素进行编辑记录,将视频资源与其中的要素信息直接关联,形成一对一的关系,保存到数据库;当相关人员查看某一事件时,会方便地获取到所有相关联的视频资源及目标信息记录内容。
本发明针对社会中的突发公共安全事件在处理时所遇到的海量视频数据提取关注目标困难、人力消耗较大、缺乏视频的时空信息等问题,结合监控视频资源与地理信息技术,采用视频拓扑关联回查方法对目标进行行进跟踪,形成一条具有时间性的关注目标的空间运动轨迹,导出统一格式的视频轨迹录像及路线。此方法很大程度上减轻了线索目标的查询检索工作,缩短了目标运动轨迹检索的时间,提高了处理效率,更好地服务于平安城市、智慧城市建设等安防行业。
附图说明
图1是本发明视频关注目标的运动轨迹的快速检索方法流程图。
具体实施方式
本发明的视频关注目标的运动轨迹的快速检索方法主要包括四部分:海量视频数据的异构接入、视频拓扑九宫格回查、要素编辑记录、目标运动轨迹的时空分析。另外,为保证算法有效实现,需要将不同监控网络的摄像头位置信息及其关联的视频数据信息储存在数据库中,而且数据库中要维护着每一个摄像头与其他摄像头的拓扑关系,并结合地理信息在地图上显示相应位置的摄像头,即实现了视频资源与地理信息的有效结合。
海量视频数据的异构接入提供了算法实现的基础资源,它是根据指令获取所关注的视频资源地址等信息,直接连接相应储存器,读取相关视频,并对其进行解码获取算法所需要的统一视频流,排除了视频格式不一的干扰。
在得到统一视频流基础上,根据数据库中维护的摄像头拓扑关系,可以同时浏览所关注摄像头视频资源与其八个方向(东、南、西、北、东南、东北、西北、西南)相邻摄像头视频资源,即所谓的九宫格关联回查。当在某方向摄像头视频中发现了可疑目标时,也可以将其选为关注摄像头资源,并同时浏览其拓扑摄像头资源,如此往复。这代替了逐一审看视频操作,提高了视频审看效率。
对发现可疑目标的视频图片进行要素编辑记录并将其保存到数据库,实现了将监控视频与视频相关重要记录的统一管理保存,便于以后查看及归档。
目标运动轨迹的时空分析是根据相应的要素记录,对其关联视频按时间排序,并在地理信息相应位置上做标记,最终形成一条具有时间性的目标运动轨迹,使人不需要去现场就能够直观地掌握目标的行动轨迹。
图1是本发明视频关注目标的运动轨迹的快速检索方法流程图。具体实现方式如下:
100整合海量视频,关联地理空间信息:建立数据库,将不同监控网络的摄像头位置信息,包括摄像头ID、视频储存器IP地址及端口号信息,及其关联的视频数据信息储存在数据库中,并建立每一个摄像头的东、南、西、北、东南、东北、西北、西南八个方向的相邻摄像头的拓扑关系表;在GIS地图上显示相应位置的摄像头;
异构数据转换:从数据库获取相关视频资源的摄像头ID、视频储存器IP地址及端口号信息,与相应的视频存储装置建立连接,读取指定视频,并对此视频加以解码生成统一格式的视频流进行视频查看;
110获取关注目标的摄像头视频:在视频查看时,将首次发现关注目标的摄像头作为当前关注摄像头;
120通过读取数据库预先设定的摄像头拓扑关系表,获取当前关注摄像头的视频及其八个方向的相邻摄像头的视频资源同时进行查看;
130九路视频同步播放;
140及150当发现某个方向的摄像头的视频中出现目标时,则将该方向的摄像头作为新的当前关注摄像头,再获取该新的当前关注摄像头及其八个方向的相邻摄像头的视频资源同时进行查看,如此循环,直至遍历所有出现目标的摄像头;
160对所有出现目标的摄像头的视频图片进行要素编辑记录,保存到数据库;所述的要素,包含时间、地点、人物、车辆、事件;对所述要素进行编辑记录,将视频资源与其中的要素信息直接关联,形成一对一的关系,保存到数据库,当相关人员查看某一事件时,会方便地获取到所有相关联的视频资源及目标信息记录内容;
170根据要素记录,按时间排序拼接关注目标的整个运动轨迹视频;
180同时在GIS地图上按时间顺序形成一条目标运动轨迹路线,实现目标运动轨迹在实际空间中的映射。
190导出统一格式的目标运动轨迹视频路线图,便于存档及以后查看。

Claims (3)

1.一种视频关注目标的运动轨迹的快速检索方法,其特征在于,包括以下步骤:
(1)整合海量视频,关联地理空间信息:建立数据库,将不同监控网络的摄像头位置信息,包括摄像头ID、视频储存器IP地址及端口号信息,及其关联的视频数据信息储存在数据库中,并建立每一个摄像头的东、南、西、北、东南、东北、西北、西南八个方向的相邻摄像头的拓扑关系表;在GIS地图上显示相应位置的摄像头;
(2)异构数据转换:从数据库获取相关视频资源的摄像头ID、视频储存器IP地址及端口号信息,与相应的视频存储装置建立连接,读取指定视频,并对此视频加以解码生成统一格式的视频流进行视频查看;
(3)获取关注目标的摄像头视频:在视频查看时,将首次发现关注目标的摄像头作为当前关注摄像头,通过读取数据库预先设定的摄像头拓扑关系表,获取当前关注摄像头的视频及其八个方向的相邻摄像头的视频资源同时进行查看,九路视频同步播放,当发现某个方向的摄像头的视频中出现目标时,则将该方向的摄像头作为新的当前关注摄像头,再获取该新的当前关注摄像头及其八个方向的相邻摄像头的视频资源同时进行查看,如此循环,直至遍历所有出现目标的摄像头;
(4)对所有出现目标的摄像头的视频图片进行要素编辑记录,保存到数据库;
(5)根据要素记录,按时间排序拼接关注目标的整个运动轨迹视频;同时在GIS地图上按时间顺序形成一条目标运动轨迹路线,实现目标运动轨迹在实际空间中的映射。
2.根据权利要求1所述的方法,其特征在于,所述方法进一步包括:
(6)导出统一格式的目标运动轨迹视频路线图,便于存档及以后查看。
3.根据权利要求1所述的方法,其特征在于:步骤(4)中所述的要素,包含时间、地点、人物、车辆、事件;对所述要素进行编辑记录,将视频资源与其中的要素信息直接关联,形成一对一的关系,保存到数据库,当相关人员查看某一事件时,会方便地获取到所有相关联的视频资源及目标信息记录内容。
CN201310475493.XA 2013-10-12 2013-10-12 一种视频关注目标的运动轨迹的快速检索方法 Pending CN104581000A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310475493.XA CN104581000A (zh) 2013-10-12 2013-10-12 一种视频关注目标的运动轨迹的快速检索方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310475493.XA CN104581000A (zh) 2013-10-12 2013-10-12 一种视频关注目标的运动轨迹的快速检索方法

Publications (1)

Publication Number Publication Date
CN104581000A true CN104581000A (zh) 2015-04-29

Family

ID=53096081

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310475493.XA Pending CN104581000A (zh) 2013-10-12 2013-10-12 一种视频关注目标的运动轨迹的快速检索方法

Country Status (1)

Country Link
CN (1) CN104581000A (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105336171A (zh) * 2015-10-16 2016-02-17 浙江宇视科技有限公司 一种摄像机位置标定方法及装置
CN105681761A (zh) * 2016-03-09 2016-06-15 国网浙江省电力公司湖州供电公司 一种基于菊花链拓扑的输变电设备状态监测系统
CN106096577A (zh) * 2016-06-24 2016-11-09 安徽工业大学 一种摄像头分布地图中的目标追踪系统及追踪方法
CN106777078A (zh) * 2016-12-13 2017-05-31 广东中星电子有限公司 一种基于信息数据库的视频检索方法和系统
CN107004271A (zh) * 2016-08-22 2017-08-01 深圳前海达闼云端智能科技有限公司 显示方法、装置、电子设备、计算机程序产品和非暂态计算机可读存储介质
CN107948586A (zh) * 2017-11-14 2018-04-20 山东大学 基于视频拼接的跨区域运动目标检测方法和装置
CN108234961A (zh) * 2018-02-13 2018-06-29 欧阳昌君 一种多路摄像机编码与视频流引流方法及系统
CN109871815A (zh) * 2019-02-26 2019-06-11 新华三技术有限公司 一种查询监控信息的方法及装置
CN110188691A (zh) * 2019-05-30 2019-08-30 银河水滴科技(北京)有限公司 一种移动轨迹确定方法及装置
CN110245268A (zh) * 2019-06-26 2019-09-17 银河水滴科技(北京)有限公司 一种路线确定、展示的方法及装置
CN110399525A (zh) * 2019-07-27 2019-11-01 广州伽利略网络科技有限公司 一种基于图像识别技术的人、物运动轨迹快速检索方法
CN110532432A (zh) * 2019-08-21 2019-12-03 深圳供电局有限公司 一种人物轨迹检索方法及其系统、计算机可读存储介质
CN110532923A (zh) * 2019-08-21 2019-12-03 深圳供电局有限公司 一种人物轨迹检索方法及其系统
CN110543457A (zh) * 2019-09-11 2019-12-06 北京明略软件系统有限公司 轨迹类文档处理方法和装置、存储介质及电子装置
WO2020135523A1 (zh) * 2018-12-29 2020-07-02 华为技术有限公司 目标对象的检索定位方法和装置
CN111372036A (zh) * 2018-12-25 2020-07-03 中国科学院深圳先进技术研究院 一种视频监测方法、系统
CN112256819A (zh) * 2020-11-13 2021-01-22 珠海大横琴科技发展有限公司 一种基于监控地图的数据处理方法、装置和可读存储介质
CN112449160A (zh) * 2020-11-13 2021-03-05 珠海大横琴科技发展有限公司 一种视频监控方法、装置和可读存储介质
CN112464757A (zh) * 2020-11-16 2021-03-09 复旦大学 一种基于高清视频的目标实时定位和轨迹重构方法
CN112511864A (zh) * 2020-11-23 2021-03-16 北京爱笔科技有限公司 轨迹展示方法、装置、计算机设备和存储介质
CN113497915A (zh) * 2020-04-01 2021-10-12 苏州迪沃蒙特智能科技有限公司 一种视频监控系统
CN113837023A (zh) * 2021-09-02 2021-12-24 北京新橙智慧科技发展有限公司 一种跨摄像头行人自动追踪方法
CN114422828A (zh) * 2022-01-26 2022-04-29 复旦大学 一种基于高清摄像头的实时定位和轨迹重构方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110302130A1 (en) * 2010-06-03 2011-12-08 Electronics And Telecommunications Research Institute Intelligent video information retrieval apparatus and method capable of multi-dimensional video indexing and retrieval
CN102436662A (zh) * 2011-11-29 2012-05-02 南京信息工程大学 一种非重叠视域多摄像机网络中的人体目标跟踪方法
CN102930061A (zh) * 2012-11-28 2013-02-13 安徽水天信息科技有限公司 一种基于运动目标检测的视频摘要方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110302130A1 (en) * 2010-06-03 2011-12-08 Electronics And Telecommunications Research Institute Intelligent video information retrieval apparatus and method capable of multi-dimensional video indexing and retrieval
CN102436662A (zh) * 2011-11-29 2012-05-02 南京信息工程大学 一种非重叠视域多摄像机网络中的人体目标跟踪方法
CN102930061A (zh) * 2012-11-28 2013-02-13 安徽水天信息科技有限公司 一种基于运动目标检测的视频摘要方法及系统

Cited By (29)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105336171B (zh) * 2015-10-16 2017-12-29 浙江宇视科技有限公司 一种摄像机位置标定方法及装置
CN105336171A (zh) * 2015-10-16 2016-02-17 浙江宇视科技有限公司 一种摄像机位置标定方法及装置
CN105681761A (zh) * 2016-03-09 2016-06-15 国网浙江省电力公司湖州供电公司 一种基于菊花链拓扑的输变电设备状态监测系统
CN106096577A (zh) * 2016-06-24 2016-11-09 安徽工业大学 一种摄像头分布地图中的目标追踪系统及追踪方法
CN106096577B (zh) * 2016-06-24 2019-05-31 安徽工业大学 一种摄像头分布地图中的目标追踪方法
CN107004271A (zh) * 2016-08-22 2017-08-01 深圳前海达闼云端智能科技有限公司 显示方法、装置、电子设备、计算机程序产品和非暂态计算机可读存储介质
CN106777078A (zh) * 2016-12-13 2017-05-31 广东中星电子有限公司 一种基于信息数据库的视频检索方法和系统
CN106777078B (zh) * 2016-12-13 2020-12-22 广东中星电子有限公司 一种基于信息数据库的视频检索方法和系统
CN107948586A (zh) * 2017-11-14 2018-04-20 山东大学 基于视频拼接的跨区域运动目标检测方法和装置
CN108234961A (zh) * 2018-02-13 2018-06-29 欧阳昌君 一种多路摄像机编码与视频流引流方法及系统
CN108234961B (zh) * 2018-02-13 2020-10-02 欧阳昌君 一种多路摄像机编码与视频流引流方法及系统
CN111372036A (zh) * 2018-12-25 2020-07-03 中国科学院深圳先进技术研究院 一种视频监测方法、系统
WO2020135523A1 (zh) * 2018-12-29 2020-07-02 华为技术有限公司 目标对象的检索定位方法和装置
CN109871815A (zh) * 2019-02-26 2019-06-11 新华三技术有限公司 一种查询监控信息的方法及装置
CN109871815B (zh) * 2019-02-26 2021-07-20 新华三技术有限公司 一种查询监控信息的方法及装置
CN110188691A (zh) * 2019-05-30 2019-08-30 银河水滴科技(北京)有限公司 一种移动轨迹确定方法及装置
CN110245268A (zh) * 2019-06-26 2019-09-17 银河水滴科技(北京)有限公司 一种路线确定、展示的方法及装置
CN110399525A (zh) * 2019-07-27 2019-11-01 广州伽利略网络科技有限公司 一种基于图像识别技术的人、物运动轨迹快速检索方法
CN110532923A (zh) * 2019-08-21 2019-12-03 深圳供电局有限公司 一种人物轨迹检索方法及其系统
CN110532432A (zh) * 2019-08-21 2019-12-03 深圳供电局有限公司 一种人物轨迹检索方法及其系统、计算机可读存储介质
CN110543457A (zh) * 2019-09-11 2019-12-06 北京明略软件系统有限公司 轨迹类文档处理方法和装置、存储介质及电子装置
CN113497915A (zh) * 2020-04-01 2021-10-12 苏州迪沃蒙特智能科技有限公司 一种视频监控系统
CN112256819A (zh) * 2020-11-13 2021-01-22 珠海大横琴科技发展有限公司 一种基于监控地图的数据处理方法、装置和可读存储介质
CN112449160A (zh) * 2020-11-13 2021-03-05 珠海大横琴科技发展有限公司 一种视频监控方法、装置和可读存储介质
CN112464757A (zh) * 2020-11-16 2021-03-09 复旦大学 一种基于高清视频的目标实时定位和轨迹重构方法
CN112511864A (zh) * 2020-11-23 2021-03-16 北京爱笔科技有限公司 轨迹展示方法、装置、计算机设备和存储介质
CN112511864B (zh) * 2020-11-23 2023-02-17 北京爱笔科技有限公司 轨迹展示方法、装置、计算机设备和存储介质
CN113837023A (zh) * 2021-09-02 2021-12-24 北京新橙智慧科技发展有限公司 一种跨摄像头行人自动追踪方法
CN114422828A (zh) * 2022-01-26 2022-04-29 复旦大学 一种基于高清摄像头的实时定位和轨迹重构方法

Similar Documents

Publication Publication Date Title
CN104581000A (zh) 一种视频关注目标的运动轨迹的快速检索方法
CN107370983B (zh) 用于视频监控系统的行踪轨迹的获取方法和装置
CN103839194B (zh) 基于输电线路和gis的无人机巡检图像检索系统及方法
CN102929903B (zh) 一种基于视频信息分层结构化描述的快速视频检索方法
CN106156199B (zh) 一种视频监控图像存储检索方法
CN103679730A (zh) 基于gis的视频摘要生成方法
CN105335387A (zh) 一种视频云存储系统的检索方法
WO2016119368A1 (zh) 目标跟踪方法及装置
TWI450207B (zh) 物件追蹤的方法、系統、電腦程式產品與記錄媒體
CN103794234A (zh) 基于海量视频的事件轨迹快速检索平台
CN103530326A (zh) 在视频文件中添加地理位置信息并建立索引的方法
CN102411765A (zh) 三维电力网络的建设方法及装置
CN104361586A (zh) 基于全景图的电力杆塔维护与预警系统
CN109698943B (zh) 一种容量可伸缩的网络摄像机接入平台
CN102595103A (zh) 一种基于gis地图推演智能视频的方法
US20140085485A1 (en) Machine-to-machine enabled image capture and processing
CN103260010A (zh) 智慧天网视频快速检索系统
CN106534784A (zh) 一种用于视频分析数据结果集的采集分析存储统计系统
CN104486585A (zh) 一种基于gis的城市海量监控视频管理方法及系统
CN113255477A (zh) 一种行人视频图像的综合管理系统及方法
CN105472333A (zh) 视频监控设备拓扑系统的建立方法及关联监控方法
Alfarrarjeh et al. Spatial coverage measurement of geo-tagged visual data: A database approach
CN104426937A (zh) 通过移动电话和云计算实现对其他设备所记录内容的定位
Lin et al. Moving camera analytics: Emerging scenarios, challenges, and applications
Wu et al. An adaptive organization method of geovideo data for spatio-temporal association analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20150429

WD01 Invention patent application deemed withdrawn after publication