Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  433 1677 (N36E59 L=166)   715    144        1 ( 394/1647- 478/1693)  33.7 E62-E55 N33-N37
      2  269 1318 (N11E59 L=166)   442    139        0 ( 238/1292- 303/1355) -12.5 E64-E55 N09-N12
      3  656 1296 (N14E27 L=199)   360    192        1 ( 589/1217- 706/1345) -32.5 E32-E21 N09-N16
      4  470 1156 (N03E38 L=188)   490     78        0 ( 437/1113- 502/1201)  63.4 E39-E35 N00-N06
      5  895  939 (S05E07 L=219)   670    167        0 ( 866/ 845- 923/1010) -60.9 E09-E04 S10-N00
      6  578  904 (S11E27 L=198)   276    100        1 ( 530/ 884- 621/ 923)   9.5 E30-E24 S12-S09
      7  990  709 (S19W01 L=227)   267     81        1 ( 972/ 668-1022/ 739)  68.2 W00-W03 S21-S16
  

    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