Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  562 1558 (N51E25 L=255)   297     68        2 ( 540/1534- 578/1580) -53.1 E27-E23 N48-N52
      2 1576 1503 (N21W50 L=330)   893    151        3 (1516/1472-1632/1526) -32.0 W45-W53 N18-N23
      3 1124 1474 (N31W18 L=298)   616    140        4 (1077/1415-1159/1514) -81.9 W15-W19 N27-N34
      4 1556 1396 (N16W43 L=324)  1130    231        4 (1474/1309-1643/1461) -59.7 W39-W47 N08-N20
      5  594 1313 (N33E23 L=258)  1499    365        4 ( 523/1173- 663/1453) -68.2 E28-E19 N23-N43
      6 1779 1389 (N09W63 L=344)  1261    169        3 (1735/1318-1826/1449) -68.2 W61-W66 N03-N13
      7 1746 1296 (N05W55 L=336)  1245    281        5 (1699/1163-1791/1401)  59.7 W47-W60 S01-N11
      8 1272 1322 (N19W22 L=302)   500    110        1 (1226/1258-1311/1362) -81.9 W19-W23 N14-N21
      9 1751 1052 (S09W48 L=329)  1030    245        0 (1694/ 997-1827/1129) -55.0 W44-W55 S13-S03
     10  907 1093 (N13E05 L=275)   201     60        0 ( 875/1072- 945/1109)  14.0 E07-E03 N12-N14
     11  189 1010 (N20E61 L=219)  3277    490        6 ( 103/ 901- 268/1086) -36.9 E71-E53 N12-N24
     12  149  903 (N13E69 L=211)   576     82        2 ( 128/ 879- 170/ 923) -11.3 E72-E66 N11-N14
     13 1786  875 (S20W48 L=329)   396     90        1 (1751/ 840-1811/ 912)  26.6 W44-W51 S21-S18
     14 1376  740 (S17W15 L=296)  1517    308        2 (1289/ 613-1432/ 903)  47.0 W07-W21 S23-S08
     15  445  828 (N06E41 L=240)   353     64        2 ( 421/ 801- 465/ 862)  71.6 E42-E39 N04-N07
     16  457  641 (S06E45 L=236)   849    217        3 ( 382/ 614- 532/ 670)  -4.8 E51-E39 S07-S04
     17  305  488 (S14E68 L=212)   748    129        2 ( 280/ 450- 329/ 512) -59.0 E71-E64 S17-S12
     18  800  454 (S23E28 L=253)   534     84        2 ( 763/ 424- 833/ 484)   9.5 E30-E24 S25-S21
     19 1220  383 (S38E02 L=278)   388     46        0 (1191/ 363-1254/ 401) -26.6 E04-W00 S40-S36
  

    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