Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  698 1561 (N48E06 L= 42)   412     64        1 ( 660/1533- 733/1585)  -9.5 E09-E03 N46-N49
      2  313 1459 (N52E52 L=356)   674    121        1 ( 273/1412- 365/1498)  20.0 E57-E45 N49-N54
      3  366  920 (N16E45 L=  4)  2800    503        5 ( 283/ 770- 513/1138) -65.2 E50-E36 N04-N30
      4  747 1007 (N13E16 L= 32)   237     52        1 ( 715/ 985- 767/1023) -45.0 E18-E14 N11-N14
      5  180  754 (N09E71 L=338)   807    102        1 ( 153/ 707- 203/ 799)  21.0 E79-E65 N06-N11
      6 1237  431 (S35E04 L= 44)   889    206        0 (1167/ 373-1286/ 515)  28.6 E10-W01 S37-S30
      7  521  467 (S13E49 L=360)   329     61        1 ( 500/ 436- 539/ 498)  45.0 E51-E46 S15-S11
      8  587  392 (S20E48 L=  0)   584    108        2 ( 567/ 355- 609/ 458)  60.3 E51-E44 S22-S15
      9  981  214 (S44E34 L= 15)  1120    175        2 ( 918/ 169-1037/ 274)  18.4 E43-E25 S46-S40
  

    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