Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  380 1605 (N34E66 L=132)   579     51        0 ( 363/1574- 396/1634) -21.8 E72-E60 N32-N36
      2 1213 1467 (N32W10 L=208)  1626    455        4 (1079/1330-1382/1566) -30.6 W00-W22 N23-N38
      3  507 1387 (N21E39 L=158)   600    160        0 ( 478/1326- 556/1458) -38.7 E44-E34 N17-N25
      4 1694 1317 (N24W49 L=246)   422     83        1 (1677/1271-1710/1353) -90.0 W47-W50 N21-N26
      5 1692 1156 (N14W46 L=244)   439     98        1 (1659/1129-1724/1182) -31.0 W43-W49 N12-N15
      6 1861 1122 (N12W65 L=263)   961    119        1 (1837/1062-1887/1178)  60.3 W61-W70 N08-N15
      7  874 1094 (N06E10 L=188)   471    192        0 ( 828/1053- 925/1148) -51.3 E12-E06 N03-N08
      8 1749  857 (S04W53 L=251)   482     22        0 (1723/ 837-1772/ 880) -11.3 W50-W55 S05-S03
      9 1112  727 (S16W07 L=205)   464    131        1 (1067/ 690-1140/ 784)  81.9 W04-W09 S18-S11
     10 1826  653 (S17W70 L=268)   611     34        1 (1799/ 637-1845/ 674)  -5.7 W65-W75 S18-S16
     11 1298  638 (S20W21 L=219)   267     64        0 (1272/ 608-1331/ 670)  53.1 W19-W22 S22-S18
     12 1680  607 (S21W53 L=251)   424     38        1 (1664/ 579-1700/ 631) -45.0 W50-W55 S22-S19
     13  357  568 (S33E52 L=146)   702     76        1 ( 313/ 554- 400/ 590)   0.0 E57-E47 S33-S30
     14 1347  570 (S25W26 L=223)   250     83        0 (1313/ 555-1387/ 588)   9.5 W22-W28 S25-S23
     15 1060  478 (S33W07 L=204)   276     61        0 (1029/ 462-1094/ 496)  36.9 W04-W08 S34-S31
     16 1192  412 (S37W18 L=216)   381     95        0 (1149/ 389-1241/ 429)  23.2 W14-W21 S39-S35
  

    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