Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  606 1674 (N58E28 L=  3)  1559    210        5 ( 533/1626- 688/1720)  13.2 E38-E17 N54-N61
      2 1109 1604 (N41W20 L= 50)  1458    398        4 ( 936/1573-1290/1641) -20.3 W06-W33 N36-N47
      3 1253 1379 (N24W21 L= 52)   224    135        0 (1196/1362-1309/1405) -23.2 W18-W25 N22-N26
      4  989 1203 (N18W01 L= 32)   394    140        0 ( 932/1153-1030/1253) -66.8 E02-W02 N14-N21
      5 1879  812 (S25W65 L= 96)  1016    194        1 (1829/ 748-1923/ 901)  28.3 W57-W73 S27-S20
      6 1398  663 (S22W18 L= 49)  1400    302        6 (1333/ 553-1475/ 790) -74.5 W13-W22 S31-S13
      7 1711  647 (S31W45 L= 75)   364     47        2 (1684/ 629-1733/ 667)  14.0 W41-W47 S32-S29
  

    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