Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1087 1808 (N57E10 L= 95)  1260    328        0 ( 933/1775-1231/1848)  -1.7 E29-W05 N54-N58
      2  341 1565 (N25E66 L= 40)  2386    365        3 ( 255/1434- 459/1682)  74.5 E72-E61 N16-N34
      3 1056 1618 (N39E07 L= 99)  1152    333        1 ( 899/1579-1196/1679)  -2.4 E20-W04 N36-N41
      4  425 1561 (N26E55 L= 51)   481    122        1 ( 404/1504- 451/1610) -74.1 E57-E52 N22-N29
      5 1631 1369 (N30W40 L=146)  2831    449        7 (1426/1301-1813/1441)  -6.2 W22-W59 N26-N32
      6 1803 1100 (N14W55 L=161)   513    128        0 (1785/1041-1838/1131)  63.4 W53-W59 N10-N16
      7  356 1073 (S05E45 L= 61)   480    110        0 ( 328/1018- 390/1116)  60.3 E47-E42 S08-S01
      8  614 1035 (S04E25 L= 80)   350    101        0 ( 587/ 994- 633/1089)  80.5 E26-E24 S06-N00
      9  889 1029 (S01E08 L= 98)   424     84        1 ( 854/ 983- 926/1068)  63.4 E09-E06 S03-N01
     10  323  906 (S15E47 L= 59)   332     86        0 ( 307/ 862- 335/ 951) -68.2 E48-E45 S17-S12
     11 1176  712 (S16W13 L=119)   745    291        0 (1087/ 664-1234/ 773)  26.6 W08-W16 S19-S12
  

    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