Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  535 1701 (N24E58 L= 23)   508     84        3 ( 504/1670- 558/1731)  90.0 E60-E55 N20-N26
      2  407 1608 (N15E63 L= 18)   620     80        1 ( 387/1571- 435/1645) -14.0 E69-E57 N12-N17
      3  278 1519 (N07E72 L=  8)  1197     89        2 ( 249/1479- 310/1548)  68.2 E77-E68 N04-N09
      4 1887 1035 (N22W63 L=144)  2867    281        4 (1815/ 913-1947/1134)  38.0 W56-W80 N12-N30
      5  534  708 (S36E24 L= 56)  2272    399        7 ( 435/ 600- 658/ 846)  54.0 E33-E17 S45-S23
      6 1570  611 (S12W47 L=127)  1379    386        2 (1452/ 489-1648/ 698)  77.5 W39-W52 S22-S04
      7 1254  494 (S28W31 L=112)  1033    250        1 (1143/ 433-1356/ 546)  53.1 W27-W36 S34-S22
      8 1556  483 (S19W53 L=134)   470    126        1 (1507/ 440-1607/ 516)  66.8 W51-W56 S22-S15
      9 1101  426 (S37W25 L=106)   522     87        2 (1070/ 391-1133/ 463) -59.0 W22-W28 S39-S34
  

    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