Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  398 1581 (N31E60 L=178)  1351    151        0 ( 343/1537- 463/1627)  -6.8 E76-E51 N28-N33
      2 1683 1431 (N33W50 L=288)  1036    161        1 (1645/1363-1715/1492)  71.6 W46-W53 N28-N37
      3  567 1253 (N12E32 L=206)   252     74        0 ( 545/1217- 600/1286) -38.7 E34-E29 N09-N13
      4 1652  802 (S07W45 L=283)   979    197        3 (1616/ 731-1673/ 896) -68.2 W41-W47 S11-S01
      5  555  833 (S14E29 L=209)   641    170        0 ( 481/ 807- 637/ 860)  -5.2 E34-E23 S15-S12
      6  867  668 (S22E07 L=231)   229    122        0 ( 849/ 621- 891/ 712) -80.5 E08-E05 S24-S18
      7 1171  584 (S25W14 L=252)   250    102        0 (1122/ 570-1216/ 606)  -8.1 W10-W17 S25-S23
      8 1377  522 (S27W30 L=268)   289     32        0 (1352/ 508-1402/ 537)  33.7 W28-W32 S28-S26
      9  395  463 (S43E53 L=185)   594     66        0 ( 378/ 425- 415/ 494)  32.5 E59-E47 S46-S39
  

    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