Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1048 1717 (N36E16 L= 89)   617    105        1 ( 999/1684-1099/1740)  -6.3 E20-E11 N34-N37
      2  391 1535 (N12E54 L= 51)   884    164        3 ( 359/1456- 426/1599) -73.3 E57-E52 N06-N16
      3 1246 1507 (N26W03 L=108)   371    100        1 (1207/1465-1286/1554)  66.8 W01-W04 N22-N30
      4 1508 1507 (N33W21 L=126)  1603    426        1 (1300/1443-1659/1546)   7.9 W05-W34 N29-N37
      5 1092 1255 (N07E01 L=104)   289     92        0 (1048/1234-1140/1275)  36.9 E03-W01 N05-N09
      6  203 1212 (S10E62 L= 44)   508    125        0 ( 175/1161- 241/1255)  50.2 E64-E58 S12-S06
      7  328 1135 (S13E47 L= 58)   494    128        1 ( 267/1108- 386/1157)   0.0 E52-E41 S13-S11
      8 1832 1092 (N16W53 L=158)   770    107        0 (1786/1048-1869/1134) -18.4 W47-W57 N14-N17
      9  258 1087 (S16E53 L= 52)   327     65        0 ( 244/1046- 273/1121)  90.0 E54-E51 S18-S14
     10 1447 1036 (N03W24 L=129)   192     98        0 (1436/ 986-1462/1083) -68.2 W22-W25 N00-N05
     11 1017 1012 (S08E00 L=105)   428    209        0 ( 946/ 960-1100/1047)  40.6 E03-W04 S12-S05
     12 1863  875 (N05W61 L=166)  1146    245        0 (1835/ 759-1879/ 984) -67.4 W56-W64 S01-N11
     13  317  953 (S24E46 L= 59)   398     72        2 ( 276/ 927- 354/ 977)  14.0 E50-E42 S25-S21
     14  776  928 (S17E13 L= 93)   390    196        0 ( 697/ 889- 863/ 965)  33.7 E16-E07 S20-S14
     15 1145  836 (S15W11 L=116)   183     59        1 (1124/ 809-1167/ 862)  71.6 W10-W11 S16-S13
     16 1303  721 (S18W23 L=128)   245     43        1 (1280/ 697-1330/ 744)  56.3 W21-W24 S19-S16
     17 1153  681 (S24W16 L=121)   334    153        0 (1083/ 662-1215/ 716)   6.3 W10-W19 S24-S21
     18 1802  599 (S11W67 L=173)   618     70        3 (1776/ 566-1823/ 630) -36.9 W63-W71 S12-S09
  

    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