Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  973 1613 (N33E19 L=296)   254     63        0 ( 944/1590-1005/1641)  68.2 E20-E17 N30-N35
      2 1429 1370 (N27W18 L=333)   745    148        2 (1364/1326-1482/1406)  -6.3 W13-W22 N25-N29
      3  340  977 (S19E45 L=270)   296     95        2 ( 306/ 949- 373/1012)  45.0 E47-E42 S21-S16
      4 1677  862 (N03W46 L=  0)   369    113        1 (1619/ 842-1719/ 884)  35.5 W41-W48 N00-N05
      5  639  819 (S22E20 L=295)   380     90        1 ( 589/ 806- 686/ 837)  15.9 E23-E16 S23-S20
      6 1004  614 (S27W08 L=323)   242     94        1 ( 965/ 595-1046/ 633)   9.5 W05-W11 S27-S25
      7 1604  471 (S21W57 L= 12)  1027    232        0 (1507/ 423-1698/ 516)  37.6 W51-W64 S26-S16
      8 1667  455 (S20W66 L= 20)   531     57        1 (1640/ 434-1693/ 482)  45.0 W63-W68 S21-S18
  

    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