Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  981 1919 (N64W47 L= 37)   502     46        0 ( 958/1891- 998/1941)  19.7 W40-W54 N60-N66
      2  839 1857 (N64W22 L= 12)  1021    111        2 ( 800/1811- 894/1899)   6.0 W15-W34 N60-N67
      3  708 1777 (N62W00 L=350)   443     31        0 ( 690/1748- 724/1804)  68.2 E01-W02 N58-N64
      4  930 1624 (N42W12 L=  2)  1398    304        0 ( 796/1593-1088/1656) -21.8 W01-W21 N35-N45
      5 1685 1552 (N16W61 L= 51)   650     62        1 (1650/1532-1721/1571)   0.0 W57-W66 N14-N17
      6 1226 1507 (N26W26 L= 16)   368     66        0 (1204/1467-1250/1540)  45.0 W23-W28 N23-N27
      7  390 1422 (N43E37 L=313)   440     53        1 ( 361/1400- 427/1440) -15.9 E40-E33 N41-N44
      8 1521 1312 (N07W37 L= 26)   670    163        1 (1480/1235-1555/1389) -90.0 W35-W38 N02-N12
      9 1804 1240 (S04W57 L= 47)   382     44        1 (1782/1217-1827/1263)  11.3 W54-W59 S05-S02
     10 1197 1099 (N04W11 L=  1)   454    120        0 (1155/1046-1240/1154) -74.1 W09-W12 N00-N07
     11  155 1031 (N23E64 L=286)  5782    300       11 (  88/ 929- 231/1142) -65.2 E77-E54 N16-N29
     12  386  954 (N14E41 L=309)  1320    338        3 ( 260/ 869- 495/1091) -48.6 E50-E34 N07-N24
     13 1194  803 (S12W04 L=354)   494    157        0 (1168/ 732-1223/ 881)  49.4 W01-W07 S15-S08
     14  149  741 (N06E73 L=277)  1641    113        2 ( 122/ 697- 179/ 804) -90.0 E79-E69 N03-N09
     15 1752  597 (S42W45 L= 34)   995    164        1 (1680/ 575-1829/ 619) -15.3 W35-W57 S45-S38
     16  568  520 (S14E41 L=309)   271     60        0 ( 555/ 488- 580/ 552)  53.1 E42-E39 S16-S12
     17 1586  516 (S42W26 L= 16)   527     86        0 (1550/ 485-1628/ 548)   7.1 W22-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