Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  612 1425 (N37E25 L=301)   247     78        1 ( 574/1397- 647/1449) -23.2 E28-E21 N34-N38
      2 1533 1280 (N13W37 L=  3)   835    189        1 (1499/1204-1577/1363)  84.3 W34-W39 N08-N18
      3  807 1218 (N21E11 L=315)   622    268        0 ( 772/1120- 828/1349) -86.2 E13-E08 N14-N29
      4 1124 1224 (N17W09 L=335)   410    188        2 (1051/1187-1209/1260) -36.9 W04-W13 N13-N19
      5  399 1162 (N22E42 L=284)   469    122        2 ( 355/1117- 449/1204) -36.9 E46-E37 N18-N24
      6 1723  984 (S08W47 L= 12)   723    186        3 (1681/ 922-1761/1048)  54.5 W42-W50 S11-S04
      7 1764  938 (S12W50 L= 16)   355     90        0 (1737/ 903-1782/ 976)  45.0 W47-W52 S13-S09
      8 1686  717 (S24W43 L=  8)  1116    235        2 (1624/ 623-1743/ 810)  45.0 W37-W47 S29-S18
      9 1550  747 (S19W31 L=356)   273     67        0 (1539/ 705-1562/ 782)  90.0 W29-W31 S21-S16
     10  484  413 (S27E53 L=273)  2419    297        4 ( 453/ 290- 523/ 543)  32.5 E76-E42 S38-S17
     11  332  443 (S24E71 L=254)  1115     83        2 ( 307/ 412- 366/ 476) -78.7 E79-E67 S26-S21
     12  953  322 (S40E18 L=307)   305     55        0 ( 925/ 306- 984/ 343)  21.8 E20-E15 S41-S38
  

    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