Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1187 1824 (N52W25 L= 64)   375     81        2 (1170/1793-1208/1850)  39.8 W22-W28 N49-N54
      2 1485 1666 (N36W41 L= 80)   443    107        1 (1443/1646-1528/1689)  -9.5 W38-W44 N34-N38
      3 1639 1442 (N20W46 L= 85)   349    135        2 (1593/1416-1685/1455)  -8.1 W41-W49 N18-N20
      4  792 1372 (N21E12 L= 27)   305     82        1 ( 753/1348- 840/1393) -31.0 E14-E09 N19-N22
      5 1742 1358 (N14W53 L= 92)   438    133        3 (1694/1324-1786/1384) -20.6 W48-W56 N11-N15
      6  377 1073 (N06E42 L=357)   326    120        2 ( 366/1023- 398/1148) -90.0 E42-E40 N02-N10
      7  275 1047 (N05E51 L=348)  1300    309        2 ( 190/ 917- 315/1143) -52.4 E59-E47 S02-N11
      8  145 1063 (N07E65 L=334)   828    210        7 ( 106/1025- 185/1093)   0.0 E71-E59 N05-N09
      9 1769 1019 (S06W50 L= 89)   305     62        0 (1752/ 984-1785/1055)  63.4 W48-W52 S08-S04
     10  224  967 (N01E56 L=342)   740    245        5 ( 187/ 882- 258/1039) -60.9 E59-E53 S04-N05
     11  141  955 (N01E67 L=332)  2270    285        6 (  93/ 904- 185/ 998) -26.6 E75-E61 S02-N03
     12 1863  789 (S21W61 L=100)   598    158        0 (1826/ 757-1907/ 829)  24.0 W57-W66 S22-S18
     13  196  767 (S11E63 L=335)   567    120        2 ( 175/ 728- 215/ 824)  35.5 E67-E59 S12-S07
     14 1293  611 (S29W15 L= 54)   207     73        0 (1250/ 592-1330/ 626)  11.3 W12-W17 S29-S28
     15  294  579 (S23E58 L=340)   786    148        4 ( 270/ 544- 319/ 621)  78.7 E61-E55 S25-S20
     16  806  365 (S42E24 L= 15)   607    133        0 ( 743/ 343- 880/ 388)   0.0 E29-E17 S44-S41
  

    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