Method Article
우리는 역동적 인 활동을하는 동안 밀리 초 이하의 시간 분해능을 가진 근육의 수축 속도, 변형 및 변형 속도를 측정하는 새로운 초음파 기반의 벡터 조직 도플러 영상 기술에 대해 설명합니다. 이 방법은 동적 인 근육 기능의 보완적인 측정을 제공하고 근골격계 질환을 기본 메커니즘에 대한 이해로 이어질 수 있습니다.
초음파는 동적 인 작업을하는 동안 영상 근육과 힘줄의 움직임에 대한 매력적인 양상과 임상 또는 실험실 환경에서 생체 역학 연구를위한 상호 보완적인 방법 론적 접근 방식을 제공 할 수 있습니다. 이 목표를 향해, 초음파 영상에서 근육 운동 역학의 정량 방법이 개발 한 이미지 프로세싱을 기반으로하고 있습니다. 이러한 방법의 시간 해상도는 드롭 방문으로 매우 역동적 인 작업을 위해 일반적으로 충분하지 않습니다. 우리는 근육의 운동 역학을 정량화하기위한 도플러 방법을 활용하는 새로운 접근 방식을 제안한다. 우리는 새로운 벡터 조직 도플러 영상 초음파를 사용하여 동적 활동을하는 동안 밀리 초 이하의 시간 해상도와 근골격계 수축 속도, 변형 및 변형 속도를 측정 할 수있다 (vTDI) 기술을 개발했다. 이 예비 연구의 목적은 근골격계를 내보 측정의 반복성과 vTDI 기술의 적용 가능성을 알아보고자 하였다건강한 주제에있는 드롭 랜딩 작업 동안 ities. vTDI 측정은 공동 기구학 및 동역학에 대한 3 차원 모션 캡처,지면 반력에 대한 근육의 활성화와 힘 플레이트의 타이밍에 대한 근전도와 같은 다른 생체 역학 기술과 동시에 수행 할 수 있습니다. 이러한 보완적인 기술의 통합은 근골격계 질환의 병인 및 병태 생리를 기본 동적 근육의 기능과 장애의 더 나은 이해로 이어질 수 있습니다.
근골격계 질환은 성인 1에서 널리 유행이다. 그들은 미국 2에서 선도적 인 만성 질환이며, 전 세계적으로 3 사람의 25 %에 영향을 미칠 것으로보고있다. 근골격계 질환은 일상 생활 (ADL), 기능 제한과 수명 4의 낮은 품질의 활동을 감소 기능과 연관된다. 경제적 부담 때문에 생산성 손실과 높은 의료 비용 4 중요하다. 이러한 장애의 여러 가지의 병태 생리는 부적절하게 이해 상태를 유지합니다. 예를 들어, 전방 십자 인대 (ACL) 부상의 재건 다음과 골관절염 (OA) 4의 병인은 대퇴사 두근의 근력과 기능 5 변경에 연결되어 있습니다,하지만 기본 메커니즘은 불분명하다. 기본 메커니즘을 설명하기 위하여, 더 나은 동적 근육 기능을 이해하기 위해 필요하다.
기능개별 근육의 평가는, 부분 또는 ADL과 활동적인 라이프 스타일 (예 : 스포츠)과 관련된 전체 작업의 성능을하는 동안 근육의 기능과 이러한 질환의 병인 및 병태 생리에서의 잠재적 인 역할에 대한 자세한 통찰력을 제공 할 수 있습니다. 또한 재활 동안 근육 기능 개선 정량화 결과 측정으로 사용될 수있다. 병원에서 근육과 관절 기능을 측정하는 종래의 기술은 운동 범위, 근력 및 / 또는 근육 그룹 지구력 등의 체력 검사를 포함한다. 현재 병원에서, 근전도 (EMG)를 근육 활성화 / 공동 활성화, 주파수 및 근육 활동의 진폭을 평가하는 데 사용된다. 그러나, EMG는 근육의 전기 활성의 측정은 반드시 근력, 수축 능력 및 근육의 다른 기능 요소에 대한 정보를 제공하지 않는다. 같은 3D 모션 캡쳐 시스템 F와 같은 다른 복잡한 생체 역학적 평가,또는 접지 반력 조인트 동역학 및 운동학 및 포스 판 보용 실험실 6-9에서 수행 될 수있다. 이들 기술에 의해 측정 조인트 레벨에 반드시 동적 또는 기능적 활동을하는 동안 개별 근육 기능의 이해를 직접 제공하지 않는다. 동적 활동을 수행하면서 동시에 근육의 이미징을 수행 할 수있는 능력은 잠재적으로 근육의 수준에서보다 효율적이고 현실적인 기능 평가로 이어질 수 있습니다.
연구의 대부분은 정적 인 경향이 위치에있는 근육의 기능에 초점을 맞추고있다,이 방법은 실시간 상황에서 근육의 행동에 대한 우리의 이해를 향상시키기 위해 새로운 길을 열 수 있습니다.
진단 초음파는 실시간으로 근육과 힘줄의 직접 이미지를 사용할 수 있습니다, 따라서 ADL 동안 근골격계 역학과 기능을 측정하기위한 매력적인 대안이다. 의 초음파 기반의 정량적 측정근육의 두께, 길이, 폭, 단면적 (CSA), 섬유 pennation 각도와 다발의 길이와 근육의 형태와 구조는, 널리 10-12을 사용하고 있습니다. 최근 몇 년 동안, 이미지 처리 방법은 동적 인 작업 13-14 동안 이러한 양적 조치를 평가하고 정량화하기 위해 사용되어왔다. 이러한 발전은 생체 근육의 기능에 대한 이해에 대한 새로운 방법 론적 접근 방식을 사용할 수있다. 그러나, 이러한 방법은 주로 기존의 그레이 스케일 (또는 B-모드) 초음파 영상을 사용에 의존, 따라서 충분히 가치있는 것으로 표시되었습니다 조직 속도, 변형 및 도플러 원리를 사용하여 변형 속도를 측정하는 초음파의 가능성을 악용하지 않은 심장 근육 기능 15-16을 평가.
우리는 높은 시간 해상도 (하위 millisecon와 수축 속도, 변형 및 변형 속도를 측정 할 수있는 벡터 조직 도플러 영상 (vTDI) 기술을 개발 한D) 동적 활동 17 ~ 18시. 특히, vTDI 기술은 높은 프레임 속도에서 매우 역동적 인 작업 (예 : 드롭 방문, 걸음 걸이 등)을하는 동안 근육과 힘의 측정을 수행 할 수 있습니다. vTDI 기술은 초음파 빔을 따른 속도 성분만을 추정 종래 도플러 초음파 비해 개선하고 insonation 각도에 따라서 좌우된다. vTDI는 서로 다른 각도로 스티어링 된 두 개의 서로 다른 초음파 빔을 이용하여 근육 및 힘줄의 속도를 추정하고, 결상면에서 insonation 각도 따라서 독립적이다. 근육 수축은 3D로 발생하기 때문에, 결상면의 만곡은 여전히 중요하다. 우리는 임상 환경에서 만들어지는 이러한 측정을 가능하게 연구 인터페이스를 시판 초음파 시스템에서이 방법을 구현했습니다.
반복성과 vTDI의 SYST의 적용 가능성을 조사하려면EM은 동적 작업하는 동안 대퇴 직근 근육 속도를 측정, 우리는 건강한 성인 지원자에 예비 연구를 시행 하였다. 이 논문은 방법론과 수축 속도를 추정하기위한 실험 장치를 보여줍니다, 변형 및 직근의 변형 속도는 드롭 랜딩 작업 중에 밀리 초 이하의 시간 분해능을 가진 근육을 대퇴.
1. 수단
벡터 TDI는 두 개 이상의 독립적 인 경로에서 찍은 도플러 속도 측정에서 얻어진 속도 벡터를 추정에 기초한다. 연구 인터페이스 초음파 시스템은 vTDI 개발에 사용 하였다. 연구 인터페이스는 소프트웨어 개발 키트 (SDK)를 사용하여 저역 빔 포밍 및 펄스 시퀀스 제어를 허용했다. 128 트랜스 듀서 요소와보기의 38mm 필드로 구성된 5-14 MHz의 선형 배열 센서가 사용되었다. 연구 인터페이스는 두 개의 전송로 배열 변환기를 분리하고 구멍을 받고 일반에 대하여 15 °로 수신 빔을 조종하기 위해 사용되었다. 송신 빔은 관심 영역 (예를 들면 근육 배)에 집중되었다. 전송 및 수신 구멍 32 요소로 설정 하였다.
여덟 과목, 4 남자와 여자 4 명 (29.7 ± 6.5 세)을 대상으로 모집했다. 의 과목에서 운동 학적 측정오른쪽 하반신는 고속 성능 및 200 Hz의 샘플링 속도로 팔 카메라 모션 캡쳐 시스템을 사용하여 포착 하였다. 실험 기간 동안 그라운드 반력 데이터는 2000 Hz에서이 힘 플레이트 샘플링을 통해 얻어졌다.
고속 카메라를 삼각대에 장착하고 피사체로부터 2m에 배치는 500 프레임 / 초에서 드롭 착륙을 캡처하는 데 사용되었다.
2. 제목 준비
3. 실험 프로토콜
4. 초음파 데이터 분석
5. 3D 모션 캡처 데이터 분석
6. 고속 카메라 데이터 분석
방법을 보여주는 우리의 이전 직장에서 대표 결과는 다음과 같다. 우리의 현재의 연구에 사용 된 방법은 영상과 모션 캡처를 통합하는 동안, 아래에 제시 대표적인 결과는 이러한 측정을 별도로 수행 된 연구에서이다.
I. 초음파 (vTDI)
3D 모션 캡처 및 고속 카메라, 피사체의 점프의 패턴의 데이터를 사용하여, 방문 및 안정화 단계는 각 시험에 대한 공부했다. 축 방향 및 외 직근이 vTDI에서 근육 속도를 대퇴은 3D 모션 캡처 및 고속 카메라에서 수집 된 데이터를 비교 하였다. 이 데이터를 사용하여, 드롭 착륙 순서에 걸쳐 축 방향 및 횡 방향 대퇴 직근 근육 속도의 시간적 특성을 연구 하였다. 긍정적 인 측면 속도는, 무릎 굴곡시 대퇴 직근 근육의 편심 수축에 대응부정적인 측면 속도는 무릎 확장하는 동안 근육의 동심 수축에 대응하면서. 이것은도 2에 도시되어있다. 모든 주제에 대한 전체 드롭 방문 순서는 약 1.45 ± 0.27 초 지속되었다.
각 과목의 경우, 축 방향 및 횡 근육 속도는 0.99의 기울기 및 R2가 = 825 (도 3)과 실험 사이의 강한 반복성을 보였다. 두 과목 (두 사람이) 더 높은 속도가 있었을 때 8 번 중 6 번 주제에 대한 속도 값은 48-62cm / sec의 유사 범위에 있었다. 각 과목의 개별 체중과 근육의 두께에 대한 조정 때 남성은 (72.96 cm / 초)이 여성보다 훨씬 높은 근육 속도 (48.71 cm / 초), P = 0.029을 발표했다.
초음파 트랜스 듀서의 위치는 고속 카메라를 사용하여 드롭 방문 순서를 생각 추적 하였다. 대전 자부와 커프스 (녹색 대시 사이의 선분 사이의 각도에드 선분) 중반 허벅지와 팔목 (자주색 점선 구간) 사이의 라인 세그먼트 (segment)를 계산 하였다. 주제 당 2 시험 16 시험의 총, (재판 1 & 2는 1 주제와 관련 등) 그림 4에서 관찰된다. 착륙시 해부학 적 마커 변환기 홀더 상대의 최소 각도 변화 (0.91 ° ± 0.54도)는 16 실험을 통해 관찰되었다. 초음파 변환기 각도 변화가 (ICC 2,1 = 0.90, P <0.05)뿐만 아니라 높은 반복성을 발표했다.이 착륙 시험 중에 변환기 운동은 최소한의 것을 보여주고 속도 측정으로 인해 모든 변환기의 움직임에 영향을받지 않았다.
II. 3D 모션 카메라 · 포스 플레이트
우리는 주로 무릎과 엉덩이 굴곡 각도, 무릎 외반 각도와 무릎의 외반 순간에 집중했다. 엉덩이 F : 우리는지면과의 초기 접촉시, 주제는 다음과 같은 운동 패턴을주었습니다LEXION 41 ° ± 13도, 무릎 굴곡 23 ° ± 9도, 무릎 외반 0.03 ° ± 6도. 그들이 착륙 단계 진행으로, 달성의 최대 각도는 : 엉덩이 굴곡 58 ° ± 19도, 무릎 굴곡 54 ° ± 24도, 무릎 외반 -4 ± 8도 (그림 5) °. 무릎 외반 모멘트가 착륙 단계 (그림 6) 동안 최대로 초기 접지에서 0.1 ± 0.1 ㎚ / km에 0.03 ± 0.03에서 감소를 발표했다.
그림 1. 직근의 vTDI 속도 측정의 표현은 근육 대퇴. 회색 빔이 두 개의 개별 전송을 대표하는 빔과 레드 라인 (무릎의 근위-말단 방향) 횡 방향 속도 성분을 나타내며, 파란색 선은 축 속도를 나타냅니다를받을 (근육의 두께에 따라) 구성 요소입니다.
그림 2. 드롭 착지시 축 방향과 횡 속도는 비디오 프레임들의 시퀀스 (상부 패널)에 비교된다. 하판은 초기 무릎의 굴곡에 대응하는 축 방향 및 측 방향 속도가이고, B는 무릎 확장자에 대응하고, C는에 대응 발가락 바닥에 충격을 가한 경우, D는 땅을 치는 발 뒤꿈치에 해당하고, E는 무릎 굴곡 포스트 착륙 해당하고 F 무릎 확장 및 안정화에 해당합니다.
그림 3. 8 과목 (과목당 2 시험)의 결과 속도 벡터의 크기의 반복은. 남자는 파란색 동그라미에 붉은 다이아몬드와 여자로 표시되어 있습니다.
그림 4. 패널 A. 오류초음파 변환기 홀더 제 선분 및 허벅지 중간 (보라색 점선 부분)에 대한 마커 및 초음파 변환기에 의해 선분과 대전 자부의 마커 사이의 각도 (녹색 선분을 점선). 패널 B. 허벅지 중간 및 초음파 변환기 및 대 전자에 마커 제 선분에 초음파 변환기 홀더 및 마커 제 선분 사이의 각도의 절대 오차.
그림 5. 그림은 드롭 랜딩 작업 중에 3D 모션 캡처를 보여줍니다. 플랫폼에서 발사 초기 무릎의 굴곡에 대응이, B는 지상 눈에 띄는 발가락에 해당하고, C가 바닥을 치는 발 뒤꿈치에 대응, D는 무릎 굴곡 포스트 착륙에 해당하고 E는 KNE에 해당전자 확장 및 안정화. 큰 그림을 보려면 여기를 클릭하십시오.
그림 6. 그림은 드롭 점프의 자세 단계에서 담당자 무릎 외반 모멘트의 변화를 보여줍니다. 무릎 외반 모멘트가 착륙 단계에서 최대 초기 접지에서 0.1 ± 0.1 ㎚ / km에 0.03 ± 0.03에서 증가를 발표했다. 큰를 보려면 여기를 클릭하십시오 그림.
초음파 영상은 3 차원 모션 캡처, dynamometry, 근전도, 및 지상 반력 측정과 같은 종래의 수단을 보완 할 수 동적 연구에서 근육 운동학의 직접적인 평가를 제공 할 수있는 능력을 갖는다. 이 방법은 기본적인 생체 역학 연구 및 임상 평가에 광범위하게 적용 할 수 있습니다. 원시 무선 주파수 (RF) 초음파 데이터 또는 봉투 감지 그레이 스케일 (또는 B-모드) 이미지 데이터에 대한 상호 상관을 사용하여 (1) 얼룩 추적 방법 : 초음파를 이용하여 조직의 움직임을 추정하는 세 가지 주요 방법이 있습니다. 이러한 기술은 광범위 골격 24-25 및 심장 (26)의 근육 움직임 추적 및 추정에 모두 사용되고, (2), 근육 다발 또는 기능 27-28 및 심장 (29) 모두에 사용되는 (3) 조직 도플러 이미징 기술을 추적하는 화상 처리 방법 -30과 골격 31 모션 추정. 공간 간 C를 기반으로 얼룩 추적orrelation 조직의 움직임을 추적하고 서브 픽셀 해상도의 움직임을 추적 할 수 있습니다 널리 사용되어왔다. 그러나, 스페 클 패턴은 큰 동작시 신속하게 상관 관계를 해제. 이미지 평면 중 운동은 또한 얼룩 추적을위한 도전을 포즈. 추적 근육 다발 길이 방법은 전체 다발이 역동적 인 작업을하는 동안 이미지의 시각화 잘 적용이 있습니다. 처리 된 화상 데이터에 의존하는 방법으로는 촬상 프레임 레이트에 의해 제한 낮은 시간 해상도를 가지고 있고, 따라서 높은 속도로 동작을 추적 할 수 없다. 또한, 이러한 다발 추적 방법은 비행기의 움직임 아웃에 매우 민감하다. 따라서 근육에 프로브 이동 상대는 추적이 실패 할 수 있습니다. 기존의 조직 도플러 영상 (TDI)의 속도 추정뿐만 아니라 작은 프로브의 움직임에 더 강력하고, 더 높은 시간 해상도를 가질 수 있습니다. 도플러 방법에 따라서 도플러 추정이 부정확 한 뒤 수 있습니다 만 초음파 빔을 따라 속도의 구성 요소를 추정 할 수있다근육의 움직임 insonation의 다양한 각도로 전자. 제안한 vTDI 방법은 서로 다른 각도로 스티어링 된 두 개의 서로 다른 초음파 빔을 이용하여이 문제를 극복한다, 따라서 속도 추정치는 결상면에서 insonation 각도 독립적이다. 또한, vTDI의 유효 시간 해상도는 약 0.1 밀리 될 수 있으므로이 방법은 동적 인 활동을하는 동안 (예 : 드롭 방문, 보행과 조깅을) 골격 근육의 움직임을 추적 할 수 있습니다.
우리의 접근의 다른 이점은 벡터 조직 도플러 촬상을 수행하기위한 임상 초음파 시스템을 기반으로 선형 어레이 영상 센서의 사용을 포함한다. 우리는 전자적으로 전송이 / 넓은 시야를 스캔, 빔 스티어링, 조리개의 크기와 위치를 포커스를받을 제어. 또한,이 방식은 동시에 실시간으로 촬상 듀플렉스 vTDI을 수행하도록 확장 될 수있다. 우리의 시스템은 우리가 L로 기존의 B-모드 영상을 수행 할 수 있습니다조직의 긴장과 운동 역학의 부량 관심 영역을 ocate. 이 방법은 임상 스캐너에 구현 된 이후, 우리는 생체 역학 연구를위한 보행 실험실에서이 vTDI 방법을 배포 할 수있게되었습니다.
이 기술의 한계는 인정해야합니다. 여러 요인들은 도플러 측정의 정확도에 영향을 미친다. (을 따라 근육 섬유에서) 두 가지 차원에서 vTDI 기반의 속도 추정 (32 요소 폭) 두 개의 송신 / 수신 서브 개구로 분할, 15 °로 빔을 조종하는 선형 배열 변환기가 필요합니다. 초음파는 높은 각도로 빔을 전송 및 수신 핸들은 엽 (叶)을 격자에 의한 속도 측정에 영향을 미칠 수 있습니다. 또한 vTDI의 빔 중첩 영역의 면적은 잠재적으로 속도 추정에 영향을 미치는 32 안쪽 가변 빔 포커스를 변경한다. 도플러 추정치의 분산은 t (1) 분석 시간 창 내에서 가속 및 감속의 조직 (2) 분산에 의존도플러 레인지 게이트 내의 문제 속도 (3) 광대역에 사용되는 개구 내의 가변 도플러 각도 송신라고도 초음파 빔을 수신 스펙트럼 형상 (33) 및 송신 초음파 펄스 (4)의 대역폭을 넓혀시기 도플러 시프트 반송파 주파수 (34)에 비례한다. 몇 가지 방법은 분산을 제한하는 데 사용될 수있다. 그런 자기 상관 위상 기반 속도 추정기는, 일반적 스펙트럼 추정기에 비해 더 작은 분석 시간 윈도우를 이용하지만, 그들은 도플러 시프트보다는 피크 시프트를 의미 추정된다. 광대역 2D 같은 푸리에 스펙트럼 추정기 (35)로 인해 펄스 대역폭 분산을 줄일 수있는 변환. 이 도플러 빔 스티어링을 이용 vTDI의 경우, 근육 빔 중첩 영역 상대적 조직 속도의 변화는 또 다른 고려 사항이다. 대퇴 직근의 근육 수축은 3D 및 수축를 내보에ITY 근육을 따라 공간적으로 변화한다. 따라서주의 깊게 관심 영역을 선택하는 것이 중요합니다.
이 연구에서, 우리는 vTDI를 사용하여 여덟 건강한 지원자에서 드롭 랜딩 작업하는 동안 대퇴 직근 근육 운동학의 반복성을 조사 하였다. 시험은 독립했다하더라도, 우리는 시련과 개인에 대한 높은 상관 반복 피크 근육의 수축 속도를 관찰했다. 우리는 현재 추가로이 패턴을 조사하는 연구에 더 많은 과목을 모집하고 있습니다. 이 연구는 비 침습적 제공하고 직근의 수축 속도의 실시간 측정 드롭 착륙하는 동안 근육 대퇴. 수축 속도의 다음과 같은 패턴을 드롭 랜딩 작업 (그림 2)의 다양한 단계에서 관찰되었다 : 1. 근육의 수축 속도는 (무릎 굴곡 (실행 단계) 및 연장 동안 축 방향에 비해 좌우 방향으로 지배적 인 공기 P하세). 대퇴 직근 근육의 공중파 단계에서 실행 단계 및 동심 수축시 편심 수축을 겪고 있기 때문에 이것은 예상된다. 2. 무시할 낮은 축 근육 속도와 세 번째 단계 (발가락이지면에 닿지) 동안 낮은 측면 근육 속도. 이것은이 단계 3시 대퇴 직근의 근육 수축을 낮출에 해당합니다. 발 뒤꿈치가 땅에 닿을 직후 축 방향 및 횡 방향의 근육 속도의 상당한 증가. 이 근육 섬유 각각 근육 섬유에 정상에 따라 속도의 증가를 초래, 아마 때문에 압축 형태로 편심 수축과 변화 모두를 겪고 근육 때문이다. 드롭 랜딩 작업이 높은 충격 작업이라는 사실에도 불구하고, vTDI는 반복 직근 근육 속도를 대퇴 보여 주었다. 이 근육이 과도한로드에서 무릎 관절을 보호 할 책임이있는 때문에이 초음파 기술은 임상 적 영향을 미칠 수 있습니다.따라서, 전방 십자 인대 재건술 환자에서 대퇴 직근 근육의 추가 평가는 OA의 초기 및 가속 발병으로 이어지는 메커니즘을 이해하기 위해 보증합니다.
이 연구의 참가자는 모두 30 ㎝ 높이 플랫폼에서 자연 드롭 랜딩 작업을 수행하도록 요청했지만, 우리는 점프 또는 실행의 높이의 차이를 발견했다. 또한, 고속 카메라 데이터를 사용하여, 모든 피사체가 다른 방울 착지 스타일을 한 것으로 관찰되었다. 이 작업을하는 동안 활성화 패턴의 가능한 차이의 결과로 근육 대퇴 직근의 피크 결과 속도 값의 주제 사이의 약간의 차이를 설명 할 수있다. 또 다른 가능한 요인은 잠재적으로 근육 수축의 다른 수준으로 이어질 및 생산을 강제 할 수있는 대퇴 직근 근육의 단면적의 차이입니다.
저자의 아무도는 재정적 공개 또는 관심과 연구의 충돌이 우리 기관의 IRB 승인을 받았다가 없습니다.
이 작품은 국립 과학 재단 (National Science Foundation)에서 부여 번호 0953652에 의해 부분적으로는 조지 메이슨 대학 도서관 오픈 액세스 출판 기금에 의해 부분적으로 지원되었다. 우리는 고속 카메라에 대한 액세스를 제공하는 닥터 존 로버트 Cressman 주니어에게 감사의 말씀을 전합니다.
Name | Company | Catalog Number | Comments |
Ultrasound System | Ultrasonix | Sonix RP | |
3D Motion Capture System | Vicon Motion Systems | Vicon T-20 | |
Force Plates | Bertec Corporation | Bertec 4060-10 | |
High Speed Camera | Photron | Photron 512 PCI 32K |
JoVE'article의 텍스트 или 그림을 다시 사용하시려면 허가 살펴보기
허가 살펴보기This article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. 판권 소유