Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  734 1347 (N30E15 L=  5)  1336    296        3 ( 592/1307- 868/1400)   0.0 E25-E04 N28-N32
      2 1364 1237 (N14W25 L= 44)   536     88        1 (1319/1202-1412/1264) -26.6 W21-W27 N11-N15
      3 1716 1022 (S06W46 L= 66)   684    152        1 (1681/ 951-1747/1092)  53.1 W42-W49 S09-S01
      4  210 1041 (N16E60 L=320)   854    204        0 ( 157/1009- 265/1083) -22.6 E65-E53 N13-N18
      5  738  777 (S05E21 L=359)  1188    280        1 ( 706/ 659- 776/ 916)  74.1 E25-E18 S11-N03
      6  267  807 (N01E57 L=322)   498     70        1 ( 256/ 771- 285/ 853)  63.4 E59-E55 N00-N04
      7 1001  816 (S06E04 L= 15)   413    116        2 ( 937/ 788-1059/ 838) -29.7 E08-E01 S08-S04
      8 1782  755 (S24W50 L= 70)  1179    157        3 (1723/ 712-1827/ 798)   0.0 W44-W55 S26-S21
      9 1615  715 (S23W35 L= 55)   362     77        1 (1600/ 678-1627/ 754) -80.5 W33-W36 S25-S19
     10  930  516 (S24E14 L=  5)   395    104        0 ( 906/ 469- 955/ 560)  56.3 E16-E12 S27-S21
     11  693  287 (S39E41 L=339)  1003    153        1 ( 650/ 253- 723/ 329)  31.0 E46-E36 S41-S34
  

    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