Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1702 1098 (N16W44 L=284)  5266    988        5 (1416/ 777-1894/1408) -27.8 W19-W75 N00-N30
      2  220 1274 (N00E64 L=175)   444    192        0 ( 174/1242- 270/1324) -39.8 E70-E59 S02-N04
      3  338  913 (S19E45 L=195)   312     90        1 ( 321/ 872- 355/ 954)  90.0 E45-E42 S22-S16
      4  449  914 (S17E35 L=204)   259     86        0 ( 429/ 896- 487/ 940) -11.3 E37-E32 S18-S16
      5 1758  689 (S07W59 L=299)  1241    281        1 (1685/ 658-1842/ 724)  12.5 W51-W70 S09-S04
      6  578  642 (S33E24 L=216)   303    119        2 ( 530/ 621- 647/ 672)  -6.3 E27-E18 S33-S31
      7  960  638 (S25W03 L=243)   429    112        2 ( 894/ 614-1017/ 665)   7.1 E01-W07 S26-S23
      8 1477  570 (S20W39 L=279)   532    133        1 (1412/ 552-1531/ 591)   6.3 W34-W43 S20-S18
      9 1407  337 (S37W46 L=286)   375     58        2 (1385/ 312-1427/ 362) -53.1 W44-W47 S38-S34
     10 1317  300 (S41W42 L=282)   355     40        0 (1296/ 283-1337/ 320)  76.0 W39-W43 S43-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