Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  670 1683 (N58E19 L=  2)  1995    220        4 ( 586/1615- 778/1726)   0.0 E30-E05 N52-N61
      2 1226 1641 (N41W31 L= 51)  2089    454        4 (1022/1608-1434/1667) -22.8 W14-W46 N34-N47
      3 1612 1519 (N24W55 L= 76)   363     51        0 (1593/1490-1636/1546) -76.0 W53-W56 N21-N25
      4  243 1334 (N36E63 L=318)   960    134        0 ( 208/1272- 287/1369)  36.9 E66-E56 N32-N38
      5  464 1309 (N33E36 L=345)   289     37        0 ( 434/1294- 493/1325)   0.0 E38-E33 N31-N33
      6 1157 1212 (N15W11 L= 32)   818    237        2 (1086/1133-1203/1299) -70.0 W08-W13 N10-N21
      7 1578  679 (S25W31 L= 52)  1570    197        1 (1545/ 605-1634/ 765) -63.4 W29-W36 S31-S19
  

    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