Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  910 1891 (N65W34 L= 10)  1009     72        1 ( 869/1854- 948/1928)   6.0 W26-W46 N61-N67
      2  826 1833 (N62W17 L=353)   351     46        0 ( 806/1812- 851/1856) -90.0 W14-W20 N59-N65
      3 1107 1689 (N41W27 L=  3)  1632    386        3 ( 949/1650-1285/1729) -26.6 W14-W37 N32-N45
      4  723 1575 (N44E05 L=330)   281     61        1 ( 700/1552- 745/1601)  26.6 E07-E02 N42-N45
      5 1660 1442 (N11W52 L= 28)   556     92        0 (1631/1393-1688/1483) -90.0 W50-W53 N07-N13
      6 1404 1209 (N05W26 L=  2)   445     91        0 (1367/1159-1439/1250) -80.5 W24-W27 N01-N07
      7  440 1130 (N23E34 L=302)   281     54        1 ( 410/1107- 468/1156)   0.0 E36-E31 N21-N24
      8  197 1064 (N24E58 L=277)  2077    214        5 ( 149/ 996- 243/1141) -55.0 E64-E52 N19-N29
      9 1759 1083 (S11W48 L= 23)   413     72        0 (1732/1046-1783/1121) -68.2 W46-W49 S13-S08
     10 1413  941 (S10W20 L=356)  1119    319        0 (1355/ 798-1446/1086)  51.8 W13-W24 S16-S02
     11  613  989 (N11E25 L=311)   896    253        2 ( 515/ 929- 713/1077) -52.4 E29-E19 N05-N18
     12  269 1011 (N19E51 L=285)   585    102        1 ( 250/ 960- 300/1058) -90.0 E52-E47 N16-N22
     13  198  946 (N17E59 L=276)   537     80        1 ( 165/ 925- 239/ 970) -18.4 E63-E54 N15-N18
     14  215  750 (N05E63 L=273)  3132    305        6 ( 147/ 653- 266/ 857)  36.3 E75-E54 N00-N11
     15 1834  649 (S41W56 L= 32)  1045    116        0 (1784/ 625-1890/ 674) -16.3 W47-W71 S44-S37
     16 1683  556 (S42W37 L= 12)  1171    250        0 (1586/ 479-1761/ 641)  18.4 W26-W45 S45-S38
     17 1067  510 (S27E11 L=325)   290     61        0 (1039/ 477-1091/ 541)  31.0 E13-E08 S28-S24
     18  646  219 (S36E59 L=277)   831     85        1 ( 610/ 188- 675/ 261)  13.2 E68-E51 S38-S33
  

    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