Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1724 1343 (N19W54 L= 34)   685     86        4 (1694/1306-1762/1382) -33.7 W51-W58 N16-N21
      2  214 1246 (N20E63 L=276)   720    124        1 ( 191/1200- 247/1301) -31.0 E68-E58 N17-N23
      3 1611 1227 (N13W41 L= 20)   269     60        0 (1572/1214-1647/1241) -11.3 W37-W44 N12-N13
      4  411 1069 (N09E41 L=299)   279    105        0 ( 397/1021- 439/1106)  68.2 E41-E38 N06-N11
      5  569  999 (N05E29 L=310)   703    186        1 ( 543/ 914- 595/1096)  84.3 E31-E27 N00-N10
      6  290  872 (S03E52 L=287)  3423    615        4 ( 213/ 662- 490/1081) -50.2 E60-E37 S15-N10
      7 1015  864 (S05E01 L=338)   275     86        0 ( 975/ 843-1053/ 896)  26.6 E03-W01 S06-S03
      8 1237  733 (S15W12 L=351)   401    130        2 (1166/ 695-1276/ 761) -39.8 W07-W14 S17-S12
      9  457  581 (S21E44 L=296)   769    175        0 ( 379/ 545- 533/ 624) -22.6 E49-E37 S23-S18
     10  711  429 (S33E27 L=312)   357     18        0 ( 688/ 410- 740/ 448)  11.3 E29-E24 S34-S31
     11  617  301 (S43E43 L=296)  1207    118        0 ( 577/ 256- 653/ 350)  26.6 E54-E36 S48-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