Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1452 1216 (N13W22 L= 12)   244     91        0 (1425/1171-1465/1250) -36.9 W19-W23 N10-N14
      2  218 1086 (S14E57 L=292)  2527    463        8 ( 147/1000- 277/1178)  65.6 E65-E51 S19-S08
      3 1885  881 (N04W65 L= 54)  1339    228        4 (1831/ 801-1920/ 931) -10.0 W56-W73 N00-N06
      4 1116  779 (S19W10 L=360)   413    132        1 (1055/ 759-1193/ 810)  36.9 W06-W14 S21-S15
      5 1476  783 (S12W32 L= 22)   332    113        0 (1425/ 766-1525/ 796)   8.1 W28-W35 S12-S10
      6 1604  742 (S11W42 L= 32)   237     58        1 (1582/ 711-1622/ 768) -33.7 W40-W43 S12-S09
      7  342  648 (S40E45 L=304)   601    183        0 ( 300/ 604- 384/ 718)  38.7 E50-E40 S43-S35
      8 1109  658 (S26W13 L=  2)   488    190        0 (1028/ 613-1209/ 713)  45.0 W08-W17 S30-S21
      9 1613  644 (S16W46 L= 36)   434    112        1 (1569/ 613-1662/ 674)  39.8 W43-W49 S18-S13
     10 1817  560 (S15W71 L= 61)   747    157        1 (1798/ 537-1835/ 588) -21.8 W68-W75 S16-S13
     11 1593  533 (S23W50 L= 39)  1040    302        4 (1486/ 477-1672/ 583)   6.7 W41-W58 S25-S19
     12 1230  510 (S33W25 L= 15)   989    249        4 (1128/ 484-1312/ 546)  18.4 W18-W30 S35-S29
  

    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