Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  985 1921 (N64W48 L= 37)   534     49        0 ( 962/1894-1002/1942)  19.7 W41-W55 N60-N66
      2  842 1857 (N64W22 L= 11)  1273    156        2 ( 787/1815- 906/1904)   9.1 W12-W37 N60-N67
      3  717 1782 (N62W02 L=351)   468     42        0 ( 702/1751- 736/1811)  50.2 E01-W04 N59-N65
      4  935 1627 (N42W12 L=  1)  1457    329        0 ( 788/1597-1094/1659) -20.0 W00-W22 N35-N45
      5 1692 1555 (N16W62 L= 52)   685     53        2 (1661/1535-1725/1574)   0.0 W58-W67 N14-N17
      6 1236 1512 (N26W26 L= 16)   383     60        0 (1217/1471-1260/1546)  45.0 W24-W28 N23-N28
      7  398 1425 (N43E36 L=313)   454     44        0 ( 369/1402- 435/1443)  -8.1 E39-E32 N41-N44
      8 1531 1316 (N08W37 L= 27)   723    187        0 (1493/1231-1567/1395)  78.7 W35-W39 N02-N12
      9 1809 1244 (S04W57 L= 46)   366     39        1 (1788/1221-1829/1266)  14.0 W54-W60 S05-S02
     10 1205 1108 (N04W12 L=  1)   442    134        0 (1166/1052-1251/1165) -76.0 W10-W13 N00-N08
     11  161 1033 (N23E63 L=286)  5918    302        6 (  92/ 935- 237/1138) -47.3 E76-E53 N16-N29
     12  283 1081 (N24E48 L=301)   385     48        0 ( 265/1050- 302/1112) -63.4 E50-E46 N21-N25
     13  426  921 (N11E39 L=310)   900    207        3 ( 349/ 877- 509/ 998) -45.0 E43-E33 N07-N17
     14 1203  806 (S12W05 L=354)   524    158        0 (1177/ 736-1232/ 888)  45.0 W01-W08 S15-S08
     15  154  748 (N06E72 L=277)  2266    117        6 ( 121/ 690- 193/ 810) -82.9 E81-E66 N02-N10
     16 1740  601 (S41W43 L= 32)   786    126        0 (1681/ 579-1805/ 623) -21.3 W35-W53 S45-S37
     17 1590  517 (S42W27 L= 16)   477     74        0 (1558/ 487-1626/ 549)  23.2 W23-W30 S43-S40
  

    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