Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1297 1788 (N51E09 L=238)   656     83        0 (1252/1762-1340/1810)  14.0 E13-E04 N49-N53
      2  401 1627 (N13E63 L=185)   875     89        2 ( 366/1596- 433/1663)  59.0 E66-E60 N10-N15
      3 1563 1621 (N48W20 L=267)   396     65        0 (1533/1595-1591/1648) -15.9 W16-W23 N47-N50
      4  270 1434 (S01E64 L=183)   594     98        0 ( 241/1390- 291/1462) -18.4 E68-E59 S03-N00
      5 1667 1439 (N38W31 L=279)   414     40        0 (1644/1418-1690/1461) -14.0 W29-W33 N36-N39
      6  225 1331 (S08E64 L=183)   718     88        1 ( 203/1288- 246/1373)  90.0 E67-E61 S10-S05
      7  199 1194 (S16E63 L=185)  2268    322        4 (  88/1128- 352/1254)  15.5 E84-E47 S20-S10
      8  643 1176 (S08E26 L=222)   583    193        1 ( 610/1096- 707/1238)  78.7 E27-E23 S12-S02
      9 1323 1188 (N11W12 L=259)   490    109        1 (1274/1144-1366/1234)  -8.1 W08-W15 N08-N12
     10  976 1084 (S04E04 L=243)   295     70        0 ( 953/1050- 998/1117) -76.0 E05-E03 S06-S02
     11  451  959 (S24E34 L=214)   765    220        1 ( 407/ 871- 534/1026)  59.7 E36-E28 S30-S18
     12 1270  789 (S12W20 L=268)   245    145        0 (1211/ 766-1337/ 812)  39.8 W17-W23 S14-S08
     13 1383  666 (S15W31 L=279)   324     81        0 (1344/ 640-1404/ 701)  90.0 W29-W32 S17-S12
     14  528  543 (S49E22 L=226)  1238    252        2 ( 407/ 512- 641/ 602)   4.6 E35-E10 S51-S44
     15 1651  561 (S12W55 L=303)   331    105        1 (1629/ 512-1679/ 597) -78.7 W52-W57 S14-S09
     16 1210  525 (S28W27 L=275)   470    127        1 (1155/ 506-1274/ 553)  20.6 W23-W31 S29-S25
     17 1587  421 (S21W60 L=308)   842    102        3 (1548/ 388-1635/ 453)  51.3 W57-W63 S23-S18
     18 1455  413 (S26W50 L=297)   436     48        1 (1431/ 382-1478/ 438) -31.0 W47-W52 S27-S24
     19 1105  261 (S48W40 L=287)  1828    370        0 ( 948/ 184-1264/ 342)   6.3 W22-W59 S49-S44
  

    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