Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  763 1759 (N53E31 L=159)   449     80        2 ( 716/1739- 799/1777)   7.1 E35-E27 N51-N54
      2 1587 1679 (N48W58 L=248)   610     53        0 (1552/1658-1617/1698)   7.6 W50-W65 N45-N49
      3 1188 1547 (N37W11 L=201)   417     94        0 (1144/1525-1231/1568)   8.1 W07-W14 N35-N38
      4  305 1524 (N31E67 L=123)   698     88        1 ( 273/1507- 347/1546)  -6.7 E76-E59 N30-N32
      5 1246 1491 (N33W15 L=205)   272     43        1 (1224/1470-1266/1522)  76.0 W13-W16 N31-N35
      6  460 1475 (N29E45 L=145)  1128    214        3 ( 390/1440- 562/1508)  -6.7 E52-E35 N26-N31
      7 1265 1325 (N22W15 L=205)   404    106        2 (1245/1269-1287/1373)  66.8 W13-W16 N18-N25
      8  285 1168 (N09E53 L=137)   810    176        0 ( 261/1096- 316/1256) -51.3 E57-E49 N04-N14
      9 1676 1048 (N05W44 L=234)  1512    335        0 (1590/ 919-1720/1181) -59.0 W37-W47 S02-N13
     10  591 1066 (N04E28 L=163)   230     48        0 ( 574/1037- 613/1090)  71.6 E28-E26 N02-N05
     11 1843  808 (S10W63 L=253)   873    109        2 (1816/ 755-1867/ 865) -74.1 W60-W66 S13-S06
     12 1580  752 (S13W38 L=228)   313     57        1 (1561/ 728-1597/ 780) -56.3 W36-W39 S14-S11
     13  265  698 (S21E58 L=132)   769    140        0 ( 212/ 660- 308/ 738)  20.6 E67-E51 S23-S17
     14  588  635 (S23E29 L=161)   243    107        0 ( 575/ 584- 608/ 683)  74.1 E30-E27 S26-S19
     15 1763  588 (S24W62 L=252)  2850    203        5 (1683/ 523-1817/ 663)  48.4 W54-W69 S28-S19
     16  870  578 (S26E09 L=181)   831    182        0 ( 802/ 545- 953/ 635) -21.8 E13-E03 S27-S21
     17  659  584 (S26E24 L=166)   403    123        0 ( 601/ 562- 710/ 620)  20.6 E28-E20 S27-S23
     18 1284  218 (S55W33 L=223)  2626    218        8 (1186/ 183-1401/ 252)  11.3 W21-W49 S59-S51
  

    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