Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1596 1573 (N45W25 L=142)   377     39        1 (1568/1550-1622/1601) -21.8 W21-W27 N42-N46
      2  902 1477 (N15E19 L= 98)   513    158        0 ( 875/1404- 935/1550) -45.0 E23-E15 N11-N19
      3 1033 1445 (N17E11 L=106)   314     87        0 ( 992/1431-1083/1456)  18.4 E13-E07 N15-N18
      4 1125 1348 (N13E03 L=114)   248     91        0 (1105/1300-1148/1390) -45.0 E04-W00 N11-N15
      5 1814 1101 (N21W50 L=167)   793    154        0 (1799/1023-1830/1177) -66.0 W48-W53 N16-N25
      6 1576  889 (N01W36 L=153)   422     78        0 (1542/ 845-1607/ 917)   0.0 W33-W39 N00-N02
      7  850  594 (S36W02 L=119)  1133    281        1 ( 729/ 561- 983/ 641)  38.2 E06-W09 S41-S29
      8 1045  485 (S37W19 L=136)   517    115        1 (1025/ 435-1070/ 538) -74.1 W16-W22 S40-S33
      9 1479  441 (S25W50 L=167)   637    153        0 (1418/ 385-1524/ 502)  90.0 W48-W51 S30-S20
  

    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