Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  736 1841 (N36E51 L=168)   564    176        0 ( 717/1806- 757/1879) -33.7 E56-E47 N33-N39
      2  379 1137 (S16E43 L=176)   529    156        0 ( 327/1076- 437/1179)  49.4 E46-E39 S20-S12
      3 1715 1085 (N18W41 L=260)   324     63        0 (1692/1045-1734/1121) -36.9 W38-W42 N16-N19
      4 1597  979 (N08W34 L=254)   879    178        1 (1547/ 902-1659/1062) -28.6 W29-W40 N05-N11
      5  508  984 (S22E30 L=189)   263    108        0 ( 480/ 956- 552/1036) -15.9 E33-E26 S22-S19
      6 1904  906 (N15W65 L=284)  2719    275        1 (1883/ 788-1928/1024) -42.7 W58-W73 N09-N21
      7 1087  922 (S10W06 L=226)   276    112        0 (1075/ 868-1100/ 993) -56.3 W03-W08 S12-S06
      8  730  905 (S21E14 L=206)   235     87        1 ( 685/ 891- 772/ 922)  11.3 E16-E11 S21-S19
      9  355  734 (S41E39 L=180)   898    290        2 ( 269/ 705- 460/ 757)  17.5 E48-E29 S44-S36
     10 1056  599 (S29W15 L=235)   317    134        0 ( 986/ 587-1123/ 619)  35.5 W11-W18 S30-S25
  

    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