Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  870 1873 (N56E26 L=287)   908    102        2 ( 833/1832- 924/1897)  60.3 E30-E21 N50-N59
      2  906 1702 (N39E14 L=299)  1067    187        3 ( 849/1613- 945/1765) -50.2 E20-E09 N32-N44
      3  620 1564 (N27E32 L=281)   528    120        1 ( 567/1517- 657/1591) -24.0 E37-E28 N24-N29
      4 1443 1450 (N25W25 L=338)   247     26        0 (1427/1428-1461/1472)  56.3 W24-W26 N23-N26
      5  663  832 (S18E22 L=291)   381    103        0 ( 646/ 781- 692/ 879)  63.4 E23-E19 S20-S14
      6  455  763 (S23E38 L=275)   966    216        2 ( 380/ 689- 548/ 816)  25.0 E45-E30 S27-S19
      7  828  733 (S23E11 L=302)   246     79        0 ( 789/ 721- 867/ 747)   0.0 E13-E08 S24-S22
      8  796  494 (S39E14 L=299)   830    114        2 ( 745/ 463- 853/ 543)  14.0 E17-E09 S41-S35
  

    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