Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  947 1675 (N50W04 L=144)   281     54        1 ( 918/1652- 970/1690) -26.6 W00-W05 N48-N51
      2 1074 1644 (N46W14 L=155)   302     78        2 (1045/1618-1112/1662) -38.7 W11-W16 N43-N47
      3  839 1615 (N47E07 L=134)   424     98        1 ( 809/1567- 858/1644) -56.3 E09-E05 N42-N49
      4  809 1554 (N42E10 L=131)   286     63        1 ( 795/1519- 829/1584)  78.7 E11-E08 N39-N44
      5 1461 1442 (N26W36 L=177)   436    130        2 (1426/1398-1487/1486)  39.8 W32-W38 N23-N28
      6 1044  998 (N04W01 L=142)   411    151        0 ( 971/ 975-1117/1040) -20.6 E03-W05 N02-N07
      7  282  937 (N06E53 L= 87)   480    165        2 ( 241/ 914- 342/ 957)   5.7 E57-E47 N04-N07
      8 1818  913 (S12W56 L=197)   465    111        1 (1794/ 883-1849/ 944)  26.6 W53-W59 S14-S10
      9  960  497 (S27E11 L=130)   352    113        0 ( 931/ 459- 990/ 541)  45.0 E13-E08 S29-S24
     10  775  402 (S32E26 L=114)   365    135        2 ( 720/ 380- 826/ 432)  12.5 E30-E21 S33-S30
  

    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