Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  904 1792 (N37E35 L= 25)  2874    716        6 ( 838/1625- 949/1939) -37.7 E54-E23 N24-N50
      2  504 1479 (N06E45 L= 15)   765    234        1 ( 452/1422- 567/1549) -26.6 E51-E38 N03-N09
      3 1786 1386 (N39W43 L=103)  1481    221        2 (1720/1325-1853/1443)  -8.1 W36-W51 N36-N41
      4  817 1400 (N09E22 L= 38)   219    103        1 ( 765/1365- 854/1426)  -8.1 E25-E18 N07-N09
      5  336 1372 (S03E54 L=  5)   429    113        2 ( 309/1345- 366/1406) -23.2 E57-E50 S04-S01
      6  366 1279 (S08E48 L= 11)   499     96        2 ( 336/1242- 391/1331)  80.5 E50-E46 S10-S04
      7  442 1288 (S06E42 L= 17)   245     55        2 ( 415/1268- 470/1306)  45.0 E44-E40 S07-S04
      8  762 1175 (S05E19 L= 41)   784    323        1 ( 637/1063- 884/1300)  64.8 E23-E14 S14-N04
      9  385 1078 (S19E41 L= 18)   269     57        1 ( 360/1061- 412/1100)   0.0 E43-E38 S20-S17
     10  333  990 (S26E44 L= 16)   864    149        3 ( 281/ 951- 372/1037) -33.7 E48-E40 S27-S22
  

    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