Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  970 1785 (N59W15 L= 48)  1378    278        3 ( 828/1740-1093/1825)   0.0 E02-W31 N57-N60
      2 1382 1514 (N30W34 L= 67)   590    247        0 (1285/1435-1462/1564) -54.0 W28-W37 N23-N34
      3  366 1353 (N35E46 L=347)  3324    560        5 ( 273/1136- 529/1520)  50.7 E53-E32 N21-N45
      4 1186 1459 (N30W18 L= 51)   317    131        0 (1129/1433-1246/1486) -39.8 W15-W21 N27-N32
      5 1840 1328 (N09W68 L=101)  1102    148        3 (1788/1286-1899/1370) -33.7 W62-W77 N04-N12
      6  183 1225 (N27E67 L=326)  2383    210        3 ( 134/1129- 239/1308) -21.0 E84-E57 N21-N31
      7  513 1223 (N25E32 L=  1)   412    205        1 ( 441/1168- 580/1297) -39.8 E38-E26 N20-N30
      8 1342 1216 (N13W22 L= 55)   288    116        0 (1304/1172-1388/1241)   9.5 W19-W25 N10-N13
      9  768 1177 (N19E14 L= 19)   261    114        0 ( 725/1146- 817/1222) -45.0 E16-E10 N16-N22
     10  431  997 (N12E39 L=354)   482    130        1 ( 411/ 949- 450/1071)  82.9 E41-E37 N08-N16
     11  380  934 (N08E44 L=349)   297     77        1 ( 362/ 899- 403/ 969) -76.0 E45-E42 N06-N10
     12 1561  822 (S14W31 L= 64)   910    201        0 (1476/ 770-1641/ 872)   8.7 W24-W37 S16-S12
     13  652  600 (S15E30 L=  3)  1507    489        1 ( 498/ 522- 769/ 719)  -3.0 E40-E21 S19-S08
     14  410  633 (S10E48 L=345)   891    128        4 ( 387/ 574- 446/ 682)  45.0 E51-E43 S14-S07
     15 1549  436 (S42W32 L= 65)   366     49        0 (1523/ 420-1576/ 453)  11.3 W29-W35 S43-S40
     16 1385  340 (S46W17 L= 50)   787    105        1 (1342/ 313-1438/ 373)   6.3 W13-W22 S48-S42
     17 1219  286 (S48W02 L= 35)   445     56        0 (1191/ 264-1259/ 309) -31.0 W00-W05 S49-S44
     18 1170  230 (S52E04 L= 29)  1159    115        4 (1120/ 200-1229/ 265) -12.5 E08-W01 S55-S48
     19 1185  164 (S62E06 L= 27)   627     54        1 (1157/ 147-1210/ 184) -56.3 E09-E02 S65-S58
  

    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