Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  888 1593 (N45E01 L=236)   308     65        0 ( 857/1575- 925/1612) -36.9 E03-W01 N42-N46
      2  445 1536 (N46E45 L=192)   701    122        1 ( 399/1493- 492/1574) -25.0 E52-E37 N42-N49
      3 1485 1463 (N25W40 L=277)   382    125        1 (1427/1445-1542/1482) -20.6 W35-W44 N23-N26
      4 1200 1437 (N29W18 L=255)   322    162        1 (1134/1423-1286/1451) -15.3 W13-W24 N26-N30
      5  735 1427 (N35E15 L=222)  1429    393        2 ( 548/1396- 911/1447) -12.1 E30-E01 N30-N37
      6 1445 1265 (N14W30 L=267)   952    283        3 (1385/1136-1474/1355) -85.9 W26-W33 N06-N20
      7 1746 1298 (N10W56 L=293)   756    167        0 (1663/1255-1802/1324) -30.3 W48-W60 N06-N13
      8  180  991 (N12E65 L=172)  1149    142        3 ( 124/ 953- 250/1015)   9.0 E75-E56 N09-N13
      9  868  844 (S03E12 L=225)   194    115        0 ( 826/ 809- 915/ 862)  18.4 E14-E08 S04-S02
     10  326  807 (N00E51 L=186)   448     91        1 ( 282/ 790- 368/ 822)   0.0 E55-E47 N00-N00
     11  797  693 (S11E19 L=218)  1049    336        4 ( 706/ 601- 873/ 805)  37.6 E25-E12 S16-S05
     12 1839  743 (S26W61 L=297)   554     77        2 (1821/ 703-1853/ 778) -51.3 W57-W64 S28-S23
     13  283  645 (S10E61 L=176)  4147    292       14 ( 200/ 528- 380/ 746) -54.5 E74-E54 S18-S03
     14  890  538 (S23E16 L=221)  1192    261        4 ( 787/ 497-1007/ 589)   0.0 E23-E07 S24-S18
     15  435  401 (S28E59 L=178)  1804    195        4 ( 405/ 308- 460/ 484)  34.3 E72-E50 S36-S21
     16 1102  169 (S59E13 L=224)   735     94        0 (1060/ 144-1138/ 195) -12.5 E18-E08 S63-S55
  

    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