Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  713 1701 (N48E32 L=342)   535    113        1 ( 689/1661- 734/1742) -54.5 E34-E27 N44-N51
      2  670 1616 (N41E31 L=342)   347     38        1 ( 649/1597- 690/1639)  71.6 E33-E29 N39-N42
      3 1049 1478 (N32W00 L= 14)   516    124        1 ( 990/1433-1095/1508) -33.7 E03-W03 N29-N34
      4 1525 1339 (N23W34 L= 48)  2147    562        6 (1293/1238-1709/1446) -23.4 W17-W48 N17-N30
      5 1146 1280 (N19W07 L= 21)  1009    301        3 (1091/1137-1201/1403) -70.6 W03-W10 N10-N27
      6  609 1342 (N21E29 L=344)   378     61        1 ( 591/1309- 630/1375) -63.4 E30-E27 N19-N23
      7 1296 1275 (N19W17 L= 31)   691    175        1 (1239/1194-1371/1353) -48.4 W13-W21 N14-N24
      8  955 1066 (N05E04 L=  9)  1323    410        2 ( 867/ 967-1103/1231) -45.0 E10-W04 N00-N15
      9  805  921 (S04E13 L=  0)   334     99        2 ( 756/ 905- 862/ 936)   0.0 E16-E09 S05-S02
     10 1131  245 (S53W14 L= 28)   342     70        0 (1083/ 232-1163/ 255) -14.0 W09-W17 S54-S51
  

    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