Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  928 1872 (N72W05 L=153)   649     69        1 ( 911/1839- 953/1899)  56.3 W00-W09 N67-N76
      2  788 1844 (N70E23 L=126)   471     46        0 ( 771/1814- 811/1868) -74.1 E26-E18 N66-N73
      3 1570 1558 (N33W52 L=200)   407     84        0 (1530/1543-1611/1572)  -7.1 W47-W55 N31-N33
      4 1105 1528 (N37W12 L=160)   384    106        0 (1058/1498-1154/1551)  14.0 W07-W15 N35-N38
      5 1177 1507 (N34W17 L=165)   337     74        2 (1157/1468-1196/1540)  68.2 W15-W18 N31-N36
      6 1694 1472 (N25W59 L=208)   872    149        2 (1650/1421-1744/1540)  82.9 W51-W64 N21-N29
      7 1539 1474 (N28W44 L=193)   729    169        0 (1500/1388-1570/1536)  79.7 W41-W45 N21-N32
      8 1774 1469 (N23W70 L=219)  1053    126        0 (1722/1435-1824/1495) -22.6 W64-W77 N20-N25
      9  786 1350 (N28E13 L=135)   497    174        1 ( 758/1278- 829/1420) -68.2 E15-E10 N22-N32
     10  378 1235 (N23E45 L=103)  3416    614        6 ( 271/1028- 556/1341) -87.1 E56-E30 N09-N29
     11 1241 1263 (N18W16 L=165)   395    130        1 (1200/1230-1292/1303) -33.7 W13-W19 N15-N20
     12  474 1001 (N08E36 L=112)   264     53        0 ( 453/ 973- 492/1033)  76.0 E38-E34 N05-N09
     13 1744  859 (S13W50 L=198)   512     60        0 (1724/ 832-1769/ 906)  90.0 W47-W52 S14-S10
     14  658  796 (S06E25 L=123)   210     92        0 ( 616/ 769- 694/ 812) -26.6 E28-E23 S07-S04
     15  235  725 (S09E64 L= 85)   793     76        2 ( 220/ 680- 254/ 767)  80.5 E67-E60 S12-S06
     16  833  398 (S35E21 L=127)   455    118        0 ( 776/ 386- 891/ 413) -12.5 E25-E16 S35-S33
     17 1230  248 (S54W10 L=159)   456     49        0 (1208/ 232-1258/ 269) -18.4 W07-W13 S56-S51
     18  942  185 (S57E23 L=125)   467     71        0 ( 907/ 171- 974/ 199)  -8.1 E26-E18 S58-S55
  

    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