Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  238 1319 (N41E61 L= 42)   440     65        2 ( 207/1303- 271/1331)  -9.5 E67-E55 N39-N41
      2  369 1286 (N36E43 L= 60)   308     81        0 ( 334/1263- 413/1301)   0.0 E46-E38 N35-N37
      3  389 1102 (N25E41 L= 62)   430    122        1 ( 357/1057- 415/1164) -74.1 E43-E38 N21-N28
      4 1712  994 (S12W43 L=146)   341     53        0 (1678/ 981-1740/1010) -21.8 W40-W45 S13-S11
      5  165  889 (N14E69 L= 34)  1422    178        4 ( 140/ 820- 203/ 982) -78.7 E72-E64 N10-N20
      6  423  838 (N08E43 L= 60)   446     80        0 ( 398/ 808- 451/ 867)  31.0 E44-E39 N06-N09
      7 1206  627 (S20W01 L=104)  1601    509        1 (1031/ 470-1348/ 790)  26.6 E13-W13 S26-S13
      8 1417  583 (S29W14 L=117)   251     84        0 (1378/ 553-1449/ 616)  18.4 W11-W17 S30-S27
      9  904  452 (S25E22 L= 81)  1550    402        2 ( 845/ 377-1014/ 610)  56.3 E26-E15 S31-S14
     10 1252  433 (S35E00 L=102)   272     63        1 (1222/ 409-1281/ 460)  14.0 E03-W01 S36-S33
     11  594  264 (S32E59 L= 44)   932     92        2 ( 568/ 218- 619/ 305)  41.6 E64-E54 S36-S28
  

    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