Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  727 1660 (N41E31 L=290)  4058    545        8 ( 487/1528- 931/1770) -11.0 E73-E11 N32-N47
      2 1142 1540 (N36W04 L=325)  1350    293        3 (1010/1496-1283/1564)   8.5 E05-W15 N32-N38
      3  426 1462 (N24E49 L=272)   373     63        0 ( 410/1441- 455/1496) -26.6 E51-E45 N22-N26
      4  238 1414 (N19E68 L=252)  1683    122        2 ( 198/1366- 268/1469) -50.2 E74-E64 N16-N22
      5 1428 1383 (N27W25 L=346)   600    203        1 (1365/1303-1490/1457) -48.4 W21-W29 N22-N31
      6 1221 1348 (N23W10 L=331)   482    142        0 (1194/1279-1237/1418)  83.7 W08-W11 N18-N27
      7 1121 1334 (N22W04 L=325)   323    132        2 (1095/1285-1154/1407) -63.4 W02-W06 N18-N26
      8  342 1139 (N03E47 L=273)  2053    336        2 ( 309/ 987- 375/1296) -86.8 E49-E44 S05-N13
      9 1833 1237 (N20W62 L= 22)  1139    145        3 (1776/1204-1878/1276)   0.0 W54-W68 N17-N22
     10 1623 1245 (N19W40 L=  1)   312     98        0 (1574/1235-1669/1255)   7.1 W35-W43 N18-N20
     11 1898 1158 (N15W70 L= 30)  1292     61        1 (1859/1125-1930/1194)  29.7 W63-W77 N13-N17
     12  182  971 (S08E63 L=258)   477    112        1 ( 149/ 929- 211/1018)  33.7 E68-E59 S11-S05
     13 1907  959 (N03W70 L= 31)   610     70        1 (1882/ 927-1936/ 981)  21.8 W66-W76 N00-N04
     14  564  769 (S17E28 L=292)   250     95        0 ( 518/ 752- 601/ 793)  -9.5 E31-E25 S18-S16
     15 1333  656 (S19W23 L=344)  1173    238        1 (1229/ 618-1454/ 704)  19.7 W16-W30 S21-S14
     16  364  679 (S25E46 L=274)   306     75        1 ( 340/ 654- 396/ 699)  26.6 E49-E42 S27-S23
     17 1029  520 (S30W04 L=325)   722    140        1 ( 963/ 498-1102/ 543)  12.5 W00-W09 S32-S28
     18  666  375 (S45E25 L=296)   302     29        1 ( 645/ 360- 687/ 396)   0.0 E27-E22 S46-S43
  

    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