Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1070 1748 (N38E21 L=101)  1237    189        0 (1005/1684-1155/1784)  71.6 E23-E16 N31-N43
      2  510 1627 (N15E51 L= 70)   448     92        2 ( 484/1589- 533/1671) -63.4 E53-E48 N12-N18
      3 1525 1484 (N34W19 L=140)   438     78        0 (1499/1450-1548/1522) -26.6 W16-W21 N31-N36
      4 1709 1284 (N28W37 L=158)   576     88        0 (1674/1239-1739/1329) -33.7 W33-W39 N25-N29
      5 1018 1306 (N08E08 L=114)   484    130        1 ( 966/1271-1091/1325)  26.6 E10-E03 N05-N11
      6  517 1204 (S09E34 L= 87)   402    124        0 ( 478/1154- 542/1261)  81.9 E35-E33 S12-S05
      7  933 1068 (S07E06 L=115)   645    181        1 ( 877/1018- 974/1139) -26.6 E09-E03 S09-S03
      8  670 1094 (S12E22 L=100)   469    110        1 ( 638/1052- 715/1134)  71.6 E23-E19 S14-S08
      9  852 1038 (S11E10 L=112)   236     41        1 ( 827/1015- 873/1061)  71.6 E11-E08 S12-S09
     10 1019  864 (S16W04 L=125)   562    219        1 ( 928/ 848-1138/ 884)  31.0 E01-W10 S18-S12
     11 1809  848 (N07W55 L=177)   471     86        1 (1788/ 815-1832/ 881) -18.4 W52-W58 N05-N08
     12 1184  779 (S16W16 L=137)   266     88        1 (1134/ 765-1229/ 794)  18.4 W12-W18 S17-S14
     13 1158  684 (S22W17 L=139)   237    111        0 (1110/ 668-1198/ 716)  45.0 W14-W18 S23-S18
     14 1675  589 (S11W55 L=177)   517    103        0 (1645/ 540-1705/ 628) -78.7 W53-W57 S13-S08
     15 1151  543 (S30W23 L=144)   215     32        1 (1123/ 524-1170/ 561)  14.0 W20-W24 S30-S28
     16 1658  463 (S18W63 L=184)  1539    258        1 (1543/ 427-1737/ 493)  24.0 W53-W72 S22-S14
     17  985  367 (S46W22 L=143)   285    104        0 ( 961/ 326-1027/ 409) -21.8 W17-W27 S47-S43
     18 1439  363 (S30W53 L=174)   620     39        2 (1409/ 336-1467/ 394) -31.0 W49-W55 S31-S28
  

    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