Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1427 1450 (N19W37 L=249)   563    133        1 (1359/1430-1492/1466) -32.0 W32-W40 N16-N21
      2 1278 1312 (N16W22 L=234)   203     43        0 (1261/1286-1295/1343)  63.4 W21-W23 N14-N18
      3 1437 1264 (N09W30 L=242)   222    107        0 (1400/1221-1456/1310)  26.6 W26-W33 N07-N11
      4  648 1190 (N25E19 L=193)   300     81        0 ( 606/1172- 700/1208) -26.6 E21-E15 N23-N27
      5  222 1099 (N28E57 L=155)  1285    222        1 ( 135/1037- 294/1170) -24.4 E72-E49 N23-N33
      6  586 1003 (N16E27 L=185)   226     94        0 ( 554/ 967- 617/1032) -59.0 E28-E25 N13-N18
      7  837  991 (N10E11 L=201)   231     68        0 ( 809/ 964- 862/1021) -63.4 E12-E09 N07-N11
      8 1217  685 (S17W02 L=214)   283     88        1 (1187/ 648-1244/ 717) -90.0 W00-W03 S19-S14
      9 1832  686 (S38W53 L=265)   887    102        1 (1785/ 657-1883/ 715)  -3.8 W46-W63 S39-S36
     10 1068  644 (S16E07 L=205)   668    218        1 (1002/ 589-1151/ 703)  15.3 E12-E01 S17-S13
     11 1544  664 (S28W23 L=235)   281     56        0 (1529/ 636-1564/ 700) -90.0 W21-W24 S30-S25
     12  966  661 (S12E12 L=199)   245     87        0 ( 931/ 638-1010/ 680)   0.0 E14-E09 S13-S10
     13 1165  590 (S21E03 L=209)   388    117        0 (1118/ 551-1194/ 640)  33.7 E06-W00 S22-S18
     14  240  572 (S04E72 L=140)   648     53        0 ( 215/ 553- 267/ 595) -26.6 E75-E66 S04-S01
     15  978  521 (S21E16 L=196)   438     87        2 ( 937/ 500-1028/ 548) -38.7 E18-E13 S23-S18
     16  610  419 (S19E44 L=168)   376     53        1 ( 586/ 398- 634/ 441)  21.8 E45-E40 S20-S17
     17 1487  389 (S47W15 L=227)  1204    170        1 (1412/ 354-1562/ 430)  -3.8 W07-W22 S48-S44
     18  617  260 (S30E56 L=156)  1445    149        3 ( 591/ 195- 640/ 344)  39.1 E66-E49 S36-S23
  

    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