Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  979 1774 (N58W14 L= 60)   397     54        1 ( 943/1753-1008/1790)  -7.1 W10-W18 N56-N59
      2  877 1757 (N59W03 L= 48)   352     74        1 ( 852/1729- 903/1777) -45.0 W00-W06 N55-N60
      3  753 1732 (N59E12 L= 33)   530    138        0 ( 685/1711- 814/1750)  -3.6 E20-E04 N57-N60
      4 1223 1491 (N32W22 L= 67)   352    161        0 (1149/1439-1273/1516) -49.4 W17-W23 N27-N34
      5  259 1334 (N34E59 L=346)  3425    420        5 ( 171/1139- 350/1492)  68.6 E66-E52 N21-N44
      6 1030 1418 (N30W07 L= 52)   384    159        0 ( 976/1386-1107/1460) -36.9 W03-W11 N26-N33
      7 1718 1319 (N12W53 L= 98)   859    187        1 (1687/1268-1778/1395) -90.0 W50-W58 N08-N16
      8 1865  953 (S13W61 L=106)   801    101        2 (1836/ 913-1894/1003) -50.2 W57-W65 S15-S09
      9 1177  936 (S01W08 L= 53)   184    101        0 (1134/ 922-1225/ 953) -21.8 W05-W10 S02-N00
     10 1429  794 (S14W22 L= 67)   571    163        1 (1369/ 744-1489/ 844)  24.0 W17-W26 S15-S11
     11  578  633 (S12E35 L= 11)   391     99        0 ( 561/ 585- 593/ 681)  74.1 E36-E33 S15-S08
     12  287  605 (S12E62 L=343)   609     65        1 ( 270/ 572- 309/ 642) -90.0 E63-E59 S14-S09
     13  387  602 (S13E51 L=354)   458     54        0 ( 357/ 578- 424/ 626) -36.9 E53-E48 S14-S10
     14  515  536 (S18E42 L=  3)   477    109        0 ( 461/ 507- 557/ 576)  20.0 E48-E37 S19-S15
     15 1423  400 (S42W20 L= 66)   423     53        2 (1392/ 376-1457/ 419)   0.0 W17-W23 S43-S39
     16 1253  308 (S46W05 L= 50)   799    101        2 (1204/ 281-1305/ 337)   5.7 W00-W10 S48-S43
     17 1081  255 (S48E11 L= 34)   320     35        0 (1062/ 228-1097/ 276) -78.7 E12-E09 S50-S45
     18 1046  202 (S53E17 L= 28)   913     84        1 (1012/ 169-1085/ 237) -26.6 E20-E12 S56-S48
  

    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