Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1340 1285 (N14W25 L=  7)  1997    514        5 (1232/1150-1520/1465) -71.6 W21-W33 N03-N27
      2  684 1287 (N30E16 L=326)  1267    387        2 ( 550/1214- 849/1371) -38.3 E25-E06 N22-N37
      3 1804 1225 (S02W58 L= 40)   468     31        1 (1780/1204-1828/1245)  -9.5 W55-W61 S02-N00
      4  949 1126 (N15E02 L=339)   443    275        0 ( 851/1049-1051/1202)  15.9 E09-W05 N12-N17
      5 1884 1011 (S17W63 L= 44)   480    104        2 (1875/ 958-1896/1065)  51.3 W60-W65 S19-S13
      6  420  987 (N17E40 L=302)   648    116        2 ( 362/ 960- 466/1020)   5.7 E45-E35 N15-N18
      7  260  970 (N18E55 L=286)  1038    160        0 ( 187/ 946- 326/1009) -14.9 E64-E48 N15-N20
      8  627  910 (N09E26 L=315)   883    207        3 ( 553/ 863- 715/ 960)   0.0 E31-E19 N05-N10
      9 1579  859 (S16W31 L= 12)   231     53        0 (1565/ 827-1595/ 885)  45.0 W28-W31 S17-S14
     10 1158  566 (S23E03 L=339)   235     78        1 (1118/ 551-1198/ 581) -45.0 E04-W00 S25-S21
     11  371  517 (S11E59 L=282)   439     37        0 ( 357/ 487- 389/ 540)  53.1 E61-E56 S13-S09
     12 1196  455 (S32E03 L=338)   386     72        2 (1175/ 422-1216/ 497)  59.0 E05-W00 S34-S28
  

    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