Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  773 1478 (N35E05 L=176)   871    167        1 ( 685/1436- 857/1507) -10.3 E11-W00 N33-N37
      2 1691 1386 (N08W52 L=233)   441     91        1 (1658/1370-1728/1402) -18.4 W49-W55 N06-N09
      3  187 1028 (N20E58 L=123)   382    156        0 ( 150/ 987- 225/1053) -29.7 E61-E54 N17-N22
      4 1043  847 (S07E03 L=178)  3189    785       11 ( 855/ 594-1169/1038)  70.2 E13-W04 S21-N04
      5  166  880 (N12E64 L=117)  1493    228        5 ( 119/ 840- 210/ 918)  76.0 E70-E59 N09-N14
      6  845  700 (S12E18 L=163)   212     85        0 ( 812/ 664- 870/ 737)  31.0 E20-E15 S12-S09
      7 1238  627 (S25W04 L=185)   490    149        1 (1167/ 586-1315/ 670)   5.7 E01-W09 S25-S23
      8  612  332 (S29E48 L=133)   667    177        0 ( 575/ 277- 634/ 407)  63.4 E50-E44 S33-S23
  

    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