Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1561 1508 (N18W49 L= 14)   690    148        2 (1487/1475-1613/1526) -41.2 W44-W53 N14-N21
      2  581 1170 (N26E24 L=301)   624    194        2 ( 517/1094- 652/1238) -61.4 E27-E20 N20-N31
      3 1822 1150 (S09W56 L= 21)   518    147        3 (1808/1084-1844/1213)  45.0 W53-W60 S13-S06
      4  162 1023 (N25E62 L=263)   454    137        0 ( 133/ 966- 194/1057) -50.2 E65-E59 N21-N27
      5 1209  820 (S09W05 L=330)   894    447        3 (1047/ 768-1390/ 890) -33.7 E03-W15 S15-S02
      6  395  859 (N12E44 L=281)  1039    326        2 ( 272/ 817- 516/ 886) -20.2 E54-E35 N07-N14
      7 1304  639 (S23W06 L=331)   236     57        1 (1283/ 615-1337/ 660)  14.0 W04-W08 S23-S21
      8  931  578 (S16E17 L=308)   515    108        1 ( 895/ 538- 975/ 631) -90.0 E18-E15 S19-S12
      9  395  518 (S08E57 L=268)   605    102        4 ( 357/ 488- 435/ 551) -38.7 E59-E53 S10-S06
     10  505  297 (S26E61 L=264)   960    225        2 ( 488/ 256- 529/ 339)  35.0 E66-E56 S29-S22
  

    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