Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  985 1457 (N17E15 L=349)   410    124        1 ( 926/1420-1046/1485)   7.1 E18-E10 N16-N18
      2 1803 1381 (N41W46 L= 49)   696    117        1 (1770/1363-1841/1403)  20.6 W42-W51 N38-N43
      3 1642 1346 (N32W30 L= 33)   711    170        0 (1587/1265-1696/1392) -21.8 W25-W35 N28-N33
      4  980 1268 (N06E09 L=354)   673    173        3 ( 946/1183-1024/1350) -53.1 E12-E04 N02-N10
      5 1251 1149 (N07W09 L= 12)   346    142        1 (1215/1085-1303/1212) -26.6 W05-W13 N04-N09
      6  326  972 (S26E44 L=319)   499    116        1 ( 286/ 921- 357/1020) -26.6 E48-E40 S28-S24
      7 1143  767 (S16W14 L= 17)   764    190        2 (1074/ 716-1230/ 843)  -8.7 W07-W20 S18-S13
  

    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