Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  992 1735 (N39E16 L=136)   440     87        1 ( 941/1718-1043/1751)  23.2 E19-E12 N37-N40
      2  695 1617 (N25E32 L=120)   449    129        0 ( 664/1578- 729/1664) -29.7 E35-E28 N23-N28
      3  449 1595 (N21E50 L=102)  1955    270        2 ( 373/1548- 531/1628)  -7.1 E57-E41 N17-N23
      4  311 1489 (N13E59 L= 94)  2836    466        9 ( 222/1434- 391/1550)  20.6 E68-E50 N09-N17
      5  502 1398 (N09E39 L=113)  1208    338        4 ( 482/1269- 535/1542) -64.8 E42-E34 N01-N18
      6  378 1414 (N09E49 L=103)   922    114        3 ( 328/1358- 433/1461)  45.0 E53-E45 N05-N12
      7  792 1361 (N10E19 L=133)   192     70        0 ( 754/1345- 836/1381)  21.8 E21-E16 N08-N11
      8 1586 1220 (N15W32 L=184)   259     68        0 (1563/1176-1602/1251) -53.1 W30-W33 N12-N16
      9 1097 1201 (N05W02 L=154)   393    141        0 (1040/1134-1128/1235) -38.7 E02-W04 N01-N06
     10  279 1021 (S15E50 L=102)  1154    342        0 ( 227/ 884- 376/1136)  58.0 E54-E43 S23-S07
     11 1755 1027 (N07W48 L=200)   364    135        1 (1709/ 982-1786/1052)   0.0 W43-W51 N05-N08
     12 1263  911 (S09W16 L=168)   578    199        2 (1192/ 828-1303/ 960) -45.0 W11-W18 S13-S06
     13 1385  761 (S16W26 L=179)   516    167        1 (1293/ 739-1461/ 793)  33.7 W21-W30 S18-S12
     14  631  749 (S28E21 L=131)   207     78        0 ( 594/ 729- 671/ 772)  -9.5 E24-E18 S29-S27
     15  512  716 (S32E30 L=122)   467    143        2 ( 465/ 679- 568/ 754)  40.6 E34-E26 S34-S28
     16  335  568 (S44E49 L=104)  1017    151        2 ( 298/ 521- 376/ 613)  35.0 E53-E43 S47-S40
     17  505  493 (S47E33 L=119)   896    184        2 ( 434/ 476- 577/ 524)  12.1 E40-E26 S48-S45
     18  862  205 (S66W07 L=160)   617     98        1 ( 839/ 178- 888/ 244) -63.4 W03-W11 S70-S62
  

    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