Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  930 1664 (N34E25 L=270)  1047    179        2 ( 851/1613- 984/1709)  -8.1 E33-E19 N30-N35
      2  551 1658 (N24E53 L=242)   502     73        2 ( 533/1620- 572/1691) -63.4 E55-E50 N20-N26
      3  681 1596 (N23E38 L=256)   894    266        0 ( 605/1502- 732/1677) -27.9 E49-E32 N16-N26
      4 1486 1588 (N45W19 L=313)   340     61        1 (1458/1561-1516/1605)  21.8 W16-W21 N42-N47
      5  953 1391 (N16E13 L=281)   225     49        1 ( 932/1364- 977/1413)  71.6 E14-E11 N14-N17
      6  210 1207 (S10E63 L=231)  2560    207        5 ( 168/1127- 262/1282) -41.6 E70-E56 S14-S06
      7  138 1150 (S15E74 L=221)  1988    184        7 ( 112/1085- 167/1240)  66.0 E79-E68 S18-S09
      8 1878 1153 (N27W63 L=358)   946     65        2 (1843/1120-1904/1192)  59.0 W57-W68 N24-N30
      9 1488 1101 (N13W26 L=321)   313     86        0 (1463/1050-1509/1139) -38.7 W23-W28 N10-N14
     10 1411  504 (S23W40 L=335)   414     81        0 (1398/ 462-1424/ 553) -51.3 W37-W43 S25-S20
     11 1506  472 (S22W49 L=344)   449    102        0 (1485/ 435-1543/ 512) -80.5 W47-W51 S25-S19
     12 1048  228 (S53W36 L=331)   477     54        0 (1031/ 206-1071/ 259) -26.6 W32-W40 S55-S50
  

    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