Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1371 (N33W29 L=124)   415     80        1 (1581/1333-1621/1414) -51.3 W26-W30 N30-N35
      2 1791 1121 (N23W49 L=145)  2013    372        5 (1668/ 992-1890/1196)  38.0 W40-W63 N12-N30
      3  367  767 (S36E40 L= 56)  1782    256        7 ( 292/ 678- 445/ 874)  46.8 E50-E33 S44-S28
      4 1740  662 (S04W59 L=155)   479     78        0 (1729/ 619-1753/ 706) -45.0 W56-W61 S06-S01
      5 1354  619 (S18W32 L=127)   510    190        0 (1283/ 581-1426/ 685)  58.0 W28-W33 S21-S12
      6 1128  621 (S24W18 L=114)   292     76        0 (1096/ 600-1162/ 640)  11.3 W15-W20 S25-S22
      7 1433  579 (S18W39 L=134)   487    135        1 (1369/ 560-1492/ 608)  45.0 W34-W42 S20-S14
      8  994  540 (S33W13 L=109)   588    199        2 ( 920/ 482-1076/ 604)  65.6 W09-W16 S38-S27
  

    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