Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1363 1784 (N46W58 L=235)   446     53        1 (1336/1765-1389/1799)  -8.1 W53-W62 N44-N47
      2 1607 1650 (N30W66 L=244)   852     87        2 (1576/1625-1646/1673) -10.3 W62-W73 N28-N31
      3 1680 1582 (N24W68 L=246)  1048    125        3 (1640/1526-1719/1639)  53.1 W64-W71 N20-N28
      4 1608 1557 (N25W58 L=235)  1272    194        5 (1515/1510-1679/1590) -36.0 W49-W64 N20-N29
      5 1205 1404 (N26W20 L=197)  1125    229        4 (1118/1327-1279/1481) -69.0 W16-W22 N19-N32
      6 1033 1441 (N32W09 L=187)   226    129        0 ( 997/1386-1078/1473)  20.6 W05-W13 N29-N33
      7 1400 1424 (N23W33 L=211)   279     88        1 (1354/1407-1435/1448)  -9.5 W30-W36 N21-N23
      8  792 1143 (N18E12 L=165)   263    129        1 ( 760/1098- 826/1196) -74.1 E13-E10 N14-N21
      9 1044 1116 (N12W03 L=180)   298    186        1 ( 980/1063-1127/1137)  -6.3 E01-W08 N09-N13
     10 1057 1040 (N07W02 L=180)   185     96        0 (1011/1025-1086/1060) -26.6 W00-W04 N06-N08
     11  133  930 (N13E75 L=102)  1380    129        3 ( 118/ 865- 149/ 982)  74.1 E80-E71 N08-N15
     12 1398  901 (S08W20 L=198)   491    180        1 (1380/ 820-1446/ 953)  54.5 W18-W24 S12-S04
     13 1148  867 (S05W04 L=182)   625    163        3 (1093/ 839-1239/ 903) -32.0 W01-W09 S07-S02
     14 1274  742 (S15W10 L=187)   518    169        1 (1233/ 677-1297/ 818)  54.5 W06-W11 S18-S10
     15 1422  624 (S26W18 L=196)  1554    478        0 (1274/ 468-1584/ 768)  26.6 W05-W31 S33-S19
     16 1590  632 (S30W31 L=209)   880    269        1 (1510/ 530-1662/ 719)  37.6 W24-W38 S35-S25
     17  447  538 (S14E50 L=128)   413     72        2 ( 412/ 509- 482/ 557) -18.4 E52-E46 S16-S12
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y