Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1741 1597 (N20W71 L=297)   775     67        0 (1706/1579-1773/1622)   3.2 W65-W83 N18-N20
      2 1777 1462 (N11W64 L=290)  1691    120        2 (1715/1413-1820/1499) -36.9 W58-W70 N08-N14
      3  692 1365 (N28E14 L=212)   451    161        0 ( 633/1315- 750/1426) -55.0 E17-E10 N23-N33
      4 1634 1294 (N05W44 L=269)   851    199        1 (1552/1245-1701/1336)   4.1 W36-W50 N02-N06
      5  411 1240 (N26E36 L=190)   336     96        0 ( 363/1214- 446/1274)  15.9 E40-E33 N24-N27
      6 1139 1237 (N12W11 L=237)   341     88        0 (1100/1214-1189/1266) -45.0 W08-W13 N09-N13
      7 1337 1171 (N04W21 L=247)   347    101        0 (1294/1139-1375/1208) -59.0 W19-W22 N01-N06
      8  403 1117 (N19E38 L=188)  1462    462        1 ( 241/1060- 608/1169) -19.0 E52-E23 N12-N24
      9  835  874 (S03E14 L=212)   345    158        0 ( 744/ 849- 895/ 891) -24.0 E19-E10 S05-S01
     10  803  740 (S10E18 L=207)  1032    308        3 ( 713/ 693- 938/ 816) -24.8 E22-E09 S13-S04
     11 1154  651 (S22W00 L=226)   275    102        0 (1127/ 609-1178/ 700)  51.3 E02-W02 S24-S19
     12  304  671 (S05E56 L=170)   451     79        1 ( 261/ 645- 332/ 698)   0.0 E60-E52 S06-S03
     13  486  659 (S09E41 L=185)   513     93        2 ( 447/ 626- 522/ 695) -18.4 E44-E38 S11-S07
     14 1453  543 (S36W19 L=245)  1624    504        2 (1252/ 402-1646/ 679)  16.4 W01-W35 S41-S30
     15  321  541 (S13E61 L=165)  6613    291       14 ( 230/ 427- 410/ 648) -75.1 E73-E51 S21-S06
     16  632  528 (S20E36 L=190)   402     69        1 ( 606/ 485- 655/ 564) -90.0 E36-E34 S22-S17
     17  861  321 (S38E30 L=196)   993    171        0 ( 793/ 263- 916/ 400)  23.6 E38-E22 S41-S33
     18 1024  256 (S47E22 L=204)   508     69        1 ( 988/ 234-1063/ 279)   0.0 E25-E17 S48-S45
  

    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