Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  639 1734 (N26E49 L=285)   780    116        1 ( 608/1697- 675/1775) -16.7 E53-E43 N23-N29
      2  450 1709 (N21E65 L=269)  1676    145        4 ( 378/1665- 506/1749)  33.7 E70-E59 N17-N24
      3  625 1619 (N18E42 L=292)   388     68        1 ( 610/1584- 649/1649) -78.7 E44-E39 N15-N20
      4 1018 1452 (N17E11 L=323)   327     61        0 (1001/1420-1035/1485) -45.0 E12-E09 N15-N18
      5 1708 1220 (N23W38 L= 12)  1368    316        2 (1656/1069-1767/1366) -53.1 W32-W46 N15-N31
      6  897  970 (S13E06 L=328)   236     51        2 ( 873/ 944- 921/ 991)  18.4 E07-E04 S14-S11
      7  183  714 (S44E70 L=264)   875     98        0 ( 155/ 668- 212/ 754)  12.8 E84-E62 S46-S41
      8 1707  646 (S08W56 L= 30)   362     39        1 (1687/ 624-1726/ 669) -71.6 W53-W57 S09-S06
      9  963  606 (S32W08 L=342)   304     79        1 ( 934/ 587-1008/ 635)   0.0 W05-W11 S33-S31
     10  564  583 (S45E21 L=313)   254     73        0 ( 522/ 573- 600/ 594)   8.1 E24-E17 S46-S43
     11  543  393 (S61E23 L=311)   449     63        0 ( 508/ 374- 569/ 424)   6.3 E28-E19 S62-S57
     12  706  328 (S61E01 L=333)   532     67        1 ( 677/ 298- 739/ 357) -14.0 E04-W04 S62-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