Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1113 1772 (N52E09 L=163)   572    102        1 (1066/1748-1169/1788)  24.0 E12-E03 N49-N54
      2  538 1752 (N39E62 L=110)   634     92        0 ( 498/1729- 581/1767)  29.7 E66-E56 N36-N40
      3  753 1563 (N29E28 L=143)  2200    407        2 ( 591/1462- 924/1646) -12.7 E44-E13 N25-N32
      4 1631 1561 (N45W44 L=215)   797     93        0 (1575/1539-1676/1581)  12.1 W36-W50 N42-N46
      5  940 1292 (N15E09 L=162)   432    151        0 ( 887/1226- 982/1336) -32.0 E13-E05 N11-N16
      6  313 1293 (N06E53 L=118)   504     81        1 ( 293/1248- 341/1334) -45.0 E56-E50 N03-N08
      7 1594 1180 (N17W35 L=207)   252     52        0 (1561/1168-1626/1193)  21.8 W32-W37 N16-N18
      8 1743 1165 (N19W48 L=220)   329     56        0 (1716/1137-1767/1187) -21.8 W45-W50 N17-N20
      9 1900 1023 (N12W68 L=239)  5051    227        5 (1866/ 904-1942/1118) -77.0 W62-W76 N05-N18
     10 1808  936 (N06W57 L=228)   807    124        0 (1792/ 863-1822/ 998) -63.4 W54-W59 N01-N09
     11  976  855 (S11E01 L=171)   256     91        0 ( 950/ 821-1013/ 887)  68.2 E01-W01 S13-S08
     12 1031  661 (S22W06 L=177)  1810    514        0 ( 798/ 617-1207/ 760)  31.3 E08-W15 S27-S13
     13 1141  226 (S53W31 L=203)   788    111        1 (1103/ 194-1173/ 263) -15.3 W26-W37 S55-S49
  

    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