Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1726 (N37E28 L=300)  1061    225        3 ( 920/1625-1018/1809) -31.4 E36-E18 N31-N42
      2  662 1680 (N26E45 L=282)   855    184        4 ( 620/1608- 709/1750) -22.4 E54-E37 N21-N29
      3 1521 1697 (N57W20 L=347)  1496    204        2 (1440/1653-1591/1729)  16.7 W10-W30 N50-N61
      4  169  980 (S27E65 L=262)   853    126        1 ( 142/ 916- 202/1032)  45.0 E70-E60 S30-S22
      5  539  746 (S33E24 L=303)   631    197        1 ( 448/ 718- 647/ 766)  23.6 E32-E16 S35-S28
      6 1053  594 (S28W14 L=341)   581    141        0 (1023/ 526-1088/ 653) -77.5 W12-W15 S32-S23
      7 1420  527 (S21W40 L=  7)   322     76        1 (1393/ 491-1447/ 560) -78.7 W38-W41 S23-S18
      8 1271  450 (S30W35 L=  2)   329    139        0 (1233/ 420-1323/ 506) -11.3 W29-W39 S31-S27
      9  476  398 (S61E32 L=295)  1768    241        5 ( 393/ 334- 550/ 471) -24.0 E48-E21 S65-S56
     10  891  217 (S61W28 L=355)   535     60        0 ( 867/ 185- 920/ 244) -13.0 W22-W35 S63-S58
  

    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