Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1498 1695 (N29W55 L= 76)  1973    241        3 (1411/1631-1596/1730)   2.5 W44-W68 N25-N32
      2 1654 1668 (N23W69 L= 90)   547     39        2 (1629/1651-1681/1684) -18.4 W65-W73 N22-N24
      3 1612 1539 (N17W53 L= 74)   664    163        1 (1541/1497-1679/1579) -48.8 W49-W57 N13-N21
      4  253 1225 (N32E51 L=330)  1258    262        1 ( 168/1123- 334/1293) -32.5 E65-E42 N23-N37
      5 1841 1187 (S08W59 L= 80)  2686    403        9 (1781/1106-1919/1288) -58.4 W52-W68 S13-N00
      6 1702 1160 (S05W45 L= 66)   684    179        4 (1664/1073-1753/1215)  31.0 W40-W50 S09-S01
      7 1433 1065 (S04W24 L= 45)   982    331        4 (1376/ 965-1493/1160)  74.7 W20-W29 S10-N01
      8 1142  936 (S04W05 L= 26)  1494    523        3 ( 972/ 754-1264/1070) -66.8 E03-W07 S17-N04
      9 1556 1009 (S10W31 L= 52)   282     76        2 (1522/ 979-1586/1037)  11.3 W28-W33 S12-S08
     10  582  971 (N10E27 L=354)   510    158        1 ( 503/ 950- 655/ 993) -29.1 E31-E22 N07-N12
     11 1611  747 (S27W31 L= 52)  1395    277        8 (1529/ 694-1716/ 830) -25.0 W24-W40 S30-S19
     12 1428  609 (S31W15 L= 36)  1553    310        3 (1369/ 503-1484/ 770) -81.0 W11-W18 S38-S19
     13  705  638 (S12E28 L=353)   919    207        4 ( 622/ 580- 795/ 698)  -4.8 E34-E22 S13-S07
     14  673  468 (S21E36 L=345)  2932    643        8 ( 457/ 322- 874/ 664)  14.3 E63-E16 S26-S13
     15 1372  487 (S37W09 L= 30)   638    125        1 (1331/ 446-1419/ 519)  14.0 W05-W13 S39-S34
     16 1323  352 (S46W02 L= 23)  2059    262        4 (1217/ 296-1432/ 427)   6.0 E07-W12 S50-S41
  

    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