Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  480 1559 (N52E37 L= 73)   550     72        1 ( 460/1517- 502/1593) -71.6 E40-E34 N49-N55
      2  710 1559 (N47E12 L= 99)   698    111        1 ( 653/1535- 764/1577) -26.6 E16-E06 N44-N49
      3 1784 1450 (N12W70 L=180)   750     93        1 (1768/1406-1801/1486)  45.0 W66-W76 N09-N14
      4  663 1273 (N29E19 L= 92)   216     70        0 ( 645/1238- 687/1299)  45.0 E20-E16 N26-N30
      5  308 1126 (N25E50 L= 61)   865    206        3 ( 239/1051- 360/1203) -37.6 E58-E44 N19-N30
      6  202 1138 (N27E63 L= 48)   944    168        1 ( 183/1061- 236/1197)  76.0 E65-E58 N22-N30
      7  133 1056 (N22E74 L= 36)   661    125        0 ( 116/ 999- 154/1114)  49.4 E78-E70 N18-N25
      8  781  815 (S01E19 L= 92)   559    197        0 ( 734/ 730- 832/ 890) -77.5 E20-E16 S06-N04
      9 1047  732 (S11E05 L=106)   226     89        1 (1012/ 695-1075/ 760)  31.0 E07-E02 S12-S08
     10 1702  705 (S29W41 L=152)   411     47        1 (1682/ 674-1720/ 736)  90.0 W39-W42 S30-S26
     11  448  387 (S23E61 L= 50)  1989    218        5 ( 368/ 316- 520/ 439) -55.0 E68-E54 S29-S18
     12 1131  373 (S37E08 L=102)   266     69        1 (1103/ 352-1165/ 399) -53.1 E10-E06 S39-S34
  

    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