Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  428 1674 (N17E65 L= 30)  1837    143        3 ( 378/1639- 473/1732) -13.4 E79-E58 N14-N21
      2 1664 1332 (N31W32 L=127)   278     40        1 (1640/1313-1684/1354)  45.0 W30-W34 N28-N32
      3 1147 1240 (N08W01 L= 96)   454    124        1 (1119/1190-1175/1298)  90.0 W00-W01 N04-N12
      4 1314 1255 (N14W10 L=105)   201     77        1 (1301/1215-1338/1295) -53.1 W08-W12 N12-N16
      5 1794 1236 (N30W46 L=141)   502    114        1 (1764/1190-1841/1258)  41.6 W43-W52 N25-N33
      6  688 1176 (S07E23 L= 72)   205     47        1 ( 674/1142- 702/1203) -76.0 E24-E22 S09-S05
      7  165 1081 (S24E66 L= 29)  1270    172        1 ( 114/1013- 232/1133)  27.9 E75-E57 S28-S19
  

    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