Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  828 1860 (N67W22 L= 41)   728     77        2 ( 804/1819- 854/1891)  23.6 W15-W31 N62-N71
      2 1080 1759 (N47W31 L= 50)  1197    336        1 ( 991/1647-1184/1867)  15.9 W17-W52 N41-N51
      3  743 1813 (N66W07 L= 26)   592     64        2 ( 721/1774- 762/1845)  48.8 W03-W11 N61-N69
      4  679 1765 (N63E04 L= 15)   371     54        0 ( 655/1747- 700/1782)   9.5 E07-E01 N61-N65
      5  854 1744 (N54W12 L= 32)   521     94        1 ( 808/1711- 901/1773)   4.8 W06-W18 N52-N55
      6 1605 1683 (N25W67 L= 86)   656     66        2 (1577/1656-1636/1718)  38.7 W63-W73 N23-N27
      7 1293 1577 (N28W33 L= 53)   616    184        1 (1243/1538-1375/1638) -66.0 W32-W37 N24-N33
      8  737 1335 (N28E10 L=  9)   337    130        1 ( 684/1282- 782/1364) -50.2 E13-E07 N24-N31
      9 1087 1322 (N19W11 L= 30)   458    189        0 ( 997/1283-1171/1353) -37.9 W06-W15 N15-N23
     10 1417 1320 (N11W31 L= 50)   379     74        0 (1393/1282-1452/1346)  21.8 W28-W33 N08-N12
     11  213 1075 (N24E57 L=322)   458     84        1 ( 196/1032- 229/1114) -78.7 E59-E55 N21-N26
     12  187  949 (N17E62 L=317)   611     65        0 ( 162/ 922- 214/ 980) -26.6 E65-E58 N15-N19
     13  171  879 (N14E66 L=313)   554     81        3 ( 152/ 839- 189/ 915) -68.2 E68-E63 N10-N15
     14 1819  891 (S25W52 L= 71)   411     62        1 (1788/ 871-1852/ 915) -20.6 W48-W56 S26-S22
     15 1336  882 (S12W14 L= 34)   337    144        1 (1272/ 862-1389/ 909) -14.0 W10-W18 S13-S09
     16  628  574 (S13E36 L=344)   394     75        1 ( 613/ 526- 650/ 616)  59.0 E37-E33 S15-S10
     17 1117  358 (S39E14 L=  6)  1401    291        1 (1016/ 266-1218/ 449)  20.0 E25-E03 S44-S35
     18 1375  237 (S61W02 L= 22)   509     64        1 (1342/ 216-1406/ 254)  -8.1 E01-W06 S63-S59
  

    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