Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  354 1552 (N16E58 L= 32)  1023    241        2 ( 309/1487- 385/1594) -69.4 E64-E52 N11-N19
      2 1708 1528 (N40W43 L=133)   399     97        0 (1688/1498-1734/1549)  45.0 W40-W46 N37-N42
      3 1525 1428 (N27W25 L=115)   824    199        2 (1429/1366-1600/1486) -12.1 W17-W31 N24-N29
      4 1464 1180 (N10W24 L=114)   499    114        1 (1434/1109-1483/1232) -66.8 W21-W25 N06-N13
      5 1687 1076 (N09W41 L=131)   459    151        2 (1641/1016-1731/1120) -51.3 W36-W44 N05-N11
      6  262 1059 (S15E52 L= 38)   387     96        1 ( 232/1033- 300/1084) -15.9 E55-E48 S16-S13
      7 1852  893 (N03W60 L=150)   699    180        0 (1810/ 846-1905/ 926)  26.6 W56-W66 N00-N06
      8 1059  750 (S22W07 L= 97)   381    131        0 ( 999/ 723-1119/ 796)  45.0 W03-W10 S24-S18
      9 1719  562 (S18W59 L=149)   896    143        2 (1665/ 509-1779/ 625)  60.9 W55-W62 S22-S13
     10  916  477 (S42W04 L= 94)   233     45        0 ( 889/ 459- 941/ 499)  45.0 W02-W05 S43-S40
     11 1064  370 (S46W19 L=109)   600    104        1 (1019/ 343-1120/ 399)   0.0 W15-W24 S47-S44
     12  862  204 (S67W11 L=101)   589     91        1 ( 839/ 174- 884/ 229) -26.6 W06-W17 S69-S64
  

    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